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COSMOGENIC  ”P  AND  ”P  IN  TOE  ATMOSPHERE 
AND  OUGOTROPHIC  OCEAN  AND  APPUCATIONS 
TO  THE  STUDY  OF  PHOSPHORUS  CYCLING 

by 

Nathalie  Waser 
ABSTRACT 

Cosmic-ray  produced  (t  =  14.28  days)  and  ”p  (t  =  25.3  days)  are  power¬ 
ful  tracers  of  upper  ocean  P  cycling,^when  coupled  with  time-series  estimates  of  the 

32  33 

atmospheric  sources.  The  cycle  of  cosnoogenic  Pand  P  in  the  atmosphere  and  euphoric 
zone  of  the  ocean  was  investigated  in  the  Sargasso  Sea.  A  method  was  developed  for  the 
determination  of  the  low-level  beta  activities  of  both  nuclides  in  rainwater  and  marine 

32  33 

particulate  matter.  The  fallout  rates  of  Pand  P  by  precipitation  were  determined  over  a 
period  of  12  monriis  at  a  marine  site,  at  Bermuda,  coinciding  with  measurements  of  the 
activities  and  activity  ratio  ^^p/^P  in  suspended  particulate  matter  and  plankton  tows  at  the 
Bermuda  Atlantic  Hme-series  Study  (BATS)  station  (31°  SO*  N,  64°  10*  W).  The  in 
situ  production  rates  of  ^  and  ^^P  in  the  upper  ocean  were  estimated  by  measuring  the 
activities  of  ^  and  ^^P  produced  in  Cl,  K  and  S  targets  exposed  to  cosmic  rays  at  Woods 
Hole,  on  Mount  Washington,  and  on  TAiguille  du  Midi.  Fallout  by  precipitation  was 
found  to  be  the  major  source  of  radiophosphorus  to  the  ocean;  in  situ  production  in  the 
upper  ocean  was  found  to  account  for  5  %  for  P  and  1  %  fOT  P  of  the  wet  deposition 
rates.  Knowledge  of  all  the  sources  of  radiophosphorus  to  the  Sargasso  Sea  allowed  the 

32  33 

cycling  of  P  and  P  in  suspended  particulate  matter  and  macrozooplankton  to  be  studied. 
The  study  was  based  on  a  new  approach  consisting  of  the  determination  of  the  activity 
ratio  ^  V/  P  in  the  different  pools  of  phosphorus.  The  activity  ratios  ^*p/^P  in  suspended 
particulate  matter  and  plankton  tows  ranged  from  1.0  to  5.0  and  were  either  equal  to  or 
higher  than  the  activity  ratio  in  rainwater.  The  noonthly  activity  ratio  in  rain  was 
remaritably  ctmstant  and  on  average  0.96.  The  increase  in  the  ratio  in  plankton  relative  to 
rain  can  be  interpreted  as  aging  of  the  phosphorus  in  that  pool  or  in  the  ^wl  from  which  it 
is  derived.  This  is  because  of  differential  radioactive  decay  rates  of  P  and  P.  The 
activity  ratio  was  higher  in  particle  collections  dominated  by  higher  levels  in  die  ocean 
food  web.  The  residence  time  of  P  in  macrozooplankton  was  estimated  to  range  from  40 
to  60  days.  A  grazing  rate  of  macrozooplankton  of  0.024  day  was  determined,  allowing 
the  determination  of  a  grazing  flux  of  110  mg  C/m^  /d  and  a  flux  of  particulate  carbon 
due  to  defecation  of  zooplankton  of  30  mg  C/m^  /d. 
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and  net  count  rates  at  the  time  of  collection  (N3  2  and  N3  3 )  in  plankton 

tows  and  suspended  particulate  samples  (Millipore  3-pm)  samples  collected 

during  BATS  41A,  24-25  February  1992  152 

Table  A.9.  Thickness  (d),  yield  (Y),  counting  efficiencies  (£3  2  and  £33) 

and  net  count  rates  at  the  'ime  of  collection  (N3  2  and  N3  3 )  in  suspended 

particles  and  in  plankton  tows  collected  during  £N  235  in  April  1992  at 

BATS  station  153 

Table  A.10.  Thickness  (d),  yield  (Y),  counting  efficiencies  (£3  2  and  £33) 

and  net  count  rates  at  the  time  of  collection  (N3  2  and  N3  3 )  of  the  samples 

exposed  to  cosmic  rays  in  Woods  Hole  53°  N),  Mt.  Washington 

(X»  55.7®  N)  and  I’Aiguilledu  Midi  (X=  47®  N)  153 
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INTRODUCTION 


In  order  to  understand  the  magnitude  and  variations  of  primary  productivity  in 
the  oligotrc^c  Sargasso  Sea,  as  in  any  other  oceanic  province,  it  is  essential  to  study  the 
factCMTS  that  control  productivity.  Primary  productivity  in  the  Sargasso  Sea  is  likely  regu¬ 
lated  by  the  supply  of  nutrients  into  the  euphotic  zone  during  deep  winter  mixing  (Menzel 
and  Ryther,  1961;  Marra  et  al.,  1990)  and  by  the  rate  of  nutrient  recycling  (Goldman, 
1988).  Time-series  measurements  near  Bermuda  reveal  correlations  between  the  export 
flux  of  carbon,  the  flux  of  particle-reactive  radionuclides  in  the  deep  ocean  and  primary 
productivity  in  the  euphotic  zone  with  a  time-lag  consistent  with  the  sinking  speed  of 
large  particles  (Deuser  and  Ross,  1980;  Bacon  et  a/.,  1985;  Deuser,  1986;  Asper  et  al., 
1992).  The  study  of  nutrient  cycles  in  the  upper  ocean  remains  a  central  issue  for  a  com¬ 
plete  undorstanding  of  the  biological  pump  and  its  effect  on  the  deep  ocean. 

There  is  recent  evidence  that  primary  productivity  in  the  Sargasso  Sea  is  higher 
than  previously  estimated  with  *  *  C  incubation  techniques  (Jenkins  and  Goldman,  1985; 
Jenkins,  1988).  Recent  primary  productivity  estimates  (Michaels  et  al.,  1992)  at  the 
Bermuda  Atlantic  Time-series  Study  (BATS)  station  are  two-fold  higher  than  the  histori¬ 
cal  rate,  possibly  reflecting  the  use  of  cleaner  techniques  (Fitzwater  et  al.,  1982).  The 
fundamental  processes  that  govern  both  the  rate  of  carbon  fixation  and  its  export  to  the 
deep  ocean  remain  poorly  understood.  Recent  studies  of  photosynthesis-inadiance 
relationships  suggest  that  there  is  a  direct  physiological  control  by  nutrients  on  primary 
productivity  in  the  Ntnrth  Sargasso  Sea  (Platt  et  al.,  1992;  Falkowski  et  al.,  1992). 
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Overall,  the  exact  relationship  between  new  production,  prinuuy  production  and  nutrients, 
particularly  which  nutrient  is  limiting  productivity,  remain  in  doubt 

Anoong  the  processes  that  are  most  impcxtant  in  nutrient  cycles  and  their  rela¬ 
tion  to  ocean  productivity  are  uptake  and  regeneration  rates  and  thus  the  fate  and 
pathways  of  nutrients  from  the  autotrophs  to  the  higher  trophic  levels.  The  amount  of  C, 
N  and  P  exptxt  to  the  deep  ocean  depends  not  only  on  the  anoount  of  nutrients  being 
supplied  to  the  euphotic  zone,  but  also  on  the  efficiency  of  the  transfer  of  C,  N  and  P  to 
higher  trophic  levels  (Eppley  and  Peterson,  1979;  Goldman,  1988;  Michaels  and  Silver, 
1988).  Therefore,  it  is  impcntant  to  determine  the  uptake  rates,  regeneration  rates  and 
turnover  times  of  nutrients  in  the  euphotic  zone,  because  they  directly  determine  the 
extent  of  new  production  versus  primary  production  (Eppley  and  Peterson,  1979). 

Phosphorus  is  one  of  the  fundamental  building  blocks  of  all  forms  of  living 
aquatic  organisms.  Its  study  has  been  neglected  in  favor  of  the  cycles  of  nitrogen  and 
carbon.  It  is  a  widespread  view  among  biologists  that  N  is  the  limiting  nutrient  in  ocean 
ecosystems  as  opposed  to  P  limitation  in  freshwater  ecosystems  (Perry  and  Eppley, 
1981).  Recently  the  role  of  P  in  the  oligotrophic  ocean  has  received  renewed  interest 
(Fanning,  1989;  Michaels  et  al„  1992;  Karl  et  al.,  1992;  Fanning,  1992).  With  the  results 
of  the  time-series  studies  of  two  oligotrophic  oceans,  the  Sargasso  Sea  and  the  North 
Pacific  Sub-tropical  Gyre,  the  role  of  P  in  the  marine  environment  has  been  reexamined 
(Karl  et  al.,  1992;  Michaels  et  al„  1992).  Fanning  (1992)  recently  reported  puzzlingly 
high  ratios  of  dissolved  inm'ganic  nitrogen  to  dissolved  inorganic  phosphorus  (DIN/DIP) 
in  the  Nmth  Atlantic  Ocean  in  contrast  to  all  other  oceans,  raising  the  question  of  the 
relative  role  of  N  and  P  as  limiting  nutrients.  DIN/DIP  ratios  below  the  euphotic  zone 
down  to  800  m  at  BATS  are  high,  sometimes  up  to  40  (mole  ratio),  and  always  higher 
than  the  Redfield  ratio  of  16  (Knap  et  al.,  1991;  Knap  et  al.,  1992;  Michaels  eial.,  1992). 
Unfortunately,  dissolved  organic  nitrogen  (DON)  and  dissolved  organic  phosphorus 
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(DOP)  data  are  not  measured  routinely  at  BATS,  so  total  dissolved  nitrogen  (TDN)  and 
total  dissolved  phosphorus  (TDP)  cannot  be  compared.  The  high  DIN/DIP  regenouted 
below  the  euphodc  zone  might  be  indicative  of  P-depleted  particles  being  regenerated  or 
another  source  of  nitrogen.  These  results  are  voy  puzzling  and  might  be  explained  by  (1) 
nitiogen  fixation  (Karl  et  al.,  1992;  Fanning.  1992),  (2)  high  atmo^heric  supply  of 
nitrogen  relative  to  phosphorus  due  to  polluted  air  (Fanning,  1989),  or  (3)  blooms 
dominated  by  coccolithophorids  (Fanning,  1992).  These  observations  have  pointed  out 
the  weakness  of  the  current  understanding  of  P  and  N  cycles  and  their  coupling  in  the 
oligotrophic  ocean.  If  nitrogen  fixation  is  truly  an  important  process  supplying  N  to 
phytt^lankton  in  the  ocean,  then  P  should  be  the  ultimate  limiting  nutrient 

The  subject  of  this  thesis  is  the  study  of  the  phosphorus  cycle  in  the  Sargasso 

32 

Sea  using  two  naturally  produced  radioisotopes  of  phosphorus.  Cosnnogenic  P  (t^^^= 
14.28  d;  E  =  1.71  MeV)  and  ”p  (t  =  25.3  d;  E  =  0.249  MeV)  arc  beta  er  .ters 

iMX  1/2  max 

produced  by  spallation  reactions  in  Ar,  in  the  atmosphere,  and  by  spallation  in  Q,  S,  Ca 

32 

and  K  in  the  ocean.  Very  recently  the  first  measurements  of  the  minute  activities  of  P 

33 

and  P  in  seawater  were  obtained  (Lai  and  Lee,  1988;  Lai  et  al.,  1988).  More  extensive 
measurements  of  P  were  made  in  TDP,  DIP,  DOP  and  plankton  tows  (Lai  et  al.,  1988; 

32 

Lee  et  al.,  1991;  Lee  et  al.,  1992)  demonstrating  the  feasibility  of  applying  P  and  per- 

33  33 

haps  P  to  the  study  of  the  P  cycle  in  the  upper  ocean.  There  were  very  few  data  on  P 

33 

because  P  is  a  weak  beta  emitter  and  thus  hard  to  detect 

32  33  33  32 

The  approach  of  this  thesis  consists  of  utilizing  P  and  P  and  the  ratio  P/  P 
to  determine  the  residence  time  of  P  in  plankton  and  in  the  euphotic  zone  and  the  export 
of  P  out  of  the  euphotic  zone  at  the  Bermuda  Atlantic  Time-series  Study  (BATS)  station 
(31°  50*  N,  64°  10'  W)  in  different  seasons.  In  order  to  determine  the  turnover  times  of 

32  33 

P  in  particles  and  in  the  total  dissolved  phosphorus,  the  mass  balances  of  P  and  P  in  the 

32  33 

euphotic  zone  need  to  be  accurately  constrained.  The  fallout  rates  of  P  and  P  by 
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pndpitation  at  Bermuda  and  the  in  situ  production  rates  in  seawater  have  thus  been 

32  33 

determined.  The  wet  deposition  rates  of  Pand  P  have  not  been  previously  estimated  at 
Bermuda,  and  it  is  impmtant  to  document  the  variability  of  the  supply  of  these  short-lived 

32  33 

nuclides  at  the  study  site.  The  in  situ  production  of  P  and  P  has  been  estimated  pre¬ 
viously  (Lai  et  1988)  but  not  measured  directly  in  all  of  the  possible  targets  in  sea 
salt 
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Chapter  One 

METHODS  OF  MEASUREMENT  OF  AND 


I  INTRODUCTION 

32  33 

Activities  of  P  (t  =  14.28  days)  and  P  (t  =  25.3  days)  in  seawater  have  only 

U2  1/2 

been  measured  very  recently  (Lai  et  al.,  1988;  Lai  and  Lee,  1988;  Lee  et  al.,  1991;  Lee  et 
al.,  1992).  In  spite  of  dieir  tremendous  potential  for  tracing  processes  involving  phos¬ 
phorus  in  the  environment,  the  very  low  activities  in  nature  have  impeded  their  use  in 
biogeochemistry.  Very  low  counting  backgrounds  and,  more  importantly,  high  specific 
activities  are  essential  to  their  detection  by  low-level  ^  counting.  The  background  can  be 
reduced  by  the  use  of  an  anticoincidence  system  surrounding  the  primary  detector.  High 
activities  can  be  obtained  by  preconcentration  in  the  field.  Preconcentration  of  the  ac¬ 
tivity  in  rainwater  is  achieved  here  by  concentrating  the  activity  of  the  dissolved  phase  on 
alumina,  a  specific  adsorbent  of  phosphate,  in  a  system  similar  in  principle  to  the  one 
used  by  Silker  (1969,  1971).  Large  volume  filtration  is  utilized  to  collect  suspended 
particulate  matter.  To  separate  ^  ^  P  and  ^  ^  P  a  nuclear  technique  has  to  be  sought  The 
two  isotqies  are  distinguished  by  their  half-lives  and  by  the  very  distinct  energy  spectrum 
of  the  electrcHis  emitted  in  the  decay.  Their  half-lives  are  not  diffnent  enough  to  enable 
their  sqiaratitNi  by  decay-curve  analysis.  In  contrast  ^  ^  P  (£  =  1.71  MeV)  is  a  much 

max 

m(»e  energetic  emitter  than  ^  ^  P  (E  =  0.249  MeV).  As  a  consequence,  the  absorption 

mix 

characteristics  of  the  two  isotopes  differ  enough  that  an  external  absorber  of  a  given 
thickness  will  absorb  all  the  radiation  from  ^  ^  P  but  only  a  small  fraction  of  the  radiation 
from  ^  *  P.  I  have  used  this  absorber  method  to  separate  ^  ^  P  and  ^  ^  P. 
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n  COLLECTION  AND  EXTRACTION  OF  RADIOPHOSPHORUS 


!•  Principle 

The  preconcentiation  of  dissolved  ^^P  and  ’^P  from  rainwater  has  been 
adiieved  by  extracting  the  radiophosphcxus  cm  alumina.  The  smption  pr(^>eities  of 
alumina  with  respect  to  orthc^hosphate  have  been  well  studied.  Sorption  studies  have 
been  carried  out  at  different  pHs,  PO4  ccmcentration,  and  in  the  presence  of  inorganic  and 
organic  species  (Chen  etal.,  1973;  Huang,  1975;  Stumm  et  al.,  1980;  Morel,  1983).  The 
process  by  which  PO4  binds  to  the  surface  has  been  described  as  a  ligand  exchange 
reaction  in  which  PO4  rqrlaces  a  hydroxide  group  at  the  surface  (Stumm  et  al.  1980).  In 
additicm,  my  own  experiments  have  shown  faster  removal  of  PO4  from  solution  cm 
freshly  precipitated  aluminum  hydroxides  than  on  freshly  precipitated  ferric  hydroxides 
in  the  absence  of  other  incn-ganic  ac  organic  species  at  pH  7.  For  this  reason,  alumina  was 
the  adsorbent  of  choice  for  this  study. 

2-  Rain  samiriing 

Rain  samples  were  collected  from  the  roof  of  a  building  at  the  Bermuda 
Biological  Station  for  Research  (BBSR)  and  from  a  collector  situated  on  the  roof  of  Qaik 
building  in  Wcxxls  Hole.  The  collectors  enabled  the  collection  of  the  10  to  60  t  of  rain¬ 
water  necessary  to  detect  radiophosphorus.  The  rainwater  was  spiked  with  100  pmoles  of 
stable  KH2  PO4 ,  which  served  as  carrier  and  enabled  the  yield  of  the  extraction  and 
purification  procedures  to  be  determined.  The  rainwater  was  left  to  equilibrate  fcH*  4-6 
hours.  The  sample  was  then  pumped  (with  a  peristaltic  pump)  through  a  plexiglass  unit 
containing  about  7  g  of  alumina  (Activated  alumina,  98%  powder,  Matheson  Coleman 
and  Dell)  packed  between  two  porous  polyethylene  filters  (3.5  cm  in  diameter).  The 
extraction  efficiency  of  PO4  was  better  than  98%  for  a  volume  of  60  1  of  rainwater  and 
flow  rates  of  5  1/min  or  less.  This  efficiency  was  confirmed  by  laboratory  experiments. 
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In  the  labcM'atory,  SO-140  I  of  distilled  water  and  filtered  seawater  (0.2  ttm)  fitnn 
Vineyard  Sound  were  spiked  with  stable  phosphate.  Phosphate  was  extracted  on  alumina 
with  efficiencies  higher  than  98%  in  the  case  of  distilled  water.  In  the  experiment  with 
filtered  seawater,  the  efficiencies  decreased  from  1(X)  to  90  %  during  the  course  of  the 
extraction.  Once  the  extraction  was  completed,  the  resin  was  sent  by  Federal  Express  to 
WHOI,  where  the  alumina  was  processed  in  the  lab(»ah>ry. 

3-  Particulate  sampling 

Different  size  fractions  of  particles  were  sampled.  Particles  were  collected  using 

TM 

293-mm  diameter  Millipore  filters  (3  or  8  pm  mesh  size).  These  samples  were  collected 
with  a  simple  filtration  system  which  consisted  of  either  an  impeller  pump  or  a  diaphragm 
pump  connected  to  a  293-mm  filter  holder.  The  pump  was  placed  upstream  of  the  filter 
holder.  For  surface  samples  a  vacuum  hose  was  lowered  to  a  depth  of  3-5  meters  over  the 
side  of  the  ship.  For  deep  samples  the  vacuum  hose  was  attached  to  a  weight  on  the 
hydrowire,  and  the  rest  of  the  hose  was  allowed  to  float  freely  at  the  surface.  With  this 
system  I  could  pump  1,(XX)  i  at  10  to  45  psi  in  about  1  hour.  The  filter  clogged  fast  and 
was  changed  as  soon  as  the  flow  rate  started  to  decrease,  which  was  about  every  hour. 
The  filters  were  immediately  frozen  to  avoid  release  of  phosphorus  compounds  (Collier 
and  Edmond,  1984).  Particles  larger  than  67,  150,  300  and  500  pm  were  harvested  with 
plankton  nets.  The  particles  were  sampled  at  the  surface  and  at  depth.  Samples  were 
recovered  from  the  cod  end  and  gravity  filtered  on  Whatman  filters  and  immediately 
frozen. 

m  RADIOCHEMICAL  PURIFICATION 

!•  Rain  samples 
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The  AU  Os  was  removed  £rom  the  filter  holders  and  slurried  in  1  N  NaOR 
NaOH  has  been  reported  to  be  efficient  at  releasing  PO4  adsorbed  by  aluminum  oxides 
and  hydroxides  (Dickman  and  Bray,  1940;  Turner  and  Rice,  1954).  NaOH  (IN)  removed 
adsorbed  P  from  AI2  O3 ,  but  large  volumes  of  it  were  necessary  to  get  good  recoveries  (> 
50%).  Typically  for  rain  samples  800  mi  to  1  1  of  IN  NaOH  was  used  and  boiled  down 
to  about  2(X)  mi .  The  alkaline  solution  was  then  poured  gradually  into  a  6  N  HNO3 
solution.  The  pH  of  the  solution  was  kept  below  2  to  avoid  formation  of  massive  amounts 
of  aluminum  hydroxide  colloids. 

2-  Particiilate  samples 

TM 

The  frozen  Millipore  filters  and  the  zooplankton  samples  were  refluxed  for  a  few 
hours  in  concentrated  HNO3  until  the  solution  was  clear.  If  refractory  particles  were 
present,  the  solution  was  filtered.  The  solution  was  then  taken  down  to  a  small  volume 
and  diluted  to  about  1  N  HNOa  with  distilled  water.  At  this  point,  an  aliquot  of  the 
solution  was  analysed  for  phosphate,  i.e.  the  yield  monitor,  by  the  classic  molybdenum 
blue  method  (Murphy  and  Riley,  1962;  Koroleff,  1983). 

3*  Radiochemical  procedure 

The  chemical  purifrcation  of  phosphorus  is  based  on  a  series  of  specific  phosphate 
precipitations,  which  have  been  described  in  the  classic  chenustry  and  radiopurification 
methods  (Hillebrand  e/  a/.,  1959;  Mullins  and  Leddicotte,  1962;  Volchock  and  De 
Planque,  1983;  Whaley  and  Ferrara,  1972).  The  sample  solution  in  HNO3  was  heated, 
and  an  excess  of  amnxrnium  molybdate  was  added.  This  solution  was  heated  to  about  30- 
40  and  stirred  until  a  yellow  precipitate  of  ammonium  phosphomolybdate 
(NH4  )3  PO4  (M0O3  )i  2  appeared.  The  amnxrnium  phosphomolybdate  precipitation  is  a 
preliminary  separation  which  reduces  the  amount  of  heavy  metals,  particularly  Fe,  Co, 
Ni,  Cr,  Ti  and  Zr  (Mullins  and  Leddicotte,  1962).  Heating  of  the  solution  was  stopped 
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while  stirring  was  maintained  and  the  precipitation  process  allowed  to  continue.  If  heat¬ 
ing  was  continued  too  long  (at  temperatures  higher  than  SO  °  C),  molybdic  acid  would 
^ait  to  precipitate,  and  the  precipitate  would  be  contaminated  with  Si.  As  and  V 
(Hillebrand  et  al.,  1959).  After  about  one  hour  the  precipitate  was  vacuum  filtered 

TM 

through  a  Millipore  HA  filter  (0.45  pm,  47  mm  in  diameter)  and  washed  with  1  N 
HNO3 .  Sometimes  the  precipitate  was  so  fine  that  a  0.2  pm  Millipore  filter  had  to  be 
used.  The  precipitate  was  dissolved  with  NH4  OH  and  the  solution  acidified  to  1  N 
HNO3 .  Ammonium  phosphomolybdate  was  precipitated  a  second  time  and  redissolved 
with  anunonia.  The  pH  of  the  solution  was  then  lowered  to  about  7  with  HQ.  The  solu¬ 
tion  was  cooled  in  an  ice  bath.  A  reagent  containing  MgOa  and  NH4  Cl  was  added  to 
the  cold  solution.  Drops  of  NH4  OH  were  added  while  stirring.  A  white  crystalline 
precipitate  of  ammonium  magnesium  phosphate  formed.  An  excess  of  concentrated 
NH4  OH  was  added.  The  NH4  MgP04  .6H2  O  precipitate  in  solution  was  stirred  for  15 
min  and  cooled  for  half  an  hour  in  an  ice  bath.  This  step  is  not  a  separation,  since  there  is 
a  large  number  of  interferences  with  the  precipitation.  The  advantage  of  the 
NH4  MgP04 .6H2  O  precipitate  is  that  it  can  be  dissolved  in  HCl  while 
(NH4  h  PO4  (M0O3  )i  2  cannot.  The  NH4  MgP04  .6H2  O  precipitate  was  vacuum  fil¬ 
tered  through  a  Millipore  HA  filter  (0.45  pm,  47  mm  in  diameter),  washed  with  dilute 
NH4  OH,  and  dissolved  in  9  N  HQ.  A  cation  exchange  column  of  AG-50W-X8  Q,  100- 
200  (Analytical  grade  resin  from  Biorad  Laboratory)  mesh  was  conditioned  with  3 
volumes  of  9  N  HQ.  The  solution  was  loaded  on  the  column  and  1  volume  of  9  N  HQ 
was  passed  through  to  rinse.  The  cation  exchange  resin  allows  separation  of  Ca,  K,  Fe, 
Al,  V,  W,  Zn,  Zr  and  Ti,  which  will  interfere  with  the  final  MgNH4  PO4  .6H2  O  precipita¬ 
tion  (Hillebrand  et  al.,  1959;  Mullins  and  Leddicotte,  1962;  Whaley  and  Ferrara,  1972). 
Concentrated  NH4  OH  was  added  to  make  the  pH  neutral,  and  the  solution  was  put  again 
in  an  ice  bath  prior  to  the  last  precipitation.  MgNH4  PO4  .6H2  O  was  precipitated  as 
described  previously  and  vacuum  filtered  on  a  preweighed  Millipore  HA  filter  (0.45  pm, 
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25  mm  in  diameter).  The  area  of  the  precipitate  was  2.7S  cm^ .  The  precipitate  was 
hygroscopic  and  it  was  thus  dried  carefully  until  the  weight  was  constant.  It  was  then 
mounted  face  down  on  a  film  of  Mylar  (thickness  =  0.9  mg/cm^  ).  The  film  of  Mylar 
was  supported  by  a  white  Delrin  ring.  A  silver  planchet  (140  mg/cm^ )  was  placed  on  top 
of  the  ^ter,  and  the  source  was  sealed  with  tape  to  ensure  no  change  in  the  weight  of  the 
precipitate.  The  source  was  counted  on  an  anticoincidence  low-level  ^  counter. 

IV  LOW-LEVEL  BETA  COUNTING 

A  thin-wall  anticoincidence  counter  was  used  for  low-level  beta  counting.  The 
basic  feature  of  low-level  ^  counting  that  distinguishes  it  from  ordinary  ^  counting  is  the 
use  of  an  anticoincidence  guard  detector  next  to  the  primary  detector.  The  primary  detec¬ 
tor  used  here  was  similar  in  design  to  that  of  Lai  and  Schink  (1960).  The  detector  is  a  gas- 
flow  counter  and  operates  in  the  Geiger-Muller  region.  The  plateau  starts  at  820-860 
Volts.  The  detector  was  operated  at  a  voltage  centered  on  the  plateau.  Q-gas  (98.7% 
Helium  and  1.3%  isobutane)  flows  continously  through  the  sealed  counter  chamber.  The 
2.54  cm  window  is  a  thin  0.9  mg/cm^  aluminized  Mylar.  The  background  of  the  detector 
is  significantly  reduced  by  the  guard  detectors.  Typically  the  background  is  about  0.25- 
0.5  cpm  with  an  error  of  0.02  cpm.  This  background  is  about  50  times  lower  than  that  of 
a  detector  without  the  guard  detectors.  The  counting  system  is  usually  stable  (no  change 
in  the  background)  for  periods  of  30  days  or  more.  The  detector  measures  beta  radiation 
in  pulse-type  mode.  The  pulses  from  the  source  and  the  guard  detectors  are  first 
amplified  by  Pelagic  Electronics  Model  7030-2  preamplifiers.  They  are  further  amplified 
by  a  Canberra  Model  2012  amplifier.  Pulses  from  the  guard  detector  are  sent  to  a  Pelagic 
Electronics  Model  7030-4  gate  driver.  Pulses  from  the  source  are  sent  to  a  Canberra 
Model  2032  dual  discriminator,  which  is  coupled  to  a  Gate  Driver  and  further  through  a 
Canberra  Model  2035 A  single-channel  analyzer  and  a  Canberra  Model  1476A  scaler. 
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V  SEPARATION  OF  AND 


A  nuclear  method  was  used  to  separate  the  two  radioisotopes  of  P  based  on  the 
difference  in  energy  of  the  beta  particles  emitted  in  their  decay.  Phosphorus-32  is  a  hard 
P  emitter  with  an  end-point  energy,  E  =  1.71  McV,  while  ”  P  is  a  soft  P  emitter  with 

DUUC 

E  =  0.249  MeV.  This  difference  in  energy  between  the  two  isotopes  results  in  different 

max 

absorption  characteristics  for  the  ^  particles.  The  interaction  of  p  particles  with  matter  is 
described  by  the  experimental  half-thickness  of  absorption  which  is  the  thickness  of  a 
solid  material  between  the  source  and  the  detector  that  will  absorb  half  of  the  emitted 
radiation.  The  intensity  of  a  source  I  as  a  function  of  the  external  absorber  thickness  x  is 
simply  approximated  by:  I  =  lo  exp(  -  In  2  x  where  lo  is  the  intensity  of  the  source 
with  no  external  abstnber  (Libby,  1947).  I  have  found  that  the  activity  of  *  ^  P  is 
decreased  by  a  factor  of  2  when  an  aluminum  foil  of  71.6  ±  2  mg/cm^  is  placed  between 
the  source  and  the  detector  (Figure  1.1).  The  mass  absorption  coefficient  p  (p  =  In  2  /x^^^) 
was  thus  calculated  to  be  9.7  +  0.3  cm^  /g.  These  experimental  values  arc  comparable  to 
some  previous  studies.  An  empiri';al  formula  relates  p  and  the  maximum  p  energy  as 
follows:  p  =  22/E  for  0.5  s  E  i  6  MeV.  For  E  =1.71  MeV  the  calculated  p  is 
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10.8  cm^  /g  (Price  et  al.,  1958)  which  is  very  close  to  the  value  of  9.7  cm^  /g  found  in 
this  study.  Libby  (1957)  reported  a  different  relationship,  i.e.  1/p  =  55  E*^,  which  yields 
a  lower  p  and  a  higher  x^^^  of  84  mg/cm^ .  For  the  range  of  absorber  used  in  my  separa¬ 
tion  method,  i.e.  6  to  40  mg/cm^  ,  the  two  extreme  half-thicknesses  72  and  84  mg/cm^ 
yield  activities  which  differ  by  1  %  to  5  %.  In  the  case  of  ’  ’  P,  a  foil  of  only  4.3  ±  0. 1 
mg/cm^  blocks  50%  of  the  P  radiation  (Figure  1.2). 
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Absorber  thickness  (ing/cin2) 


Figure  1.1.  Net  count  rate  of  P  as  a  function  of  external  absorber  thirimeg^ 


This  half-thickness  is  in  close  agreement  with  of  4.9  mg/cm^  found  for  ^  ’  Ca 
(E  =  0.2SS  MeV)  in  A1  (Libby,  1957).  Thus  from  the  absorption  curves  of  ’  ^  P  and  ’  *  P 
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one  can  choose  a  thickness  of  aluminum  that  will  block  essentially  all  the  radiation  fhxn 
^  ^  P  and  decrease  the  radiation  from  ^  ^  P  by  a  known  small  anraunt.  The  thickness  of 
aluminum  was  chosen  so  that  the  estimated  activity  of  ^  ^  P  in  die  source  would  be 
decreased  to  less  than  0.02  cpm  when  the  source  was  counted  with  die  absorber. 

In  (nrder  to  separate  ^  ^  P  and  ^  ^  P  the  samples  were  counted  repeatedly  in  time 
with  and  without  external  absorber.  The  time-dependent  curve  generated  without  absenrb- 
er  represents  the  total  activity  of  the  sample.  The  curve  with  absorber  represents  the 
activity  of  ^  ^  P,  decreased  slighdy  due  to  the  absorption  of  P  particles  in  the  absorbing 
material.  The  initial  net  count  rates  ^  ^  N®  and  ^  ^  N®  of  ^  ^  P  and  ’  ’  P,  i.e.,  the  net  count 
rates  at  the  time  the  sample  was  first  counted,  are  then  determined  by  the  difference  of  the 
two  time-dependent  curves.  The  activities  at  the  time  th"  iample  was  collected  were 
calculated  from  the  known  decay  rates  of  ^  ^  P  and  ^  ^  F  The  equations  describing  the 
total  gross  count  rates  without  absorber  and  with  absorber  at  a  given  time  are  the 
following : 

(1.1)  G^  =  N^  +  B  =  ^*N  +  ”N  +  B  =  N®  cxp(-A.t)  +  B 

(1.2)  G  =N  +B  =  ^^N  ”p  +  B  =  ^^N®  exp(-’H  t)  +  B 

&  ft 

where 

N^:  total  net  count  rate  in  cpm 
N^:  net  count  rate  with  absorber  in  cpm 
B  :  background  count  rate 
^  ^  N:  net  count  rate  of  ’  ^  P  in  cpm 
’  ^  N:  net  count  rate  of  ^  ^  P  in  cpm 
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’  ^  P:  fraction  of  ’  ^  P  tnuismitted  diiou^  the  abswber. 

’*P«cxp(-In2x/71.6). 
x:  thickness  of  the  abs(»ber,  in  rng/cm^ . 

’  ’  fractionof  ’  ’  P  transmitted  thitMigh  the  absofber. 

”p-cxp(-ln2x/4.3). 

’  ^  X:  decay  constant  of  ^  ^  P. 

In  practice  each  sample  was  counted  long  enough  to  obtain  8%  or  less  counting 
error  both  with  the  external  absorber  and  without.  Typically  the  thickness  of  the  absorber 
ranged  frtMn  6  to  40  mg/cm^  for  the  low  count  rate  (plankton  samples)  to  the  high  count 
rate  samples  (rainwater  samples),  respectively.  N  (i.e.,  +  was 

ft 

usually  i^proximated  by  ^  ^  N  ^  ^  since  the  absorber  was  chosen  to  allow  absorbtitm  of 
essentially  all  the  radiation  of  ^  ^  P  so  that  ’  ^  N  ^  p  <  0.02  cpm.  This  approximati(»i 
was  verified  when  ’  *  N®  and  ^  ’  N®  were  determined,  by  calculating  ®  ’  N®  ”  p.  The 
two  curves  N  (t)  and  N  (t)  were  then  fitted  by  least  squares.  The  source  was  counted 

t  ft 

repeatedly  over  30  to  60  days,  allowing  a  minimum  of  10  points  on  and  5  on  N^.  The 
genoral  equations  used  for  the  fit  are: 

=  a  exp(-  a  t) 

N^  =  b  exp(-’2Xt) 

The  values  of  the  initial  net  count  rates  ^  ^  N®  and  ^  ^  N®  were  determined  by  the  inter- 
cq)ts,  a  and  b,  of  the  two  curves  with  y  axes.  At  this  point  the  term  ^  ^  N®  ^  ^  p  was 
computed  to  verify  that  ^  ^  N®  ^  ^  P  was  smaller  dian  0.02  cpm.  For  samples  of  ex¬ 
tremely  low  total  activity,  there  was  sometimes  a  remaining  activity  of  ^  ^  P  in  the 
absoiber  curve  (^  ^  N  ^  ^  p  >  0.02  cpm)  which  was  corrected  for  by  substracting  the 
estimated  activity  of  ^  ^  P  from  the  curve  with  absorber.  The  correction  was  iterated  until 
the  calculated  values  ^  ^  N®  and  ^  ^  N®  converged  to  constant  values. 
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VI  DETERMINATION  OF  ABSOLUTE  ACTIVmES 


The  initial  net  count  rates,  ^  ^  and  ^  ^  ,  were  converted  into  absolute  ac¬ 

tivities  ’  ^  A**  and  ^  ^  A*’  by  taking  into  account  the  counting  efficiency,  e,  and  the  yield 
ci  the  procedure,  y.  The  absolute  activities  were  calculated  as  follows : 

A®  =N®  /y  e 

The  absolute  activity  at  the  time  of  collection  is: 

A^®  »A®  exp(X  t) 

where,  t  is  the  time  elapsed  between  collection  and  initial  counting  time,  and  X  is  the 
decay  constant  The  activity  ratio  r  of’^Pto’^Pis: 
r  =  ”A®/’^A®=^’N® 

The  counting  efficiency,  e,  for  a  given  radioisotope  varies  with  the  thickness  of  the 
source  according  to  the  following  relationship  (Libby,  1947): 

(1,3)  e  /  eo  =  (1  -  exp(-  p  d))  /  (p  d) 

where 

p:  self-absorption  coefficient  in  cm^  /mg 
d:  thickness  of  the  source  in  mg/cm^ 
e:  counting  efficiency  for  a  source  of  thickness  d 
Co  :  counting  efficiency  for  an  infinitely  thin  source 
Typically,  the  thicknesses  of  the  samples,  d,  varied  from  2  to  12  mg/cm^ .  Self- 
absoiption  and  scattering  processes  in  the  source,  absoiption  in  the  air  and  window  and 
backscattering  processes  on  the  backing  material  were  all  included  in  one  parameter:  eo , 
the  counting  efficiency  for  an  infinitely  thin  source.  The  thickness  of  the  backing 
material  is  chosen  to  be  the  thickness  of  saturation  (or  close  to  it)  to  allow  maximum 
backscattering.  Experiments  have  shown  that  saturation  is  reached  for  a  thickness  of 
about  20%  of  the  range  of  a  nuclide  (Price  et  al.,  1958;  Wang  et  aL,  1975;  Choppin  and 
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Rydberg,  1980).  The  ranges  in  aluminum  are  800  mg/cm^  fen*  ^  ^  P  and  60  mg/cm^  f(»- 
^  ^  P,  and  the  saturation  thicknesses  are  12  mg/cm^  for  ^  ^  P  and  160  mg/cm^  for  ^  ^  P 
(Choppin  and  Rydberg,  1980).  The  fraction  of  the  elections  backscatteied  increases  with 
the  atomic  number  of  the  backing  material.  Silver  (A  =  47)  was  chosen  and  gives  a  facun- 
of  1.6S  increase  in  the  counting  efficiency  at  saturation  backscattering  for  ^  ^  P  (Choppin 
and  Rydberg,  1980).  For  ^  ^  P,  absorption  in  air  and  in  the  window  of  the  detector 
eliminates  scmie  of  the  backscatteied  particles,  and  therefore  the  factor  is  lower;  about 
1.35.  The  thickness  of  the  silver  planchet  used  for  backing  material  was  about  140 
mg/cm^ ,  which  is  well  beyond  the  saturation  thickness  of  ^  ^  P  and  very  close  to  the 
saturation  thickness  of  ^  ^  P. 

!•  Counting  efficiency 

The  ovCTall  counting  efficiencies  ^  ^  e  and  ^  ^  e  were  experimentally  obtained  by 
counting  sources  of  constant  activity  and  of  varying  thickness.  Previous  studies  by 
Nervik  and  Stevenson  (1952)  have  shown  that  the  efficiency  curves  are  not  simply 
described  by  equation  (1.3),  which  models  the  effect  of  self-absorption  only,  because  self¬ 
scattering  in  thin  sources  can  also  be  an  important  process.  Their  study  showed  that  the 
shape  of  the  efficiency  curve  is  a  function  of  the  composition  of  the  source  (nature  of  the 
precipitate)  and  the  eno'gy  of  the  beta  emitter.  For  the  weak  emitters  ^  ’  S  (E  =  0.167 
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MeV)  and  ‘  ^  Pm  (E  =  0.244  MeV)  the  efficiency  curves  exhibit  a  maximum  for  very 

mu 

thin  sources  (about  1-2  mg/cm^ )  followed  by  an  exponential  decrease  with  thickness. 
For  the  high  energy  emitters  ’  ^  P  and  ’  °  Y  (E  =  2.28  MeV),  the  curves  increase  slightly 

mu 

initially  and  reach  a  plateau  (Nervik  and  Stevenson,  1952).  The  transition  in  the  curve  for 
both  hard  and  soft  beta  emitters  is  due  to  the  relatively  higher  importance  of  self¬ 
scattering  versus  self-absorption  in  the  thin  sources.  The  transition  is  also  dependent  on 
the  (imposition  of  the  source.  The  exponential  decrease  is  due  to  self-absorption  only 
and  can  be  described  by  equation  (1.3). 
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a)- Counting  effldeiicy  of  P 

Pho^liorus-32  was  amuneicially  available  DuPont  in  a  solution  of 
KHa  PO4 .  The  ’  ^  P  purchased  was  not  a  standard  and  had  to  be  standardized. 
Verification  of  the  purity  of  the  tracer  was  obtained  by  determining  the  half-life  and  the 
half-thickness  of  absorption  on  the  low  level  beta  counting  system.  The  half-life  was, 
within  the  eircn*,  14.28  days  and  the  half-thickness  was  71.6  mg/cm^ .  The  solution  was 
then  standardized  by  counting  a  thin  source  of  ^  ^  P,  evaporated  on  a  stainless  steel  plan- 
chet  on  a  2n  gas-flow  proportional  counter  (Nuclear  Measurements  Corporation,  model 
PCC_11T).  The  counter  had  been  calibrated  by  the  manufacturer  with  sources  of  beta 
emitters  of  different  E  on  stainless  steel,  and  I  recalibrated  it  with  standards  of  ^  ^  ^  Th 
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and  ^  ^  Pm  (Amersham).  A  standard  of  ^  ^  ^  Hr  was  prepared  by  separation  from  a  high 
purity  U3  Os  standard  solution  (NBS  950A)  in  which  ^  ^  ^  Th  had  grown  into  radioactive 
equilibrium  (Fleer,  1991).  My  calibration  agreed  with  that  of  the  manufacturer  and  a  72 
%  counting  efficiency  was  determined  for  ^  ^  P  on  the  2n  gas-flow  proportional  counter. 
The  standardization  of  the  ^  ^  P  in  KHz  PO4  solution  from  DuPont  could  then  be  done. 

A  known  aiiwunt  of  activity  of  ^  *  P  was  precipitated  as  NH4  MgP04  .6H2  O  with 
varying  amounts  of  stable  PO4  .  The  sources  were  dried  until  the  weight  was  constant 
then  sealed  and  counted.  Phosphorus-32  is  a  hard  beta  emitter  and  the  efficiency  was 
almost  constant  in  the  range  of  thicknesses  encountered  in  all  the  samples.  The  efficiency 
for  an  infinitely  thin  source,  Co  >  was  calculated  by  evaporating  drops  of  a  standard  solu¬ 
tion  of  ^  ^  P  on  an  area  of  2.75  cm^  centered  on  the  Mylar  film.  This  was  to  ensure  that 
all  the  activity  was  distributed  on  an  area  equal  in  size  to  the  area  of  the  precipitate,  be¬ 
cause  the  counting  efficiency  is  not  the  same  in  all  the  area  of  the  source  (Lai  and  Schink, 
1960).  A  Millipore  filter  and  a  silver  planchet  were  placed  on  top  of  the  source.  The 
efficiency  fw  an  infinitely  thin  source  and  for  all  thicknesses  from  2  to  8  mg/cm^  was 
49%. 
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b)*  Counting  effidency  of  P 

A  solution  of  KH2  PO4  of  unknown  activity  of  ^  ^  P  was  obtained  from  DuPont. 
A  standard  of  '  ^  ^  Pm  was  available  from  Ameisham  and  was  used  as  an  analogue  of 
*  '  P.  ‘  ’  Pm  (t^^^  =  2.6234  y)  has  an  end  point  energy  of  0.224  MeV  which  is  extremely 
close  to  that  of  ^^P  (E  »  0.249  MeV).  It  is  therefore  expected  that  the  two 

aux 

radioisotopes  will  have  a  similar  half-thickness  of  absmption.  A  '  ^  ^  Pm  standard  was 
evi^rated  on  a  silver  jdanchet  and  its  absorption  curve  was  generated  (Figure  1.3).  The 
exponent  of  that  curve  is  0.156  cm^  /mg  and  the  half-thickness  of  absorption  is  4.4  ±0.1 
mg/cm^.  Within  experimental  errors  it  is  identical  to  4.3  ±0.1  mg/cm^ ,  the  half¬ 
thickness  of  *  ^  P.  Thus  ‘  ’  Pm  and  ^  ^  P  were  considered  to  be  equivalent  with  respect 
to  their  properties  of  interaction  with  matter.  A  known  activity  of  ’  ^  ^  Pm  was 
eviqrorated  on  a  Mylar  film.  A  filter  and  a  silver  planchet  were  put  on  top  so  that  this 
standard  was  as  similar  to  the  sample  as  possible.  The  ratio  of  the  net  count  rate  to  the 
activity  of  this  standard  was  the  counting  efficiency  of  ’  ^  P,  ^  ^  Co ,  for  an  infinitely  thin 
source.  To  calculate  ^  ^  e  fw  all  thicknesses,  a  carrier  free  unknown  activity  of  ^  ^  PO4 
was  precipitated  as  MgNH4  PO4  .6H2  O  with  increasing  amounts  of  stable  PO4  .  A  curve 
was  generated  o'  ^iie  net  count  rate  as  a  function  of  thickness.  The  expected  decreasing 
curve  was  regressed  widi  a  least  square  method  and  the  intercept  at  0  mg/cm^  was  deter¬ 
mined.  That  intercept  was  assumed  to  be  the  net  count  rate  for  an  infinitely  thin  source, 
^  ^  Co .  Since  from  our '  *  ’’  Pm  standard  ’  ’  eo  was  already  known  (i.e.  26%)  the  rest  of 
the  curve  could  then  be  calibrated.  The  curve  did  not  exhibit  any  maximum,  possibly 
because  (1)  I  could  not  make  sources  thinner  than  1  mg/cm^  or  (2)  self-scattering  is  small 
in  MgNH4  PO4  .6H2  O.  I  concluded  that  the  efficiencies  might  be  underestimated  by  up 
to  3  %  when  using  equation  (1.3)  to  calculate  the  efficiency  of  ^  ’  P.  Figure  1.4  shows  the 
counting  efficiency,^  ^  e,  as  a  function  of  sample  thickness.  The  counting  efficiencies  for 
an  infinitely  thin  source,  ’  ^  eo  ,  were  the  same  within  1  %  for  all  4  counters.  The  relative 
efficiency  curve  was  fitted  by  the  function  given  in  equation  (1.3). 
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Figure  U.  Net  count  rate  of  Pm  as  a  function  of  external  absorber  diickness. 
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Figure  1.4.  Relative  counting  efficiency  of  P  as  a  function  of  source  thickness. 
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The  exponent  p  tlM  curve  was  determined  for  each  detecuxr.  It  varied  fn»n 
0.145  to  0.160  and  averaged  O.ISO  +  0.01  cm^  Ang.  These  values  all  agreed  within  the 
error  of  the  regression.  Figure  1.5  shows  the  effect  of  a  change  in  p  on  the  relative  ef¬ 
ficiency.  The  error  on  the  absolute  efficiency  is  estimated  to  be  less  than  2  %  for  an  error 
of  0.01  cm^  Ang  on  p.  The  counter  efficiency  e  (expressed  in  %)  for  a  sample  of  thick¬ 
ness  X  (expressed  in  mg/cm^ )  is  calculated  from  die  experimental  relationship: 

(1.4)  €=26  (l-exp(-0.150  x))/ (0.150  x) 

2*  Yield  of  procedure 

The  chemical  yield  was  calculated  as  the  rado  of  the  amount  of  stable  P  in  the 
source  relative  to  the  amount  of  stable  P  initially  present.  The  determination  of  the  yield 
was  different  according  to  the  sample  type.  For  rainwater  samples  100  pmoles  KHa  PO4 
were  added  to  the  rainwater.  The  amount  of  PO4  present  in  the  final  NH4  MgP04 .6H2  O 
precipitate  was  determined  colorimetricaly  with  the  classic  molybdenum  blue  method 
(Murphy  and  Riley,  1962;  Koroleff,  1983).  This  procedure  was  used  because  the  amount 
of  phosphate  in  the  source  could  not  be  precisely  determined  from  the  weight  of  the 
precipitate.  This  was  because  the  precipitate  is  hygroscopic  and  the  stoichiometry  is  not 
always  6H2  O.  For  plankton  and  suspended  particulate  samples,  the  initial  amount  of 
PO4  was  measured  colorimetrically  after  digestion  with  nitric  acid.  In  the  case  where  the 
amount  of  phosphate  exceeded  100  pmoles,  the  sample  was  split  so  that  thin  sources 
could  be  obtained  to  minimize  self-absorption  of  ^  ^  P  in  the  source.  Experience  showed 
that  50-100  pmoles  of  stable  P  in  the  source  optimized  the  total  net  count  rate  and  mini¬ 
mized  self-absorption  of  ^  ^  P  for  samples  collected  in  the  Sargasso  Sea.  Usually,  for 
plankton  samples  obtained  from  15-20-minute  tows,  the  initial  amount  of  stable  P 
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Figure  1.5.  Relative  counting  efficiency  as  a  function  of  source  thickness  calculated 
for  different  values  of  the  absorption  coefficient  ti. 
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exceeded  100  (imoles  and  the  sample  was  split  For  suspended  particles,  the  sanq>les  had 
typically  30-60  pmoles  of  PO4  initially  and  were  dius  never  split 

The  data  for  all  the  rain  samples  and  plankton  samples  are  listed  in  Appendix  A. 
The  data  include  the  equations  of  the  regressitms  of  the  curves  Nt  and  Na  in  tables  A.1 
(rain  sanq>les)  and  in  table  A.2  (plankton  san^les).  The  thickness,  d,  the  yield,  y,  the 
counting  efficiencies,  ^  ^  e  and  ^  ^  e,  and  the  counting  rates  at  the  time  of  sampling  ^  ^  N 
and  ^  ^  N  are  given,  for  all  the  samples,  in  appendix  A  in  tables  A.3  for  the  rain  samples, 
in  tables  A.4  to  A.9  for  the  plankton  saiiq)les  and  in  table  A.IO  for  the  salt  samples  (see 
chapter  3). 

Vn  COUNTING  STATISTICS  AND  ERROR  ANALYSIS 

The  errors  on  the  absolute  initial  activities  ^  ^  A**  and  ^  ^  A°  were  computed  fcv 

the  time  the  sample  was  put  in  the  counter.  These  activities  are  derived  from  the  two 

experimental  curves  G  (t)  and  G  (t)  of  the  gross  count  rate  of  the  source  as  a  function  of 

t  ft 

time  expressed  in  equations  (1.1)  and  (1.2). 

!•  Error  on  the  net  count  rates 

The  curves  N^(t)  and  (t)  were  fitted  with  a  least  square  method  and  the  2  inter¬ 
cepts  N  ®  and  N  °  were  determined.  The  errors  were  given  by  the  non-weighted  least 
fl  t 

squares  frt  (samples  are  counted  in  time  typically  until  each  count  in  time  is  300  or  500 
for  low  count  rate  samples  and  1000  or  more  for  rain  samples).  ^  ^  N®  is  equal  to  (G®  - 
B)  /  ’  *  P  or  N^®  /  ’  ^  p.  The  error  on  ’  ^  N® ,  03  2  was  deduced  from  the  propagation 
theory  (Friedlander  era/.,  1981).  ’^N®  is  a  function  01  ^  ^  P,  Ga  and  B.  Thus  the  stan¬ 
dard  error  is  given  by: 

(1-5)  o,  ,  ^  -  a\  ,  p*  O’  ’  N«  O”  p)’  +  ‘(a’  ’  N’  )* 


37 


+  o^g*0’^N®/aB)^ 

Both  the  background  and  ^  ^  p  are  well  known,  and  the  errors  associated  with  these  quan¬ 
tities  are  small.  It  follows  that: 

(1.6)  03  2  jjo  "  =  <y  qo  *0’ '  N"  /dG°  *0’ '  /aN® 

a  a 

or 

(1-7)  =  1/”P‘0(^ 

a 

SinceN®  +”N®,thenN®  =N®  /  ^  ^  p +  ”  .  Therefore: 

t  t  a  ^ 

(1.8)  “N®  =N®  /^^p. 

t  a 

Applying  the  theory  of  propagation  of  error  to  ^  ^  N®  we  obtain: 

(1.9)  03  3  *  (3'  ’  N®  /  aN^  /  +  *  (a”  /  aN" 

a  t 

Thus: 

(1.10)  aj,(jo"=  (!/’ ’ P)’ •  +0"^ 

a  t 

The  standard  errors  ®3  2  and  03  3  were  therefore  calculated  from  equations  (1.7)  and 

( 1.10)  and  depend  only  on  the  errors  andc^^  which  were  determined  by  the  fit  of  the 

a  t 

time-dependent  curves  and  N^. 

2-  Error  on  the  absolute  activities 
”A°=”N^/y*”e 

R  =  ”A°  /^^A®  =(”N®  /^^N®)*(’2£/^^  e) 

RxMn  the  propagation  theory  it  is  straightforward  to  derive  the  errors  on  ^  ^  A®  ,  ^  ’  A° 
andR.  The  errors  are  given  by: 

(1.11)  03  2  ^0  "  =  d/y*’ '  e)"*03  2  jsjO  '  +  (’ '  N®  /y'*' '  e)"*Oy'  + 

(3  2N«/y*3  2eY*a3  2^^ 
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(1.12)  o>  j  ^0  ’  -  d/y*’  ’  e)’*Oi  j  ^  +  (”  N”  /y’*”  £)'*o^’  + 
(”N»/y”eY*o,,^' 

(U3)  o^’-  (”e/’'N"*”e)^*o,,|^'+  (”)(•/”  e*”N®)**ojj 

(”  N«  ♦”  e/>  >  N»  •”  cV*03 ,  N*  •”  e/”  £•”  N*  V*o, ,  ' 

In  practice,  the  errors  were  given  by  equations  (1.11),  (1.12)  and  (1.13).  For  the  rain 
samples  the  errors  associated  with  the  counting  statistic  were  small  (usually  <  2  %),  and 
thus  the  errors  on  the  initial  activities,  ^  ^  A°  and  ^  ^  A° ,  and  the  ratio,  R,  were  domina^d 
by  the  errors  on  the  counting  efficiencies.  For  plankton  samples  the  errors  were  often 
equally  dominated  by  both  counting  statistics  and  counting  efficiencies. 

vra  CONCLUSION 

The  method  described  in  this  chapter  allows  the  measurement  of  cosmogenic  ^  ^  P 
and  ^  ^  P  in  rainwater  (10  to  60  i ),  suspended  particulate  matter  filtered  from  2,000  to 
5,000  i  of  seawater,  and  15  to  20  minute  plankton  tows  in  an  oligotrophic  ocean.  The 
activities  of  ’  ^  P  and  ^  ’  P  in  rainwater  can  be  determined  within  10%-20%.  The  specific 
activities  in  both  suspended  particles  and  plankton  tows  can  be  determined  within  10%- 
30%.  The  activity  ratio  ’  ^  P/*  ^  P  was  measured  within  10%  in  rainwater  and  20-30%  in 
particles.  The  successful  measurement  of  the  weak  emitter,  ^  ^  P,  in  particulate  phases 
was  due  to  the  low  amount  of  stable  P  in  the  samples  collected  off  Bermuda.  Low 
amounts  of  stable  phosphate  are  essential  to  minimize  self-absorption,  in  the  source,  of 
the  low  energy  P  particles  emitted  by  ^  ^  P.  The  accuracy  of  the  measurements  is  limited 
by  the  accuracy  of  the  determination  of  the  counting  efficiencies.  The  precision  of  the 
measurements  is  limited  by  the  precision  of  the  counting  efficiencies,  the  chemical  yield 
and  the  counting  statistics  for  the  low  activity  samples. 
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Chq>ter  Two 


DETERMINATION  OF  THE  FALLOUT  RATES 
OF  AND  AT  BERMUDA 


I  INTRODUCTION 

In  ofder  to  understand  the  pathways,  concentrations  and  transformations  of  ^  ^  P 
and  *  *  P  in  the  upper  ocean  it  is  necessary  to  know  the  supply  to  the  ocean.  Following 
their  production  by  q)allation  in  Ar,  ^  ^  P  and  ^  ^  P  are  believed  to  be  rapidly  oxidized  to 
friio^hate  and  scavenged  by  sub-micrometer  aerosols.  The  aerosols  are  ultimately 
removed  frcmi  the  troposphere  by  precipitation.  Wet  deposition  is  believed  to  be  the 
primary  removal  mechanism  of  ^  *  P  and  ’  ’  P  (Lai  et  al.,  1957;  Goel  et  al.,  1959;  Lai  et 
al.,  1960).  Therefore,  in  mder  to  quantify  the  deposition  rates  of  ^  ^  P  and  ^  ^  P  at  a  par¬ 
ticular  site,  it  is  necessary  to  study  the  activities  of  ^  ^  P  and  ’  ^  P  in  rainwater  at  that  site. 

A  year-long  record  of  the  activities  of  ’  ^  P  and  ^  ^  P  and  of  the  activity  ratio 
3  3  py3  2  p  in  individual  precipitation  events  was  acquired  at  Bermuda  (32°  30'  N;  64° 
40*  W).  In  ccmjunction  widi  those  data  a  few  individual  rain  samples  were  collect  at 
Woods  Hole  (41°  31’  N;  70°  40'  W).  The  activities  of  ’  *  P  and  ’  ^  P  measured  in  rain- 
wata-  and  the  measured  rainfall  rates  allowed  computation  of  the  wet  deposition  rates  of 
’  *  P  and  ’  ^  P  from'March  1991  to  March  1992  at  Bermuda. 


n BACKGROUND 


Cosmogenic  ^  ^  P  and  ^  ^  P  are  ccmtinuously  produced  by  cosmic-ray  q)allation 
of  aigon  within  the  earth's  atmosphere.  The  production  rates  vary  notably  with  altitude 
and  latitude  (Lai  et  al.,  1958;  Lai  and  Peters.  1967).  The  production  rates  bodi  ^  ^  P 
and  ^  ^  P  increase  by  two  orders  of  magnitude  from  ground  level  to  the  top  of  the  tropo- 
q>here  and  are  much  higher  in  the  stratosphere  than  in  the  troposphere  (Lai  and  Peters, 
1967).  The  global  tropospheric  production  rates  of  ^  ^  P  and  ^  ’  P  are  almost  independent 
of  latitude.  This  is  in  ctmtrast  with  the  global  stratospheric  production  rates,  which  are  an 
order  of  magnitude  higher  at  a  latitude  of  90°  than  at  0°  (Lai  and  Peters,  1967).  The 
ratio  of  the  production  rates  of  the  two  isotopes  is  approximately  constant  with  time, 
latitude  and  altitude  (Lai  and  Peters,  1967).  The  atom  ratio  ^  ^  P/^  ^  P  of  the  production 
rates  has  been  determined  experimentally  in  argon  to  be  0.8  (Lai  et  al.,  1960a). 

Following  their  production  ^  ^  P  and  ’  *  P  are  believed  to  be  rapidly  oxidized  to 

phosphate.  Cbsroogenic  ’  *  P,  ^  ^  P  as  well  as  ’  Be  (t  =  53.3  d),  another  cosmogenic 

1/2 

nuclide  with  similar  source  and  pattern  of  production  as  ^  ^  P  and  ’  ^  P,  are  believed  to  be 
rapidly  scavenged  by  sub-micrometer  particles  in  the  atmosphere.  Tropospheric  aerosols 
become  condensation  nuclei  and  the  radionuclides  are  ultimately  deposited  by  precipita¬ 
tion.  Dry  deposition  of  both  short-lived  and  particle-reactive  radionuclides,  i.c.  ^  ^  P,  ^  ^  P 
and  ^  Be,  has  generally  been  assumed  to  be  small  relative  to  wet  deposition  (Lai  and 
Peters,  1967;  Young  and  Silker,  1974;  Turekian  et  al.,  1983).  Recent  studies  have  shown 
that  dry  deposition  of  ^  Be  is  only  10%  of  the  total  deposition  rate  (Brown  et  al.,  1989). 
Wet  deposition  rates  and  total  deposition  rates  of  ^  Be  have  been  measured  extensively 
(Rama  Thor  and  Zutshi,  1958;  Goel  et  al.,  1959;  Lai  et  al,  1960b;  Walton  and  Fried, 
1962;  Schumann  and  Stroeppler,  1963;  Young  and  Silker,  1974;  Lai  et  al,  1979;  Young 
and  Silker,  1980;  Kiishnaswami  et  al.,  1980;  Crecelius,  1981;  Turekian  et  al,  1983; 
Olsen  et  al.,  1985,  1986;  Dibb,  1989).  This  is  in  contrast  with  the  small  amount  of  data 
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on  the  ftUout  of  ’  ^  P  (Lai  et  cd.^  1957;  God  et  al.,  1959;  Lai  et  al.,  1960)  and  even  fewer 
data  for  ^  ’  P  fallout  rates  (Lai  et  al.,  1957;  Goel  et  al.,  1959).  The  lack  of  data  on  the 
deposition  rates  of  ^  ^  P  and  of  ^  ^  P  is  due  primarily  to  the  difficulty  of  the  measurements. 
It  results  in  a  poor  ladtudinal  coverage  of  the  wet  deposition  rates  of  ^  ^  P  and  ^  ^  P. 
There  have  been  mcxe  extensive  measurements  of  ’  ^  P  activities  in  surface  air  relative  to 
rainwater,  and  it  has  been  used  as  a  tracer  of  atmo^heric  circulation,  usually  in  adjunc¬ 
tion  with  ^  Be  and  sometimes  vdth  ^  ^  P  (Rama  and  Honda,  1961;  Aegerter  et  al.,  1966; 
Luyanas  et  al.,  1970;  Bhandari  et  al.,  1970;  Reiter  et  al.,  1971;  Rdter  et  al.,  1975;  Sanak 
etal.,  1985). 

Cosmogenic  ^  Be  has  been  widely  studied.  The  studies  of  its  distribution  in  air, 
scavenging  rates,  and  seasonal  and  latitudinal  variability  can  be  useful  for  the  study  of 
’  ^  P  and  ^  ^  P  and  those  studies  are  thus  reviewed.  The  strong  dependence  of  ’  Be  con¬ 
centrations  in  surface  air  on  latitude  and  season  is  now  well  established.  It  has  been 
attributed  to  latitudinal  and  se^.sonal  variations  in  the  exchange  between  the  troposphere 
and  the  stratosphere  (Parker,  1962;  Schumann  and  Stroeppler,  1963;  Peirson,  1%3; 
Rangarajan  and  Gopalakrishnan,  1970;  Reiter  et  al.,  1971;  Viezee  and  Singh,  1980; 
Dutkiewicz  and  Husain,  1985;  Dibb,  1989).  Maximum  ’’  Be  concentrations  in  surface  air 
and  maximum  deposition  rates  have  been  reported  at  mid-latitude  (40-50  °  N)  in  late 
spring/summer  (Parker,  1962;  Schumann  and  Stroeppler,  1963;  Viezee  and  Singh,  1980; 
Dibb  1989);  at  low  latitudes  (20-30  ^  N)  in  spring  (Rangarajan  and  Gopalakrishnan,  1970; 
Viezee  and  Singh,  1980);  and  at  polar  latitudes  during  austral  summer  (Feely  et  al.,  1977; 
Maenhaut  et  al.,  1979;  Sanak  et  al.,  1985).  Anomalously  high  deposition  rates  of  ^  Be  in 
the  intertropical  convergence  zone  (ITCZ)  have  been  reported  north  of  the  ITCZ  (10- 
20**  N)  by  Young  and  Silker  (1974).  They  explain  these  high  deposition  rates  not  by 
stratospheric  intrusion  per  se,  but  by  the  intrusion  of  tropospheric  cumulus  and 
cumulonimbus  clouds  into  the  stratosphere  and  scavenging  of  the  stratospheric  aerosols. 
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At  low  latitudes  (20-30^ )  and  mid-latitudes  (40-50^  )  the  exchanges  have  been  attributed 
to  low-pressure  troughs  associated  with  a  folding  of  the  tropopause  (Danielson  and 
Mohnen,  1977;  Viezee  and  Singh,  1980).  These  areas  of  exchange  are  found  in  the 
vicinity  of  the  jet  streams  which  are  typically  located  at  30°  and  60°  in  both  hemispheres 
(Hasse,  1983).  The  process  of  exchange  in  the  jet  stream  is  due  to  large-scale  eddy 
transport  (Reiter,  197S).  At  polar  latitudes  the  stratospheric  contribution  could  be  con¬ 
siderable  because  of  (1)  the  much  higher  production  rates  of  ^  Be  in  the  stratosphere  at 
high  latitudes  than  at  low  latitudes  (Lai  and  Peters,  1967)  (2)  substantial  and  continual 
fluctuations  of  the  height  of  the  tropopause  at  the  poles  (Sanak  et  al.,  1985)  and  (3) 
breakup  of  the  polar  vortex  in  the  winter  and  spring  (Martell  and  I>revinsky,  1960). 

During  the  time  of  atmospheric  testing  of  nuclear  weapons,  the  distributions  of 
e  bomb-produced  radioisotopes  ’  °  Sr,  *  ”  Cs,  ^  ^  *  Pu,  ^  ^  ®  Pu  and  ^  ^  °  Pu  were  similar 
to  the  distribution  of  ^  ^  P,  ^  ^  P  and  ’’  Be,  because  they  were  injected  primarily  into  the 
stratosphere  and  at  high  latitudes.  The  tropospheric  contribution  of  the  tests  was  small, 
because;  (1)  the  tropospheric  tests  were  low  yield  (Perkins  and  Thomas,  1980),  and  (2) 
the  tix^spheric  radioactivity  so  produced  was  removed  rapidly  after  the  tests.  The  bomb- 
produced  radioisotopes  were  thus  good  tracers  of  stratospheric  air.  The  concentrations  of 
the  bomb-produced  nuclides  in  surface  air  and  their  transport  to  the  troposphere  exhibited 
strong  latitudinal  and  seasonal  variations  (Rangarajan  and  Gopalakrishnan,  1970;  Silker, 
1972;  Perkins  and  Thomas,  1980;  Staley,  1982;  Dutkiewicz  and  Husain,  1985).  The 
concentrations  of  ’  ^  ^  in  surface  air  were  highest  in  the  spring  at  low  and  mid-latitudes 
(Rangarajan  and  Gopalakrishnan,  1970;  Silker,  1972).  The  fallout  rate  of  °  °  Sr  and  its 
concentration  in  surface  air  were  found  to  be  highest  (1)  in  late-winter  and  spring  at  low 
latitudes  (Machta,  1959;  Staley,  1982)  and  (2)  in  sununer  at  mid-latitudes  (Schumann  and 
Stoeppler,  1963;  Staley,  1982). 

Most  of  the  exchange  processes  between  the  stratosphere  and  the  troposphere 
evidenced  with  radionuclides  are  due  to  large-scale  eddy  transpon  across  the  tropopause 
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in  die  jet  stream  legum.  This  process  has  been  widely  reported  to  have  a  very  strong 
seasonality.  Nevertheless,  there  are  also  other  exchange  mechanisms,  in  particular  the 
adjustment  of  the  height  of  the  tropopause  (Reiter,  1975).  This  adjustment  is  continuous 
during  the  year  and  occurs  at  all  latitudes.  It  is  characterized  by  a  general  drc^  of  the 
trc^H^use  fnxn  spring  to  fall  and  a  general  rise  from  fall  to  spring. 

m  METHODS 

Rain  water  was  collected  from  the  roof  of  Hansen  Hall  at  the  Bermuda 
Biological  Station  for  Research.  The  roof  (area  of  about  30  m^ )  thus  collected  particles 
in  addition  to  rainwater.  Two  60-liter  barrels  were  connected  to  the  gutter  and  could  be 
alternately  filled.  Usually  only  major  rain  events  (rainfall  rate  >  1  cm/d)  were  sampled 
because  of  the  large  minimum  volume  required  (10-20 1 ).  The  rainwater  was  spiked  with 
stable  P  (KHi  PO4  )  as  carrier  phase  and  yield  monitor.  Phosphate  was  efficiently  ex¬ 
tracted  on  alumina,  and  the  sample  was  analysed  as  described  in  (Chapter  One.  In 
addition,  rainwater  samples  were  collected  from  a  4  m^  rain  collector  placed  on  the  roof 
of  the  Clark  Building  at  Woods  Hole  in  the  winter/spring  of  1990  and  1991.  The  collec¬ 
tion  procedure  used  was  identical  to  the  one  described  above. 

IV  RESULTS 

A  year-long  record  of  the  concentrations  of  P  and  P  and  of  the  activity 
ratio  ’  ’  P^  *  P  in  single  precipitation  events  was  acquired  at  Bermuda  (32°  30’N;  64° 
4(y  W).  Figtire  2.1  shows  the  record  of  the  ratio  fi’om  March  1991  to  March  1992. 
Figure  2.2  shows  the  daily  rainfall  rates  measured  from  January  1991  to  March  1992  at 
Saint  Davids,  situated  at  the  eastern  tip  of  the  island.  The  distribution  of  the  monthly 
averaged  rainfall  rate  during  that  period  shows  high  rates  in  February,  March,  September 
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of  1991  and  Januaiy,  February  of  1992  (Figure  2.3).  The  activities  of  ^  ^  P  and  ^  ^  P  and 
the  activity  ratios  ^  ^  P/^  ^  P  of  all  the  individual  rain  samples  collected  in  1990-92  at 
Woods  Hole  and  Bermuda  are  presented  in  tables  2.1  and  2.2.  The  activities  of  ^  ^  P  and 
’  ’  P  vary  from  0.15  to  3.5  dpm/1  for  ^  ^  P  and  from  0. 14  to  3.0  dpm/1  for  ’  ’  P  at  Bermuda 
and  from  0.4  to  0.9  dpm/l  for  ^  ^  P  and  from  0.5  to  0.98  dpm/1  for  ’  ’  P  at  Woods  Hole.  In 
contrast  to  the  wide  range  in  activities,  the  range  in  the  ratio  is  smaller  0.66  to  1.20  in  the 
Bermuda  samples  and  clustering  around  1.10  at  Woods  Hole.  The  activities  of  ^  ^  P  and 
^  ^  P  in  rainwater  are  extremely  variable,  because  they  depend  on  the  residence  time  of 
aerosols,  the  activities  of  ^  ^  P  and  *  ^  P  in  air  and  the  rate  of  scavenging  of  ’  ^  P  and  ^  ’  P 
by  aerosols.  In  contrast,  the  activity  ratio  ’  ’  P/*  ^  P  in  rain  depends  only  on  the  history 
and  pathways  of  the  air  mass.  The  highest  activities  of  ^  ^  P  and  ^  ^  P  in  Bermuda  rain  are 
found  in  March  1991,  and  they  are  associated  with  high  ratios  averaging  1.0  in  March 
1991.  The  June-November  1991  period  is  characterized  by  a  higher  frequency  of  low 
activity  ratios  than  the  March-April  1991  and  January-March  1992  periods.  At 

the  same  time,  ratios  as  high  as  1.0- 1.2  are  observed  during  the  period. 

An  important  feature  of  the  Bermuda  data  set  is  the  systematic  variation  in  the 
activity  ratio  ^  ^  p/^  *  p  about  the  annual  mean  of  0.96  (Figiue  2.4).  Higher  ratios  were 
observed  in  April  1991  and  January-February  1992,  and  lower  ratios  were  observed  in 
September-October  1991.  The  activity  ratio  was  on  average  1.01+0.04  in  the  late 
winter/spring  period  (i.  e.,  March  to  June  1991),  0.88+0.03  and  0.90+0.06  in  summer  and 
fall  1991  and  1.03+0.05  in  the  winter  of  1992  (Table  2.3).  The  variations  might  be  due  to 
stratospheric  input  in  spring  and  winter  at  Bermuda.  It  might  also  be  due  to  increased 
mixing  between  the  upper  and  lower  troposphere  in  sununer  and  fall  relative  to  spring 
and  winter. 
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Figure  2.1.  Activity  ratio  ”p/  P  in  major  precipitation  events  collected  at  Bermuda  Biological 
Station  f<w  Research  (BBSR)  from  March  1991  to  March  1992. 
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Table  2.1.  Activities  of  ^  and  and  A)  and  activity  ratio  (R) 

in  rain  samples  collected  during  1991  and  1992  at  the  Bennuda  biological  Station. 


dale 

t^nn/l 

”a 

t^nn/l 

R 

dpm/(4im 

0S>Mar-91 

0.41±0.07 

0.37±0.06 

0.91  ±0.14 

lO-Mar-91 

1.92±p.30 

2.12±0.29 

1.10±0.12 

lO-Mar-91 

3.46±Q.16 

3.02±0.46 

0.87  ±0.10 

20-Mar-91 

1.05±0.19 

1.13±0.14 

1.08  ±0.12 

20-MV-91 

1.55±0.3 

1.69±0.20 

1.09  ±0.11 

22-Apr-91 

1.02±0.11 

1.06±0.14 

1.04  ±0.12 

22-Apr-91 

0.55±0.07 

0.67+0.11 

1.22  ±0.14 

Ol-May-91 

0.75±q.l0 

0.73±0.10 

0.97  ±0.12 

ll-Jun-91 

0.60±0.09 

0.58±0.10 

0.96  ±0.13 

ll-Jun-91 

038±0.06 

039±0.06 

0.78  ±0.14 

20-Jun-91 

038±0.06 

0.64±0.08 

1.10  ±0.12 

27.Jun-91 

0.49±0.06 

0.50±0.07 

1.02  ±0.12 

07-Jul-91 

U9+P.17 

0.98±0.14 

0.71  ±0.11 

08-Jul-91 

1.31±0.13 

0.86+0.11 

0.66  ±0.10 

18-Jul-91 

0.72±0.07 

0.73±0.10 

1.02  ±0.14 

19-Jiil-91 

0.60±0.07 

0.70±0.10 

1.16±0.16 

07-Aug-91 

1.05±0.09 

0.86±0.14 

0.81  ±0.10 

26- Aug-91 

0.77±0.09 

0.75±0.10 

0.97  ±0.11 

OS-Sep-91 

1.04+0.13 

0.78+0.14 

0.70  ±0.09 

ll-Sep-91 

0.67±0.09 

0.59+0.09 

0.87  ±0.10 

14-Sep-91 

026±0.03 

026±0.04 

1.00  ±0.15 

16-Sep-91 

0.90±0.09 

0.72±0.09 

0.80  ±0.09 

29-Sep-91 

0.47±0.05 

0.48±0.04 

1.03  ±0.13 

28-Oct-91 

0.79±0.09 

0.62±0.08 

0.78  ±0.10 

28-Oct-91 

0.26±0.04 

0.22±0.03 

0.84  ±0.15 

07-NOV-91 

0.79±0.09 

0.77±0.09 

0.98  ±0.13 

17-NOV-91 

0.88±0.10 

0.81±0.09 

0.92  ±0.13 

16-Jan-92 

0.31±0.03 

0.31±0.05 

1.00  ±0.14 

21-Jaii-92 

1J6±0.15 

1.42±0.17 

1.13  ±0.12 

2S-Jan-92 

0.36±0.04 

0.39±0.06 

1.08  ±0.15 

31-Jan-92 

0.67±0.09 

0.74+0.09 

1.09  ±0.13 

02-Rb-92 

0J1±0.03 

0.38±0.04 

130  ±0.14 

14-Feb-92 

0.69±009 

0.63±0.1 

0.91  ±0.12 

17-Mar-92 

0.15±0.02 

0.14±0.03 

0.94  ±0.18 

26-Mar-92 

0.54±0.08 

0.56±0.09 

1.02  ±0.19 

26-Mar-92 

0.65±0.12 

0.6Q±0.12 

0.93  ±0.15 
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Table  U.  Aoiviiies  of  ^  and  and  and  activity  ratio  ^^P/*^(R) 

in  rain  samiries  ooUected  at  Woods  Hole. 


date 

”a 

dpm/l 

”a 

dpm/l 

R 

dpm/(4)m 

27-Dec-90 

0.82±0.09 

0.98±0.10 

\2l  ±0.09 

16-JM-91 

0.41±0.04 

O.Slip.06 

1.24  ±0.11 

21-Apr-91 

O.9Qt0.1O 

0.74i0.10 

0.82  ±0.13 

Table  13.  Seasonal  averages  of  the  activity  ratio  ^^P/^P  (R)  in  rain  and  associated  residence  times  of 
tropospheric  aerosols  at  Bermuda  for  two  different  values  of  Ro  . 


Season 

R 

dpm/dpm 

T 

days 

Ro«0.7 

T 

days 

Ro  bO.6 

Spring  and  March  1991 

1.01±0.04 

5Q±10 

11Q±2S 

Spring  1991 

1.02±O.OS 

S2±10 

115±35 

Summer  1991 

0.88i0.03 

25±5 

52±10 

Fhlll991 

0.90±0.06 

30±7 

57±10 

Winter  1992 

1.03±0.0S 

55±12 

95-290 

4Qt7 


Annual  average 


0.96±0.02 
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The  average  wet  deposition  rates  were  determined  by  multiplying  the  monthly 
average  activities  of  ^  ^  P  and  ’  ^  P  (Figure  2.S)  by  the  monthly  average  rainfall  rates. 
Monthly,  seasonal  and  yearly  averages  of  the  wet  deposition  rates  are  shown  in  tables  2.4 
and  2.5.  The  monthly  mean  wet  deposition  rates  of  ^  ^  P  and  ^  ^  P  show  significant  dif- 
feieiKes  frcmi  one  month  to  another,  mostly  due  to  the  difference  in  rainfall  rates.  The 
nxmthly  averages  vary  from  0.059  to  0.266  dpm/cm^  ^  for  ^  ^  P  and  from  0.058  to  0.264 
dpm/cm^  /yr  for  ^  ^  P.  The  extreme  upper  and  lower  values  occurred  in  March  aiul  May 
1991  (Table  2.4).  Fn*  the  9  other  months,  the  average  wet  deposition  rates  woe  much 
more  constant  and  averaged  0.074+0.04  dpm/cm^  /yr  for  ’  ^  P  and  0.073  ±  0.05 
dpm/cm*^  for  ^’P.  The  extremes  in  March  and  May  were  caused  by  (1)  a  low 
precipitation  rate  in  May  1991,  resulting  in  a  low  fallout  rate  for  that  month,  and  (2)  very 
high  concentrations  of  ^  ^  P  and  ^  ^  P  in  rain  and  high  rainfall  rates  in  March  1991.  The 
seasonal  averaged  wet  deposition  rates  range  from  3.6  10  to  10.2  10  dpm/cm  /yr  for 
’  ^  P  and  3.7  10  ^  to  1 1.3  10  ^  dpm/cm^  /yr  for  ’  ’  P  (Table  2.5).  The  annual  averages  are 
8.6  10  and  8.2  10  dpm/cm^  /yr  or  4.9  10  and  8.2  10  atom/cm^  /min.  Assuming  a  40 
day  residence  time  of  tropospheric  aerosols,  Lai  and  Peters  (1967)  estimated  global  mean 

*3  *3 

tropospheric  fallout  rates  of  5.8  10  atom/cm^  /min  for  ^  ^  P  and  6.9  10  atom/cm^  /min 
for  ^  ^  P  with  uncertainties  of  20  to  30%.  These  values  agree  within  the  uncertainties  of 
each  estimation. 

V  MODEL  OF  AEROSOL  RESIDENCE  TIMES 

In  this  section,  I  explain  how  the  residence  time  of  sub-micrometer  aerosols  in 
the  troposphere  can  be  derived  from  the  ratio  ^  ^  P/’  ^  P  in  individual  rain  events.  The 
ratio  ^  ^  P/’  ^  P  in  air  increases  between  the  time  the  nuclides  are  produced  and  the  time 
the  aerosols  are  scavenged  by  precipitation  because  of  the  differential  decay  rates  of  ^  ^  P 

4  A 

and  P.  If  one  assumes  that  there  is  no  fractionation  during  the  condensation  process. 
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thra  the  increase  in  the  ratio  in  the  air  will  be  exactly  reflected  in  the  rainwater.  The 
determination  of  the  residence  time  is  nnodei-dependent.  Both  a  non-steady-state  and  a 
steady-state  model  were  investigated.  A  non-steady-state  model  has  been  previously 
applied  to  cosmogenic  ^  Be  (Lai  and  Peters,  1967;  Shapiro  and  Forbes-Resha,  1976; 
Bleichrodt,  1978)  and  to  the  ratio  of  two  cosmogenic  radioisotopes:  (1)  ’  ^  P/  *  P  (Lai  et 
al.^  1957;  Goel  et  al.  1959;  Luyanas  et  al.,  1979;  Rama  and  Honda,  1961)  and  (2) 
’  Be^  ^  P  (Goel  et  al.,  1959;  Lai  et  al.,  1960;  Rama  and  Honda,  1961;  Walton  and  Ried, 
1962;  Bhandari  et  al.,  1970). 

a)  Non  steady-state  model 

The  mass-conservation  equation  for  a  radionuclide  in  a  parcel  of  air  is  the  fol¬ 
lowing: 

(2.1)  dOdt  =  aaat  +  u  aaax + v  aaay  +  w  acyaz  =  f  -  x  c 

Where  C  is  the  mean  concentration,  X  is  the  decay  constant,  u,  v  and  w  are  the  mean 
adveedve  velocities  in  the  x,y  and  z  dimensions,  and  F  is  the  mean  production  rate  in  the 
troposphere.  Equation  (2.1)  assumes  that  the  radionuclides  are  continually  produced  and 
lost  by  decay  in  the  air  mass.  No  removal  term  by  precipitation  is  included,  because 
precipitation  is  a  discontinuous  process  and  will  be  treated  as  such.  The  model  assumes 
no  mixing  of  the  air  mass  with  surrounding  air.  The  solution  of  (2.1)  is: 

(2.2)  C  =  F/X-(F/X-C°)  exp( -  X t) 

Where  C  =  C°  at  t  =  0.  Furthermore  it  is  assumed  that  at  t  =  0  both  radioisotopes  were 
efficiently  removed  by  washout.  Thus  =  0  and  equation  (2.2)  becomes: 

(2.3)  C  =  F/X  (1 -exp(-Xt)) 
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Apr.  Jim.  Aug.  Oct.  Dec.  Feb.  Apr. 


ratio  in  rain 
annual  average  0.96 


Figure  2.4.  Activity  ratio  *p/*P  in  rain  at  BBSR  averaged  monthly  from  March  1991 
to  March  1992. 
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Table  2^.  Mon^  ■venues  ^  the  a^viiies  of  and  (^A  and  ^A).  lainfall  laie  (r)  and  wet 
deposition  fates  of  Tand  fc  wand”  W)atBennuda  in  1991  and  1992. 


Month 

”a 

dpm/l 

”a 

dpm/i 

10  cm/min 

”w 

10**  dpta/aa^  fyr 

Mar-91 

1.46 

1.45 

3.47 

2.66 

2.64 

Apr-91<t> 

0.78 

0.86 

1.23 

030 

0.55 

May-91* 

0.75 

0.73 

0.15 

0X159 

0X158 

Jun-91 

0J3 

0.52 

1.97 

035 

034 

Jul-91 

0.99 

0.82 

1.85 

0.96 

020 

Aug-91 

0.91 

0.81 

1.42 

0.68 

0.75 

Sep-91 

0.67 

0.57 

339 

126 

1.07 

Oct-9 1+ 

0.52 

0.42 

2.46 

0.67 

0.54 

Nov-91 

0.83 

0.79 

0.94 

0.41 

0.39 

Jan-92 

0.70 

0.77 

3.14 

1.15 

127 

Feb-92 

0.50 

0.51 

3.49 

0.92 

0.93 

Mar-92 

0.74 

0.72 

1.48 

0.58 

0.56 

*  Avenges  based  on  one  data  poinL 
+  Avenges  based  on  duplicates  of  the  same  rain  event. 


Table2.5.  Seasonalwetdepositkxintesof^  Pand  Wand  W),  seasonal  activities^  A  and  A)and 
seasonal  ninlall  ntes  (r).  Spring  refers  to  April,  May  and  June;  summer  refers  to  July,  August  and 
September,  fall  refers  to  October,  November  and  December,  and  winter  refers  to  January  and  February. 


season 

“a 

I  ”w 

dpm/i 

dpm/l 

*2  2 

cm/mon  10  dpm/cm  fyr 

*3  2 

10  atom/cm  Anin 

Spring  and 
Mttch91 

0.94 

0.94 

738 

8.4±12 

8.4±2.2 

4.7±12 

8.4±22 

Spring  91 

032 

033 

4.80 

3.6£02 

3.7±1.2 

2.Qt03 

3.7±13 

Summer  91 

0.85 

0.72 

9.86 

102±2.2 

82±13 

S.8±12 

8.2±13 

FaU91 

0.73 

037 

5.77 

5.1±1.2 

4.7±1.3 

2.9±03 

4.7±1.4 

Winter  92 

0.58 

0.65 

14.33 

10.1±2.0 

11.3±23 

5.7±1.1 

11.3±22 

Annual 

avenges 

022 

0.77 

8.76 

8.6±13 

8.2t2.1 

4.9±1.4 

8.2±2.2 
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Equation  (2.3)  describes  the  evolution  of’^Por’^P  attached  to  aerosols  that  are  part  of 
an  air  parcel  and  thus  move  with  it  The  ratio  R  is  then  simply  given  by: 

(2.4)  R  *  C3  3  /C3  2  *  (F3  3  /Fa  a )  (^3  2  A3  3  )  (1  -  exp(-X3  3  t))/(l  -  exp(-X3  2  t)) 

R  is  equal  toRo  -  F33/F32  at  time  to  =  0  and  reaches  a  limiting  value  of  R  «  1.77 

c 

Fs  3  /ES  2  in  approximately  300  days.  In  the  troposphere,  since  the  residence  time  of  air 
masses  is  shorter  than  300  days  (Lai  and  Perters,  1967;  Poet  et  al.,  1972;  Moore  et  al., 
1973;  Martell  and  Moore,  1974;  Shapiro  and  Fcnbes-Resha,  1976;  Bleichrodt  1978; 
Holloway  and  Hayes,  1982)  the  ratio  ^  ^  P/^  ^  P  in  the  air  mass  will  increase  until  the 
aerosols  are  removed  by  condensation  and  ultimately  by  precipitation.  Therefore,  the 
ratio  ^  ^  P/^  ^  P  in  an  individual  rain  event  is  identical  to  the  ratio  in  the  tropospheric  air 
scavenged  by  the  rain,  and  equation  (2.4)  becomes  the  following: 

(2.5)  R^=(F33  /F32  )  1.77(1  .CXP(.X33  T))/(1-CXP(.X3  2  t)) 

Where,  R^  is  the  ratio  in  a  tain  event,  and  t  is  the  average  residence  time  of  tropospheric 
aerosols  in  an  air  mass  that  is  being  washed  out  by  a  particular  rain  event 

To  summarize,  the  assumptions  of  this  nnodel  are:  (1)  the  air  masses  are  con- 
rined  to  the  troposphere  only;  (2)  the  troposphere  is  well  mixed;  (3)  no  mixing  occurs 
during  the  movement  of  an  air  mass;  (4)  at  to ,  C=  0;  and  (5)  no  fractionation  of  ^  ^  P  and 
^  ^  P  occurs  during  condensation  on  the  aerosols.  The  residence  time  calculated  with 
equation  (2.5)  can  be  in  error  if  any  of  the  assumptions  is  violated.  This  will  be  the  case, 
for  example,  if  stratospheric  air  intrudes  and  mixes  with  the  tropospheric  air,  or  if  the 
lower  and  upper  troposphere  are  isolated  or  partly  isolated  from  one  another. 

The  most  in^Kutant  feature  of  the  Bermuda  data  set,  as  stated  previously,  is  the 
small  variability  of  the  activity  ratio  ^  ^  P/®  ^  P,  with  values  ranging  from  0.66±010  to 
1.22±0.14  (Table  2.1).  According  to  the  model,  the  ratio  in  rainwater  should  vary  from 
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an  initial  ratio  Ro  to  an  equilibrium  ratio  Re  equal  to  1.77  times  Ro  •  The  range  of  ratios 
observed  at  Bennuda  is  0.66*1.22.  The  highest  and  lowest  ratios  observed  are  in  a  ratio 
of  1.8  which  is,  within  the  error,  identical  to  1.77.  Thus  the  data  exhibit,  within  the  uncer¬ 
tainties  of  the  measurements,  the  entire  theoretical  range,  i.e.  Ro  to  R  .  In  view  of  the 

c 

fact  that  the  emnrs  on  the  measurements  are  around  0.1,  the  theoretical  range  of  the  ac¬ 
tivity  ratio  could  be  0.6  to  1.06  or  0.7  to  1.24  ch*  any  range  in  between  that  satisfies  the 

relationship R  /Ro  =1.77. 
c 

Equation  (2.5)  allows  the  calculation  of  the  residence  time  x  of  tropospheric 
aerosols  for  each  activi^  ratio  measured  in  individual  rain  events.  Figure  2.6  is 

a  plot  of  equation  (2.5)  for  different  values  of  Ro :  0.60,  0.65  and  0.70.  It  is  seen  in 
Figure  2.6  that  the  residence  time  of  tropospheric  aerosols  cannot  be  determined  precisely 
if  it  exceeds  80  days.  The  residence  time  calculated  for  the  two  extreme  values  of  Ro  is 
shown  in  Figure  2.7  and  Table  2.6.  It  is  high  on  average  in  March-April  1991  and 
January-February  1992  and  lower  on  average  in  August-October  1991.  The  activity 
ratios  are  also  averaged  over  a  season  and  the  corresponding  residence  time  calculated 
(table  2.3).  From  March  to  June  of  1991  and  January  to  February  1992,  the  residence 
time  averaged  52+10  days  while  from  July  to  November  1991  it  averaged  27±7  days, 
assuming  Ro  =  0.7.  If  one  assumes  that  Ro  =  0.6  then  the  residence  times  are  about  1 10 
days  and  55  days  for  the  two  periods.  A  value  of  Ro  =  0.7  gives  residence  times  that  are 
more  consistent  with  previous  estimate  of  30-40  days  (Beck  and  Kuroda,  1966;  Lai  and 
Peters,  1967;  Shapiro  and  Forbes-Resha,  1976;  Bleichrodt,  1978). 

by  Global  steady-state  model 

Another  method  of  estimating  the  residence  time  consists  of  establishing  a 
steady-state  mass  balance  for  each  isotope.  The  approach  assumes  that  the  mean  produc¬ 
tion  rate  of  an  isotope  in  tropospheric  air  masses  is  balanced  by  the  sum  of  its  decay  and 
its  renooval  by  precipitation: 
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(2.6)  F«W+XC 

(2.7)  x«(F-W)/(XW) 

Where 

W:  mean  wet  dqx>sitx>n  rate  of^^Por^^Pin  aumi  cm  min  ; 

*2  *1 

F:tr(^x>sphericpiodiictionrateof^^Por^^Pinatomcm  min  ; 

X:  decay  constant  of^^Por^^Pin  min 

C:  mean Par  ^’Pconcentration  in  tropospheric  air  in  atom/cm^ ; 
t:  mean  reridence  time  of  aerosols  in  the  troposphere,  t  ^C/W. 

The  monthly  mean  wet  deposition  rates  of  ’  *  P  and  *  ^  P  (i.e.,  W3  2  and  W3  3 ) 
are  given  in  table  2.4,  and  the  seasonal  means  are  given  in  table  2.S.  Equation  (2.7) 
allows  computation  of  x  fcM'  each  of  the  three  particular  periods  late  winter/spring  of  1991 
and  late  winter  1992  and  summer/fall  of  1991.  It  is  assumed  that,  on  these  time-scales, 
there  is  a  balance  between  the  production,  decay  and  removal  by  precipitation  of  ^  ^  P  and 
*  ’  P.  The  average  wet  deposition  rates  are  525  10  ^  and  9.65  10  ^  atom/cm^  /min  for  ’  ^  P 
and  *  ^  P,  reflectively,  in  late-winter/spring  of  1991  and  winter  of  1992.  For  summer  and 
fall  the  wet  deposition  rates  are  4.86  10  ^  and  7.3  10  ^  atom/cm^  /min  for  ^  ^  P  and  ’  ’  P, 
respectively.  The  global  tropospheric  production  rates  have  been  estimated  by  Lai  et  al. 
(1988)  to  be  1.6  10  ^  and  1.3  10  ^  atom/cm^  Anin  for  ^  ^  P  and  ^  ^  P,  respectively  cor- 
refxmding  to  the  years  1948-1949.  Since  the  years  1992  and  1948/49  correspond  to 
periods  of  low  production  rates  in  the  1 1-year  cycle,  the  estimates  of  Lai  et  al.,  (1988)  are 
used  unmodified.  The  residence  times  are  thus  calculated  for  each  isotope.  For  the 
period  March  to  June  1991  and  January  to  March  1992,  the  residence  time  calculated  for 
’  *  P  and  ^  *  P  are  42±12  days  and  13±5  days,  respectively.  For  the  period  July  to 
November  1991,  the  residence  times  are  47±15  days  and  28+10  days  for  ^  *  P  and  ’  ’  P, 
respectively.  The  residence  times  calculated  with  die  sready-state  model  are  similar  to  the 
residence  times  calculated  with  the  non  steady-state  model  applied  to  the  ratio  ^  ^  P/^  ^  P, 
suggesting  that  steady-state  is  a  good  approximation  on  a  seasonal  time-scale. 
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VI  DISCUSSION 


In  this  section  diree  aspects  of  the  data  will  be  discussed.  First,  the  observations 
made  both  at  Woods  H(de  and  at  Bermuda  will  be  compared  to  previous  data  on  ^  ^  P  and 
^  ^  P  in  air  and  rainwater  samples.  Second,  the  contrast  between  winter/spring  and 
summer^all  of  the  wet  deposition  rates  and  the  ratios  will  be  compared  to  the  general 
patterns  found  in  previous  studies  of  short>lived  cosmogenic  nuclides  and  artificial 
radionuclides.  Third,  the  nxne  detailed  structures  in  the  data  on  the  ratio  ^  ^  P/^  ^  P  in  rain 
will  be  discussed  as  well  as  the  validity  of  the  assumptions  made  in  the  non-steady  state 
model. 


A  range  of  0.66  to  1.24  was  observed  for  the  activity  ratio  ^  ^  P/^  ^  P  in  rainwater 
in  the  samples  collected  both  at  Bermuda  and  at  Woods  Hole.  The  observed  ratios  were 
somewhat  lower  than  the  ratios  of  0.8  to  1.6  measured  in  rainwater  by  Goel  et  al.  (1959) 
and  higher  than  the  ratios  of  0.42-0.56  determined  in  3  rainwater  samples  by  Lai  et  al. 
(1957).  Rama  and  Honda  (1961)  determined  the  ^  ’  P/^  ^  P  activity  ratio  in  air  at  different 
altitudes,  both  in  the  troposphere  and  stratosphere,  and  found  ratios  ranging  from  0.46  to 
1.1,  which  are,  within  the  uncertainties  of  the  measurements,  veiy  close  to  my  observed 
range  of  ratios.  Luyanas  et  al.  (1970)  reported  ratios  that  are  significantly  lower,  ranging 
from  0.18  to  1.01  in  surface  air  samples  collected  in  1967-1969.  They  argued  that  their 
lowest  values  0.18-0.26  were  probably  due  to  artificial  sources  of  ^  ^  P  from  weapon  tests 
and  that  they  should  be  neglected,  bringing  their  range  to  0.3-0.4  to  1.0.  Lai  and  Peters 
(1%2)  estimated  activity  ratios  in  tropospheric  air  ranging  from  0.46  to  0.81  using  a 
different  approach.  They  measured  a  ratio  of  the  production  rates  ^  ^  P/*  ^  P  in  Ar  of  0.82 
(atcnn  ratio)  or  0.46  (activity  ratio)  (Lai  et  al.,  1960a).  Their  results  suggest  that  the 
expected  range  of  ratios  is  lower  than  what  I  have  found  at  Woods  Hole  and  Bermuda. 
The  measurements  of  the  production  rates  of  ’  ’  P  and  ^  ^  P  made  by  Lai  et  al.  (1960a) 
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differ  from  die  calculation  done  earlier  by  Lai  (1958)  that  inedicted  a 
producdon-rate  ratio  of  1.1  (atom  ratio).  This  discrepancy  might  be  due  in  part  to  the 
difficulty  of  measuring  ^  ^  P,  since  very  few  direct  measurements  of  the  fHoduction  rates 
of’^Pand^’P  were  made.  It  appears  that  the  measured  ratios,  either  in  air  w  rain,  have 
a  wider  range  than  would  be  predicted  by  the  non-steady-state  model  (Ro  to  1.77  Ro ).  It 
is  probably  due  to  difficuldes  in  measuring  ^  ^  P  and  ^  ^  P.  It  might  possibly  be  because 
the  tnc^iospheie  is  not  always  well  mixed  but  divided  into  upper  and  lower  zones.  It  is 
conceivable  that  an  air  mass  could  spend  a  l<»ig  time  in  the  lower  troposphere  (which  has 
small  production  rates  of  ^  ^  P  and  ^  ^  P)  before  the  aerosols  are  being  scavenged  by 
precipitation.  Once  the  air  mass  is  confined  to  the  lower  troposphere  the  ^  ^  P/^  ^  P  ratio 
would  exponentially  increase  and  increase  much  faster  than  if  the  air  mass  was  irradiated 
in  the  whole  troposphere.  This  might  explain  why  high  ratios  (i.e.,  higher  than  the  ratios 
predicted  from  secular  equilibrium)  are  observed. 

A  contrast  that  is  possibly  seasonal  has  been  observed  in  the  activity  ratio 
3  3  py3  2  p  Bermuda.  High  ratios  were  observed  in  March- June  1991  and  January- 
March  1992  at  Bermuda.  Ratios  significantly  higher  than  the  annual  average  were 
observed  in  April  1991  and  January-February  1992.  Significantly  lower  ratios  were 
observed  in  September-October  1991.  Residence  times  based  on  the  ratio  were  high  in 
March-June  1991  and  Januaiy-March  1992  period  averaging  50!±10  days  while  they 
averaged  27±7  days  in  die  July-November  1991  period.  The  long  apparent  residence 
time  of  SO  days  was  higher  than  the  estimates  of  30-40  days  (Beck  and  Kuroda,  1966;  Lai 
and  Peters,  1967;  Shapiro  and  Forbes-Resha,  1976;  Bleichrodt,  1978).  The  observations 
might  be  explained  by  a  stratospheric  input  during  that  period,  "apparently"  increasing  the 
residence  time.  Strato^heric  air  is  charactoized  by  high  ’  ’  P/’  ^  P  ratios  because  both 
radioisotopes  are  in  equilibrium  with  their  production  rates  due  to  the  long  residence  time 
(about  14  months)  of  stratopheric  air.  The  data  suggest  that  the  activity  ratio  ^  ^  P/  ^  P  in 
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straiospheric  air  is  about  1. 1-1.2.  Intrusion  of  high  ^  ^  P/^  ^  P  air  would  explain  the  ap¬ 
parent  old  age  of  tn^XMpheric  air  masses  in  March-June  1991  and  January-Maich  1992. 
It  might  also  explain  the  high  ratios  observed  in  individual  rain  samples  in  summer  1991. 
The  observed  high  ratio  in  April  1991  is  consistent  with  evidence  of  a  seasonal  stratos¬ 
pheric  intrusion  of  ’  °  Sr  in  the  spring  at  Bermuda  (Health  and  Safety  Laboratcuy,  1977). 
Figure  2.8  shows  a  6-year  rectnd  of  ’  °  Sr  in  precipitation  with  a  maiked  peak  in  the 
firing  of  every  year.  The  high  ratios  and  deposition  rates  observed  in  the  winter  and 
spring  at  Bermuda  are  due  possibly  to:  (1)  stratospheric  intrusion  at  high  latitudes  (45-50 
**  N)  and  subsequent  advection  of  that  air  to  lower  latitudes  or  (2)  stratospheric  intrusion 
at  low  latitudes  (20-30  °  N).  The  southward  advection  of  high  latitude,  high  altitude  air 
masses  is  markedly  visible  in  the  distribution  of  ozone  in  the  atnK>spheie  of  the  northern 
hemisphere  (Johnson  and  Viezee,  1981).  It  is  furthermore  confirmed  by  trajectory 
analysis  which  shows  that  in  the  winter  most  of  the  air  masses  traced  10  days  before 
reaching  Bermuda  come  from  the  American  continent  at  both  high  latitudes  and  altitudes 
(J.  MerrUI;  AEROCE  meeting,  1991). 

A  change  in  the  convective  mixing  in  the  troposphere  can  explain  some  of  the 
contrast  observed  between  the  two  periods.  In  the  sununer,  convection  is  trxrre  efficient 
than  in  the  winter,  leading  to  more  extensive  mixing  between  the  lower  and  upper  tropo¬ 
sphere.  It  is  evidenced  by  ^  ‘  °  Pb  and  ^  ^  ^  Rn  determined  on  aerosols  at  the  AERCKE 
site  at  Bermuda  (western  tip  of  the  island).  The  ratio  ^  ®  Pb/*  ^  ^  Rn  is  higher  in  the 
summer  than  in  the  winter  because  of  increased  convection  in  the  summer  delivoing  high 
activities  of  ^  °  Pb  to  the  lower  troposphere  (W.  Graustein,  personal  comnununication  at 
AEROCE  meeting  1992).  The  magnitude  of  vertical  mixing  in  the  troposphere  affects 
the  activity  ratio  ^  ^  P^  ^  P.  Increased  mixing  between  the  lower  and  upper  troposphere, 
following  a  period  of  stable  stratification,  would  decrease  the  residence  time  determined 
from  the  activity  ratio  ^  ^  P/^  ^  P  in  rainwater. 
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The  lesidence  tune  27±7  days  estimated  in  the  summa-  and  fall  1991  at 
Bennuda  is  ccmsistent  with  previous  studies  of  shcfft-lived  radionuclides  (Bleichrodt, 
1978;  ShapiiD  and  Forbes-Resha.  1976;  Lai  and  Peters,  1967)  and  artificial  nuclides 
(Beck  and  Kuroda,  1966).  It  is  in  disagreement  with  most  of  the  studies  the  radon 
daughter  nuclides  which  tend  to  show  readence  times  on  the  cuder  of  a  week  (Poet  et  al., 
1972;  Mo(»e  et  al.^  1973;  Maitell  and  Moore,  1974;  W.  Graustein,  personal  oxnmunica- 
don  at  AEROCE  meeting,  1992).  The  disagreement  might  be  in  part  due  to;  (1) 
stratospheric  input  that  tends  to  apparendy  increase  the  residence  time  as  determined  by 
the  cosnwgenic  nuclides  and/on  (2)  poorly  mixed  troposphere  with  increasing  residence 
times  with  altitude.  The  determination  of  t,  as  explained  previously,  is  very  sensitive  to 
the  choice  of  the  initial  ratio  ^  ^  P/^  ^  P.  The  estimate  of  27  days  assumes  that  Ro  is  0.7. 
That  estimate  was  also  calculated  from  the  average  ratio  ^  ^  P/^  ^  P  in  the  rain  in  surruner 
and  fall  1991.  It  is  observed  in  figure  3.1  that  the  ratio  exhibits  sonte  high  values  in  the 
summer  and  fall  and  that  they  could  correspond  to  air  originating  in  the  stratosphere. 
Therefore,  the  residence  time  of  27  days  determined  in  the  summer  could  be  overes¬ 
timated  if  the  peaks  in  the  ratio  are  truly  associated  with  stratospheric  air  intrusion. 
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Figure  2.6.  Predicted  ”p/“p  activity  ratio  in  rain  as  a  function  of  the  residence  time  of 
tropospheric  aerosols. 
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Table  2.6.  Resktence  time  (t)  derived  firom  the  activity  ratio  (R)  in  each  piecipitatioa  event 
sampled.  Samples  collected  the  same  day  were  averaged,  x  is  calculated  assuming  that;  (1)  Ro  «  0.7; 
and  (2)  Ro- 0.6. 


date  R 

dpat/d^ 


Residence  time  t 
days 


Ro^.7 

range* 

Ro-0.6 

range* 

OS-Mar-91 

0.9110.14 

31 

10-57 

61 

30-160 

lO-Mar-91 

0.9910.08 

45 

30-62 

95 

60->290 

20-Mar-91 

1.08i0.08 

65 

47-95 

290 

>290 

22-Apr-9l 

1.1310.09 

83 

56-150 

290 

>290 

01*May*91 

0.9710.12 

41 

21-70 

85 

45->290 

ll-Jun-91 

0.8710.11 

25 

8-43 

50 

27-90 

20-Jun-91 

1.1010.12 

72 

43-150 

290 

>290 

27-JUII-91 

1.0210.12 

52 

28-85 

no 

57->290 

07-Jul-9l 

0.7110.11 

14 

0-32 

19 

0-37 

08-JU1-91 

0.66t0.10 

0 

0-8 

10 

0-27 

18-JU1-91 

1.0210.14 

52 

26-95 

no 

53->290 

19-JU1-91 

1.1610.16 

95 

47.>290 

290 

>290 

07-Aug-91 

0.81l0;10 

16 

1-31 

37 

18-60 

26- Aug- 91 

0.9710.11 

42 

23-65 

85 

48->290 

OS-Sep-91 

0.7010.09 

0 

0-12 

17 

1-32 

ll-Sep-91 

0.8710.10 

25 

10-42 

50 

30-85 

14.Sep-91 

1.0010.1S 

47 

22-90 

100 

45->290 

16-Sep-91 

0.8010.09 

14 

2-27 

35 

18-55 

29-Sep-91 

1.0310.13 

54 

28-95 

130 

58->290 

28-Oct-91 

O.8I1O.IO 

15 

4-27 

37 

18-60 

07.Nov.91 

0.9810.13 

43 

22-75 

90 

45->290 

17-NOV-91 

0.9210.13 

32 

12-58 

64 

32-160 

16-Jan.92 

1.0010.14 

47 

23-85 

100 

48->290 

21.Jan.92 

1.1310.12 

83 

50->290 

290 

>290 

2S.Jan-92 

i.08i0.1S 

65 

35-180 

290 

>290 

31.Jan.92 

1.0910.13 

70 

40-150 

290 

>290 

02.Fd>.92 

12010.14 

120 

60->290 

290 

>290 

14.feb.92 

O.9I1O.I2 

31 

12-53 

60 

32-130 

17.Mar.92 

0.9410.18 

37 

8-80 

70 

27->290 

26.Mar.92 

0.9810.12 

43 

18-85 

90 

48->290 

*  range:  Rq  *0  to  Ro  +o. 
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2.5 


90 

Figure  2A  Monthly  Sr  concentrations  in  pCi/1  in  precipitations  at  Bermuda  from  1967  to  1974 
(Health  and  Safety  Laboratory,  HASL-315, 1977). 


CONCLUSIONS 


1.  A  contrast  was  observed  between  Maich-June  1991/  Januaiy-March  1992 
and  July-Novonber  1991,  possibly  suggesting  a  seasonal  change  in  the  atmo^heric 
processes  at  Bermuda.  Activity  ratios  ^  ’  P/  ^  P  higher  than  tiie  annual  average  of  0.96 
were  observed  in  April  1991  and  January  and  February  1992  at  Bermuda.  Lower-than- 
average  activity  ratios  were  observed  in  September  and  October  1991.  The  contrast 
might  be  due  to  seasonal  changes  in  mixing  of  stratospheric  air  and/or  changes  in  mixing 
between  upper  and  lower  troposphere.  Injection  of  stratosph  r  :  air  due  to  tropopause 
folding  might  explain  the  high  ratio  observed  in  April  1991  at  Bermuda. 

2.  The  calculated  residence  time  of  tropospheric  aerosols  of  27±7  days  could  be 
overestimated  if  the  high  ratios  observed  in  July  to  November  1991  are  due  to  stratos¬ 
pheric  air  mixing  vrith  tropospheric  air.  The  residence  time  would  also  be  overestimated 
if  the  troposphere  is  not  well  mixed. 

3.  The  range  of  activity  ratios  in  rainwater  suggests  that  the  ’  ’  p/’  ^  p 
production-rate  ratio  is  1.06-1.24  (atom  ratio). 
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Chapter  Three 


IN  SITU  PRODUCTION  OF  COSMOGENIC 
AND  ^P  IN  SEAWATER 


I  INTRODUCTION 

Cosmic  rays  are  capable  of  inducing  nuclear  transmutations  in  the  atmosphere, 
in  the  oceans  and  in  rocks.  These  interactions  result  in  the  production  of  a  suite  of 
radionuclides.  Spallation  reactions  on  atnK>spheiic  Ar  are  believed  to  be  the  major 

32  33 

mechanisms  for  the  production  of  P  and  P  on  earth  (Lai  and  Peters,  1967).  In  the 

ocean,  it  is  believed  that  an  additional  source  derives  from  spallation  reactions  on  Q,  S 

and  K  in  seawater  (Lai  et  al.,  1988),  but  direct  measurements  of  the  production  rates  of 

both  P  and  P  are  scarce,  preventing  accurate  assessment  of  the  importance  of  this 

32  33 

process  relative  to  the  wet  deposition  of  P  and  P  produced  in  the  atmosphere. 

32  33 

The  purpose  of  this  study  was  to  investigate  the  in  situ  production  of  P  and  P 
in  seawater  relative  to  their  fallout  rates  via  wet  deposition  at  Bermuda  (geomagnetic 
latitude  X  =  43.6  ^  N).  NaCl,  KNO3 ,  MgS04  and  CCI4  targets  were  exposed  to  cosmic 
rays  at  sites  of  different  altitudes.  The  sites  chosen  were  Mount  Washington,  New 
Hampshire  (X  =  55.7  ®  N,  altitude  =  1917  m);  1’  Aiguille  du  Midi,  France  (X  =  47.4  ®  N, 
altitude  =  3840  m);  and  Woods  Hole,  Massachusetts  (X  =  53  ®  N,  altitude  =  0  m). 
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n  BACKGROUND 

Primary  cosmic  rays  interact  strongly  with  atmospheric  nuclei  giving  rise  to 
secondary  cc»mic  rays.  Pomerantz  (1971)  described  3  types  of  cosmic  ray  cmnponents: 
(1)  the  nucleonic  or  nuclear*active  ccMnponent  (neutrons  and  protons),  (2)  the  meson  or 
hard  ccxnponent  (pions  and  muons),  and  (3)  the  electromagnetic  or  soft  conqxment 
(electrons  and  ynys).  Muons  are  weakly  interactive  particles  and  do  not  generally 
contribute  significantly  to  the  production  of  cosmogenic  nuclides  in  the  atmosphere, 

39  40 

except  in  the  case  of  Cl,  for  which  muon  capture  by  Ar  is  the  major  production 
mechanism  (Winsbeig,  1956).  In  fact,  the  majority  of  the  nuclear  disintegrations  are 
induced  by  neutron  capture  or  nucleonic  spallation  reactions  (Rossi,  1952;  Lai  and  Peters, 
1967).  The  nuclear  disintegration  rates  and  the  production  rates  in  the  atmosphere  are 
approximated  by  an  exponential  increase  with  altitude.  In  the  atmosphere,  the  attenuation 

^  32  33  36 

length  of  the  production  rates  of  most  cosmogenic  nuclides,  e.g.  Be,  P,  P,  Cl  (Rama 
and  Honda,  1961;  Mabuchi  et  al„  1971;  Nakamura  et  al.,  1972;  Zreda  et  al.,  1991),  is 
identical  to  the  attenuation  length  of  neutrons  in  air  (Rossi,  1952;  Simpson  and  Fagot, 
1953;  Hess  et  al.,  1959).  The  attenuation  length  of  slow  neutrons  (neutrons  responsible 
for  the  spallation  reactions,  typically  with  E  >  40  MeV)  is  a  constant  at  a  given  latitude  in 
the  lower  troposphere  and  was  determined  to  be  145  ±  5  g/cm^  at  geomagnetic  latitudes 
of  40-55  °  N  (Simpson  and  Fagot,  1953;  Yamashita  et  al.,  1966;  Lai,  1991).  As  a  result 
of  the  neutron  dependence,  the  production  rates  of  nucleonic  spallation  products  vary 
significantly  with  latitude.  The  magnitude  increases  substantially  toward  the  poles, 
reflecting  the  higher  neutron  fluxes  at  high  latitudes  due  to  decreased  shielding  of  cosmic 
rays  at  the  poles  (Rose  et  al.,  1956;  Lai  and  Peters,  1967). 

32 

The  magnitude  of  nuclear  disintegration  rates  and  of  the  production  rates  of  P 

33 

and  P  have  been  estimated  previously  by  two  nuclear  physics  methods  (Lai,  1958;  Lai 
and  Peters,  1967)  and  by  direct  measurements  of  the  activities  produced  in  a  particular 
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target  (Lai  er  al.,  1960,  Bhandari  et  al.,  1968).  The  nuclear  physics  methods  refer  to  two 
distinct  methods.  The  first  method  is  a  theoretical  one  and  relies  on  the  knowledge  of  the 
neutron  flux  and  the  cross  section  of  interacdtm.  The  second  method  is  experimental  and 
consists  of  determining  the  nuclear  disintegration  rate  and  the  yield  of  the  spallation 

32  33 

reaction.  The  direct  measurements  of  the  production  rates  of  P  and  P  in  Ar  gas  have 
allowed  independent  verification  of  the  calculations  (Lai  et  al.,  I960;  Bhandari  et  al., 
1968).  The  latitude  and  altitude  coverage  of  the  production  rates  of  cosnaogenic  nuclides 
determined  by  direct  measurement  is  very  scarce.  In  contrast,  the  neutron  flux  has  been 
measured  extensively  with  neutron  monitors  and  is  relatively  well  known  at  all  latitudes 
and  altitudes.  The  wide  coverage  of  the  neutron  nx>nitors  has  allowed  extrapolation  of 
the  few  direct  measurements  of  production  rates  to  all  latitudes  and  altitudes  (Lai  and 
Peters,  1967).  The  nuclear  disintegration  rate  due  to  neutrons  has  been  parameterized  in 
the  form  of  polynomials  expressed  as  a  function  of  altitude  for  particular  latitudes  be¬ 
tween  0  and  90°  N  (Lai,  1991).  These  polynomials  serve  as  altitude  and  latitude  scaling 
factors  for  any  direct  measurement. 
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In  seawater,  production  of  P  and  P  can  occur  via  nucleonic  spallation  reac¬ 
tions  in  Cl,  S  and  K  due  to  the  abundance  of  these  targets.  The  in  situ  production  rates 
have  been  estimated  recently  (Lai  et  al.,  1988)  based  on  the  data  of  Rama  and  Honda 

32  32 

(1961).  The  nucleonic  spallation  in  S  inducing  P  has  a  well  known  cross  section  and 
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thus  the  production  rate  of  P  can  be  calculated  fairly  accurately.  Otherwise,  there  have 

32 

been  a  few  direct  measurements  of  the  production  rate  of  P  in  S  and  Cl  exposed  to  cos¬ 
mic  rays  (Rama  and  Honda,  1961,  Husain  and  Kuroda,  1968;  Mabuchi  et  al.,  1971)  but 
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none  at  all  in  K.  In  the  case  of  P,  there  has  been  one  measurement  of  in  situ  production 
in  Cl  by  Rama  and  Honda  (1961).  The  nucleonic  spallation  reactions  which  are  believed 
to  occur  in  each  target  are  given  in  table  3.1. 
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Although  neutrons  are  the  inaj(»‘  cosmic  ray  component  re^xMisible  for  nuclear 
disintegrations  in  the  atmo^here,  at  the  earth's  surface  and  below  the  neutron  flux  is 
sutetantially  diminished.  In  contrast,  the  muon  attenuation  length  is  an  order  mag¬ 
nitude  larger  than  that  of  neutrons  (Kurz,  1986;  Lai,  1987;  Bilokon  et  al.,  1989). 
Although  muons  are  weakly  interactive  particles,  they  can  become  the  primary  source  of 
cosmogenic  nuclides  at  depth  below  the  surface  (Rama  and  Hcmda,  1961;  Tanaka  et  al.., 
1968;  Takagi  and  Tanaka,  1968;  Mabuchi  et  al.,  1971;  Hampel  et  al.,  1975;  Lai,  1988; 
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Nishiizumi  et  al.,  1989).  For  instance,  Mabuchi  et  al.  (1971)  detected  P  in  CS2  exposed 
at  sea  level  and  underground  and  reported  that  a  significant  flaction  of  the  production  of 
at  depth  is  due  to  muons  according  to  the  reaction  ^^S(p-,  y)  ^P. 

Variations  in  the  neutron  flux  can  occur  over  time,  and  it  is  important  to  ex- 
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amine  the  time  scale  of  these  variations  with  respect  to  the  half-life  of  P  and  P.  It  is 
well  known  that  changes  in  the  production  rates  occur  due  to  changes  in  the  solar  mag¬ 
netic  activity.  It  is  well  established  that  solar  activity  modulates  the  intensity  of  cosmic 
rays,  and  by  far  the  most  pronounced  effect  is  the  1 1-year  solar  cycle  (O'Brien,  1979;  Lai 
and  Peters,  1967).  Other  shorter  time-scale  variations  occur  as  well,  like  the  reduction  in 
cosmic  ray  intensity,  i.e.  Forbush  effect  (Lockwood,  1962)  or  the  27-day  variations  (Lai 
and  Peters,  1967).  Changes  in  the  neutron  flux  associated  with  the  27-day  cycle  are 
small,  and  its  effect  on  the  production  rates  has  usually  been  neglected  (Lai  and  Peters, 
1967).  Solar  flares  can  produce  dramatic  increases  in  the  neutron  flux  by  factors  of  1  to  2 
(Lockwood,  1962),  but  this  effect  is  limited  in  time  to  a  few  hours  up  to  a  day  and  is  thus 
generally  neglected.  Similarly  the  Forbush  effect,  which  produces  a  decrease  of  the 
neutron  flux,  is  limited  in  time  to  a  few  days  and  is  also  neglected  (Lockwood,  1962). 
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Table  3.1.  Nuckv  feactioBs  occuning 
i^lpiir  Innnic  qMlIation  of  O.  S  and  K. 

V 


Tvgei 

Reaction 

”a 

”a(ii,o)’*p 

“a 

”a(n^pn)”p 

”a(ii.a2n)’^> 

’’a 

”a(n,on)”p 

’’k 

K(n.2a)  P 

39k 

39„  .  .  .  33_ 

K  (n,a2pn)  P 
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in  Experimental  procedure 


Salts  of  NaQ.  MgS04 ,  KNO3  as  well  as  CCI4  were  exposed  on  Mount 
Washington  (geomagnetic  latitude  55.7  °  N,  altitude  1917  m).  on  I'Aiguille  du  Midi  (47.4 
^  N,  3840  m)  and  at  Woods  Hole  (53  °  N,  10  m).  The  salts  were  packed  in  polyethylene 
bags  and  heat  sealed.  The  carbon  tetrachloride  was  exposed  in  2.5  liter  glass  bottles.  The 
targets  were  placed  on  die  roofs  of  buildings  at  each  site  and  away  from  obstacles  (Le., 
walls)  to  allow  a  2n  geometry.  The  targets  were  exposed  for  a  minimum  of  2  months, 
after  which  the  activities  produced  should  have  been  more  than  95%  and  83%  of  satura- 

32  33 

don  exposure  for  P  and  P  respectively.  At  Woods  Hole  and  on  1'  Aiguille  du  Midi  a 
combination  of  salts  was  exposed.  The  salts  consisted  of  NaCl,  KCl,  CaCl2  ■H2  O  and 
MgS04  mixed  in  the  proportions  as  they  occur  in  seawater  (i.e..  Cl:  19.35  g/kg,  SO4  ; 
2.71  g/kg,  K:  0.40  g/kg  and  Ca:  0.41  g/kg).  Chlorine  is  the  most  abundant  target  in 
seawater,  and  more  experiments  were  pursued  with  that  target  than  any  other.  The 
samples  were  generally  processed  within  4  days  after  the  end  of  the  exposure.  The  salts 
exposed  on  I'Aiguille  du  Midi  were  sent  by  air  to  WHOl  (flight  was  approximately  6 
hours  at  10-11  km  at  50-55  °  N).  On  Mt.  Washington,  NaCl,  KNOj ,  CCI4  and  MgS04 
were  exposed  separately  for  68  days  on  top  of  the  Observatory.  The  targets  were  brought 
quickly  to  WHOI  and  analyzed. 

No  corrections  for  initial  activities  were  made  for  the  high  altitude  exposure 
experiments,  because  exposure  to  cosmic  rays  was  long  enoug’.:  (68  and  82  days)  to 
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ensure  virtually  complete  decay  of  any  initial  activity  of  P.  The  production  rate  of  P  is 

expected  to  be  very  low  at  sea  level,  so,  in  the  small  ouantities  of  laboratory  chemicals 
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exposed  at  high  altitudes,  P  would  be  undetectable  pnor  to  exposure  (activity  below  the 
errcH*  on  the  background  of  0.01-0.02  cpm).  For  the  exfiosure  experiment  at  sea  level,  a 
blank  of  1  kg  of  NaCl  was  run.  No  activities  were  detected  in  the  blank  within  an  error  of 
0.01-0.02  cpm  associated  with  the  background.  At  most  the  initial  activity  was  0.02  cpm 
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32 

in  1  kg  of  NaCl  or  0.2  qpm  in  10  kg.  After  64  days  of  exposure  the  activity  of  P  would 
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be  4%  of  the  initial  activity  (i.e.,  0.006  cpm),  and  the  activity  of  P  would  be  17%  of  the 
initial  activity  (i.e.,  0.008  cpm).  It  is  thus  assumed  that  no  corrections  for  initial  activity 
have  to  be  made  for  any  of  the  experiments. 

After  exposure  the  salts  were  dissolved  in  distilled  water.  The  salts  exposed  at 
Woods  Hole  were  dissolved  in  a  140-liter  drum,  while  all  the  others  were  dissolved  in  20- 
liter  cubitainers.  Each  solution  was  spiked  with  KH2  PO4  .  An  FeQj  solution  and  then 
an  excess  of  ammonia  were  added  to  each  solution  to  precipitate  Fe(OH)3 .  The 
precipitate  of  ferric  hydroxide  was  allowed  to  settle.  The  supernatant  was  decanted  and 
the  precipitate  was  centrifuged,  washed,  and  finally  dissolved  in  HCl.  The  solution  was 
then  boiled  to  a  small  volume.  Nitric  acid  was  added  and  the  solution  boiled  until  HCl 
was  all  removed.  Finally  the  nitric  acid  solution  was  made  1  N.  The  rest  of  the  chemical 
separation  procedure  is  explained  in  Chapter  One. 

Phosphate  in  the  carbon  tetrachloride  phase  was  extracted  into  the  aqueous 
phase  in  a  two-liter  separatory  funnel.  The  4  I  of  CCI4  were  extracted  1  I  at  a  time.  For 
each  extraction,  1  i  of  distilled  water  was  spiked  with  KH2  PO4 .  The  two  phases  were 
shaken  vigorously  for  3  to  4  minutes  and  allowed  to  separate  for  half  an  hour.  After  each 
separation,  the  PO4  in  the  aqueous  phase  was  measured  by  the  classic  molybdenum  blue 
method  (Murphy  and  Riley,  1962).  100%  of  the  KH2  PO4  added  was  recovered  in  the 
aqueous  phase  after  each  separation.  Phosphorus  was  co-piecipitated  with  Fe(OH)3 . 
The  test  of  the  procedure  is  identical  to  the  one  described  for  the  salts. 

IV  RESULTS  AND  DISCUSSION 
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The  production  rates  of  P  and  P  were  calculated  for  each  target  at  each  site 
according  to  the  following  equation; 
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(3.1)  dC(t)/dt  =  F-XC(t) 

The  solution  of  equation  (3.1)  is: 

(3.2)  C(t)  =  F/X-(F/X-C®)exp(-Xt) 

where 
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C(t):  concentration  of  For  P  at  time  t  (atom/mol  target); 

C?  :  ctmcentration  of  or  at  the  beginning  of  exposure  (atom/mol  target); 
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F:  production  rate  of  For  P(atozn/inin/nK>l  target); 

32  33  -I 

X:  decay  constant  of  For  P(min  ); 
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It  was  assumed  that  C?  «  F  /  X  (no  initial  concentrations  of  P  or  P),  thus  equation 

(3.2)  becomes: 

(3.3)  A(t)  =  C  X  =  F  (1  -  exp  (X  t)) 

where 
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A(t):  activity  of  P  or  P  at  time  t 
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The  activities  of  P  and  P  produced  in  each  target  at  the  end  of  the  exposure  period  were 
computed  from  the  net  count  rates,  the  counting  efHciencies  and  the  chemical  yield 
(Table  A.  10).  Corrections  for  saturation  exposure  were  done  using  equation  (3.3).  The 
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production  rates  of  P  and  P  in  Cl,  S  and  K  at  each  site  are  presented  in  table  3.2.  For 
the  sample  exposed  on  1'  Aiguille  du  Midi  and  sent  by  air,  no  corrections  were  made  to 
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account  for  the  in  situ  production  of  P  and  P  during  the  6  hour  flight.  The  production 
rates  at  10  km  for  an  open  exposure  (not  shielded  by  the  aircraft  structure)  are  estimated 
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to  add  a  10%  and  7%  contribution  to  the  production  rates  at  3840  m  for  P  and  P, 
respectively.  In  reality  there  is  shielding  and  the  contribution  during  the  flight  is  some¬ 
what  smaller  than  the  calculation  for  open  exposure. 
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The  logarithm  of  the  production  rate  of  P  and  P  in  chlorine  was  plotted  as  a 
function  of  altitude  (Hgure  3.1).  The  altitudes  or  atmospheric  depths,  expressed  in 
^cm^ ,  were  calculated  using  the  polynomial  fit  of  the  altitude/pressure  curve  of  the 
Manual  of  the  ICAO  standard  atmosphere  given  by  Lai  (1991).  The  curve  fits  of  the 
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logarithm  of  the  production  rate  of  P  and  P  yield  slopes  which  allow  the  determination 
of  the  exponent  and  thus  of  the  attenuation  length.  The  production  rates  as  a  function  of 
altitude  can  be  approximamd  by  the  following  relationship: 

(3.4)  P  =  Po  exp((1033  -  z)  /  L  ) 

n 

where 

P  is  the  production  rate  at  altitude  z  (g/cm^ ). 

Po  is  the  production  rate  at  sea  level,  i-C-  zs  1033  g/cm^  . 

L^:  attenuation  length  of  the  production  rate. 

The  logarithm  of  the  production  rates  is  thus  simply  approximated  by: 

(3.5)  LnP  =  LnPo +  1033/L  -z/L 

n  n 

A  value  of  L  =  ISS  ±  15  g/cm^  is  calculated  for  the  attenuation  length  of  the  production 
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rate  of  P.  This  value  agrees  well  with  the  previously  determined  neutron  attenuation 
length  of  145  g/cm^  at  40-50  °  N  in  the  lower  troposphere  (0-5  km).  The  good  agree- 
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ment  between  the  attenuation  length  of  the  production  rate  of  P  and  the  previously 
estimated  value  of  145  g/cm^  confirms  that  neutrons  are  the  major  cosmic  rays  inducing 
in  chlorine. 

For  P  the  least-squares  fit  of  the  logarithm  of  the  production  rate  as  a  function 
of  altitude  gives  an  attenuation  length  of  195  ±  50  g/cm^ .  The  errors  on  each  individual 
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production  rate  measurement  of  P  in  Cl  are  large,  which  implies  that  the  inverse  of  the 
slq;>e  of  the  fitted  line  is  not  precisely  determined  with  only  three  measurements.  In  table 
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3.3,  the  production  rates  of  P  and  P  in  Cl  and  S  obtained  in  this  study  are  compared 
with  previous  estimates.  They  differ  form  the  results  of  Rama  and  Honda  (1961),  and 
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Husain  and  Kuroda  (1968),  but  the  production  rate  of  P  in  S  is  in  good  agreement  with 
the  estimate  of  Mabuchi  et  al.  (1971). 
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Table  3^.  Productioa  ntes  of  ^  and  in  Cl,  S  and  K  targets  at  1'  Aiguille  du  Midi,  Mount 
Washington  and  Woods  Hole. 


Sample,  locatkn. 

Target  Production  rates 

Atom 

weight,  target,  exposure  time 

lO'^atom/minAnol 

ratio 

«p  Mp 

”P/^^ 

Woods  Hole  (0,53'*  N) 

Exposure:  64  days 


Cl-WH* 

10  kg:  NaG,  KQ,  MgS04  .  CaCh 

a 

3.0  ±0.6 

1.4  ±05 

0.5  ±0.2 

Mt.  Washington  (1917m,  55.7°  N) 
Exposure:  68  days 

a-MW 

4Jkg:NaCl 

Cl 

32±7 

7±2 

0.2  ±0.1 

CU-MW 

6.4  kg:  ecu 

a 

10±1 

4±1 

0.4  ±0.1 

K-MW 

4.2kgofKN03 

K 

10±3 

n.d. 

n.d. 

S-MW 

1.9kgofMgSQ4 

s 

27±7 

n.d. 

n.d. 

Aiguille  du  Midi  (3840  m,  47*^  N) 

Exposure:  82  days 

a-MA* 

2.4  kg:  NaCL  KCL  MgS04  .CaCh 

a 

37  ±4 

10±2 

0.27  ±0.1 

*  In  these  samples  Cl  is  the  overwhelmingly  dominant  target  and  thus  the  amount  of  target  refers 
to  CL 

nxL:  Not  determined. 
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Figure  3.1.  Logarithm  of  the  production  rates  of  P  and  P  in  chlorine  target  The 
least-squares  of  the  logarithm  of  the  production  rates  are  computed. 
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Tabic  33.  CompMuon  of  published  estimates  of  the  productioD  mes  of  Pnd  PinandSatsea 
level 


Reference 

Geomagnetic 

Target 

Production  rate 

latitude 

atom/  min/  kg  target 

“P 

”P 

Rama  and  Honda  (1%1) 

44N 

a 

0.2410.03 

0.09  ±0.04 

44N 

s 

0J2±0.06 

luL 

Husain  and  Kufoda  (1968) 

49  N 

Cl 

0.76  ±0.08 

n.d. 

49  N 

s 

1.5  ±0.1 

n4. 

Mabuchi  era/.  (1971) 

25N 

s 

0.18  ±0.04 

n.(L 

(44  N)* 

s 

(0.20))* 

Thisworic 

47-55.7  N 

a 

0.07  ±0.01 

0.026  ±0.007 

55.7  N 

s 

0.16  ±0.04 

n.d. 

55.7  N 

K 

0.05  ±0.01 

n.d. 

n.d.:  Not  detennined. 

*  Results  extrapolated  to  44  ^  N  by  Mabuchi  et  al.  (1971). 
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To  assess  the  relative  importance  of  the  in  situ  production  rates  of  P  and  P  in 
the  upper  ocean  versus  the  wet  deposition  rates  of  ^^P  and  ^^P  at  Bermuda  (43.6  °  N),  the 
measured  production  rates  at  the  three  sites  were  scaled  to  the  same  altitude  (i.e.,  sea 
level)  and  to  a  geomagnetic  latitude  of  43.6  ^  N.  The  altitude  and  latitude  scaling  were 
done  using  the  polymnnials  of  the  nuclear  disintegration  rates  given  by  Lai  (1991). 
Linear  interpolations  of  the  coefficients  were  made  to  determine  the  particular  polyno¬ 
mials  for  the  latitudes  of  the  study  sites.  The  latitudinal  variations  turn  out  to  be  very 
small,  and  thus  the  latitudinal  corrections  to  scale  the  production  rates  to  43.6  °  N  were 
snudl. 
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The  production  rates  of  P  and  P  (in  atom/min/mol)  were  then  determined, 
taking  into  account  the  concentrations  of  the  target  in  seawater.  The  production  rates  of 
^and^  *P  in  the  upper  ocean  (in  atom/min/cm^ )  were  then  calculated  by  integrating  from 
0  m  to  infinity,  the  following  relationship: 

(3.6)  P  =  Po  exp(-  z  /  L  ) ,  the  integration  gives: 

z  n 

(3.7)  P*=Po  L 

n 

where  P  is  the  integrated  production  rate  in  the  upper  ocean. 

P  is  the  production  rate  at  a  depth  z  below  sea  level, 
z 

Po  is  the  production  rate  at  sea  level,  i.e.  z=  0. 

L^;  attenuation  length  of  neutron  in  seawater. 

The  attenuation  length  of  neutrons  in  seawater  below  the  surface  is  somewhat 
different  than  the  attenuation  length  above  the  surface.  Authors  have  reported  attenuation 
lengths  of  155  +  15  g/cm^  for  the  in  situ  production  of  cosmogenic  nuclides  in  rocks 
exposed  at  the  surface  (Kurz,  1986;  Lai,  1987;  Brown  et  al.,  1992).  It  is  assumed  in  the 
calculation  that  the  attenuation  length  of  neutrons  below  the  earth's  surface  is  the  same  in 
seawater  as  in  rocks.  An  average  value  of  155  g/cm^  is  chosen  for  the  calculation.  The 
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results  are  presented  in  table  3.4.  The  in  situ  production  rates  of  P  and  P  in  Cl,  K  and  S 
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are  then  calculated  by  adding  the  production  rates  of  P  and  P  measured  for  each  target. 
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Unfortunately,  tfie  production  rate  of  in  K  could  not  be  determined  because  of  the 
thickness  of  the  source.  The  in  situ  production  rates  of  P  and  P  in  the  upper  ocean  at 
bermuda  (43.6  ®  N)  are  thus  deteimined  to  be  2.2  +  0.4  10*^atoin/cm^  /min  for  ”p  and  8  ± 
3  10  attHn/cm  Anin  for  P.  In  table  3.5,  the  in  situ  oceanic  production  rates  of  P  and  P 
at  other  latitudes  are  then  computed,  using  the  polynomials  of  the  disint'^qration  rates  as 
scaling  factors  (Lai,  1991).  These  production  rates  are  a  factor  of  three  lower  than  the 
mean  in  situ  oceanic  production  rates  estimated  by  Lai  et  al.  (1988). 

The  length  of  the  salt  exposure  experiments  (July  to  September  1992)  is  ex¬ 
amined  with  respect  to  the  longer  time  record  of  the  neutron  flux  to  investigate  time 
changes  in  the  production  rates.  Figure  3.2  shows  the  neutron  count  rate  given  by  the 
neutron  monitor  at  the  Mt.  Washington  Observatory  since  1954.  It  is  observed  that  1991 
corresponds  to  a  year  of  minimum  neutron  flux  in  the  11 -year  solar  cycle.  Although  the 
major  variations  in  the  neutron  count  rate  are  due  to  the  1 1-year  cycle,  there  are  neverthe¬ 
less  variations  of  up  to  10%  on  time  scales  of  60-80  days  (see  years  1982  and  1991  in 
figure  3.2).  Therefore  the  in  situ  production  rates  of  ^^P  and  ^^P  determined  in  1992  on 
time  scales  of  64-80  days,  are  possibly  offset  from  the  mean  of  the  low  neutron  flux 
period  by  10%.  That  correction  will  be  made  when  the  neutron  data  become  available  for 
1992.  The  mean  of  the  neutron  flux  varies  by  factors  of  20-25  %  between  a  low  and  a 
high  flux  period  of  the  11 -year  cycle.  It  is  thus  expected  that  the  production  rates  of 
radiophosphorus  will  also  vary  by  20-25%  between  a  low  and  a  high  flux  period. 
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Table  3^.  In  situ  oceanic  production  rates  of  Pand  PmCl.SaiidK 
targets.  The  production  rates  are  caiculatert  for  the  upper  ocean  and 
geomaptetic  latitude  of  Bermuda  (X  ■  43.6  “n). 


Target 

Producdoo  rate 

2 

aiom/cm  /mm 

“P 

”P 

a-wH 

a 

2.4*  lO” 

1.1*10'^ 

Cl-MW 

a 

5.1*10"* 

DiQM 

CU-MW 

a 

1.6*10"* 

K-MW 

K 

2.8*10"^ 

luL 

S-MW 

s 

2.2*10'* 

djL 

a-MA 

a 

1.9*10"* 

n.d.:  Not  determined. 


Table  3^.  In  situ  oceanic  inxxluction  rates  of  ^^P  and  ^^P 


Geomagnetic  Production  rate 

latitude  atom/  cm^  /  min 


33„ 


0 

uio"* 

0.5  10"^ 

10 

1.4 10*^ 

03  10*^ 

20 

1.6 10"* 

0.6 10*^ 

30 

1.9  10*^ 

0.7  10"^ 

40 

2.1  lO"^ 

0.8  10"^ 

43.6 

2.2  10*^ 

0.8  10"^ 

50 

23  lO"^ 

0.8  10*^ 

60-90 

2310"^ 

0.8  10*^ 

90 


V  CONCLUSION 

32  33 

The  total  in  situ  production  rates  of  P  and  P  due  to  nucleonic  spallation  are 
estimated  to  be  2.5  ±  0.4  10  ^  atom/cm^  /min  few  ^P  and  8  ±  3  10  ^  atom/cm^  /min  for  ^^P 
at  Bermuda  (43.6  °  N).  These  production  rates  are  small  compared  to  the  annual  wet 
deposition  rates  of  S  10  atom/cm^/min  for  P  and  8  10  atom/cm^/min  few  P  at 
Bermuda  (5  and  1  %,  respectively). 

This  investigation  has  not  included  the  quantification  of  the  production  rates 
due  to  muons.  The  above  production  rates  are  thus  lower  limit  estimates.  It  is  known 
from  previous  studies  that  muons  can  contribute  to  the  production  rate  of  cosmogenic 
nuclides,  but  the  contribution  is  only  a  fraction  of  the  production  rate  due  to  neutrons, 
which  is  in  turn  a  small  fraction  of  the  wet  deposition  rates.  It  is  thus  believed  that  the 
muon  contribution  does  not  significantly  affect  the  conclusion  that  in  situ  production  is 
small  compared  to  wet  deposition.  There  have  been  no  studies  of  the  muon  production  of 

v>  33 

P  and  P  in  Cl,  and  that  needs  to  be  pursued  to  fully  assess  the  importance  of  that 

32  33 

process  on  the  in  situ  production  of  P  and  P. 
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Chapter  Four 


COSMOGENIC^^P  AND  RECYCLING 
IN  PLANKTON  OFF  BERMUDA 


I  INTRODUCTION 

In  spite  of  the  tremendous  potential  of  cosmogenic  radioisotopes  of  phosphorus 
for  tracing  short-term  processes  involving  phosphorus,  their  extremely  low  natural  abun¬ 
dance  in  the  marine  environment  has  impeded  their  use  in  biogeochemistry.  Early  studies 
utilizing  artificial  ^  ^  P  took  advantage  of  the  release  of  the  tracer  from  the  water  cooling 
system  of  nuclear  reactors  (Chakravarti,  1964).  In  a  later  study  ^  ^  P  was  released  in  a  salt- 
marsh  to  trace  the  P  cycle  in  that  ecosystem  (Pomeroy  et  al.,  1967).  It  is  only  recently 
that  the  first  measurements  of  naturally  produced  cosmogenic  ^  ^  P  were  obtained  in  the 
total  dissolved  phosphorus  (TDP)  pool  (Lai  et  a/.,  1988;  Lai  and  Lee,  1988,  Lee  et  al., 
1991),  in  the  dissolved  inorganic  phosphorus  (DIP)  pool  (Lee  et  al.,  1992),  and  in 
plankton  (Lai  et  al.,  1988;  Lai  and  Lee,  1988),  opening  new  avenues  for  studying  the  P 
cycle  in  the  environment  Although  it  is  very  desirable  to  utilize  both  radioisotopes,  only 
four  measurements  of  cosmogenic  P  in  seawater  have  been  reported  (Lai  and  Lee, 
1988). 
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The  difficulty  in  measuring  ^  ^  P  and  ^  ^  P  is  due  both  to  their  low  activities  in 
seawater  and  to  the  relatively  large  amount  of  stable  phosphorus,  which  hinders  detection 
of  the  radioisotopes,  particularly  the  soft  ^-emitting  ^  ^  P.  All  previously  published  data 
(Ml  ^  ^  P  in  plankton  were  collected  from  eutrophic  provinces  of  the  ocean,  and  ^  ^  P  could 
hardly  be  detected  because  of  self  absorption  in  the  source  due  to  the  large  amount  of 
stable  P  (Chapter  One).  In  that  respect  the  Sargasso  Sea,  with  its  low  DIP  and  padculate 
total  phosphorus  (PTP)  concentrations,  is  an  ideal  ocean  province  to  pursue  a  study  of 
*  ^  P  and,  especially,  ^  ^  P. 

In  this  chapter,  the  activity  ratio  ^  ^  P/^  ^  P  is  used  as  a  means  of  assessing  the 
residence  time  of  phosphorus  in  different  size  fractions  of  particulate  matter  and  trophic 
position.  The  principle  of  the  use  of  the  activity  ratio  ^  ^  P/^  ^  P  is  that  the  ratio  increases 
with  time  due  to  the  differential  radioactive  decay  rates  of  ^  ^  P  and  ^  ^  P.  An  increase  in 
the  ratio  in  plankton  relative  to  rain  or  total  dissolved  phosphorus  can  be  interpreted  as 
aging  of  the  phosphorus  in  that  plankton  pool  or  in  the  pool  from  which  it  is  derived, 
assuming  that  isotopic  fractionation  is  small.  Finally  the  in  situ  grazing  rate  of 
z(x>pIankton  is  estimated  from  the  specific  activities  of  ^  ^  P  and  ^  ^  P  in  both  suspended 
matter  and  plankton  tows. 

n  BACKGROUND 

While  there  has  been  limited  use  of  natural  cosmogenic  ^  ^  P  and  ^  ’  P,  tracer 
techniques  have  been  rq)plied  extensively  to  study  both  marine  and  freshwater  ecosys¬ 
tems.  Uptake  rates  by  phytoplankton  and  turnover  times  of  P  in  plankton  have  been 
consequently  determined  in  the  marine  system  (Perry,  1976;  Harrison  et  al.,  1911  \  Perry 
and  Eppley,  1981;  Herbland,  1984;  Sorokin,  1985;  Harrison  and  Harris,  1986).  The 
studies  showed  substantial  recycling  of  phosphorus  during  the  short  incubation  periods 
suggesting  an  underestimation  of  uptake  rates  by  these  methods  (Harrison,  1983; 
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Harrison  and  Harris,  1986).  In  addition,  tracer  techniques  have  been  applied  to  the  study 
of  exchange  rates  between  DIP,  dissolved  organic  phosphorus  (DOP)  and  PTP  in 
seawater  (Watt  and  Hayes,  1963).  Assimilation  efficiencies  of  P  of  single  species  of 
zoq)lankton  or  natural  assemblages  have  been  detennined  using  ^^P  in  laboratotry 
experiments.  Some  of  the  studies  have  been  reviewed  by  Comer  and  Davies  (1971). 

The  Sargasso  Sea  is  an  oligotrophic  oceanic  province.  Typically  DIP  con¬ 
centrations  in  the  upper  ocean,  as  determined  by  the  standard  colorimetric  method 
(Murphy  and  Riley,  1962;  Koroleff,  1983),  are  below  20-30  nM  (Boyle  et  al.,  1986;  Knap 
et  cd.,  1991  and  1992;  Michaels  et  cd.,  1992).  With  the  recent  in:q)rovement  of  analytical 
techniques,  it  has  been  possible  to  measure  DIP  concentrations  as  low  as  1-2  nM  at  sta¬ 
tion  S  (32°  10'  N;  64°  30*  W)  (Ormaza-Gonzalez,  1990).  On  average  DIP 
concentrations  are  about  10  nM  while  at  station  ALOHA  in  the  central  Pacific  gyre,  they 
are  typically  about  100  nM  in  the  upper  ocean  (Karl  and  Tien,  1992).  Particulate  total 
phosphorus  (PTP)  is  also  low,  ranging  from  a  few  nM  to  60  nM  in  the  Sargasso  Sea  and 
is  typically  around  10  nM  (Bishop  et  al.,  1986a;  Onmaza-Gonzalez,  1990;  J.  Ammerman, 
personal  communication  at  BATS  workshop  1992).  Particulate  nitrogen  (PN)  is  also  very 
low,  ranging  up  to  0.5  pM  (Altabet,  1989;  Knap  et  ai,  1992).  While  DIP  and  PTP  pools 
have  similar  sizes,  DOP  is  by  far  the  major  pool  of  P  at  BATS  station  (31°  50'N,  64°  10' 
W)  and  at  station  S,  with  typical  concentrations  of  l(X)-250  nM  (Ormaza-Gonzalez,  1990; 
J.  Ammerman,  personal  communication  at  BATS  workshop  1992). 

Despite  the  importance  of  P  in  the  ocean,  the  study  of  the  P  cycle  has  been 
neglected  in  favor  of  C  and  N,  partly  because  of  the  difficulties  of  conventional  tracer 
techniques,  measurement  and  sampling  techniques,  and  the  perception  by  biologists  that 
N  is  the  limiting  nutrient  in  the  ocean  (Eppley  et  al.,  1973;  Eppley  and  Peterson,  1979; 
Goldman  et  ai,  1979;  Perry  and  Eppley,  1981;  McCarthy  and  Carpenter,  1983). 
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Phosphorus  is  an  essential  element  for  living  organisms  and  is  regarded  by  geochemists 
as  the  ultimate  nutrient  limiting  productivity  in  the  ocean  because  of  nitrogen  fixation 
(Redfield,  1958;  Redfield  et  al.,  1963;  Broecker  and  Peng,  1982;  Smith  1984;  Smith  and 
Atkinson,  1984;  Fanning,  1989).  In  addition,  very  few  studies  have  succeeded  in  quan¬ 
tifying  the  flux  of  phosphorus  to  the  deep  ocean  due  to  large  uncertainties  in  the  C/P  ratio 
in  sinking  particulate  matter.  The  observed  ratios,  i.e  150-900,  are  found  to  be  generally 
higher  than  the  predicted  ratio  of  106/1  by  Redfield  or  of  127/1  predicted  by  Peng  and 
Broecker  (Bogdanov  and  Shaposhnikova,  1971;  Bishop  et  al.,  1977;  Knauer  et  al.,  1979; 
Bishop  et  al.,  1980;  Knauer  and  Martin,  1981;  Bishop  et  al.,  1986b;  Martin  et  al.,  1987; 
Liebezeit,  1991).  Phosphorus  recycling  in  the  upper  ocean  appears  to  be  nwre  efficient 
than  cycling  of  carbon  and  nitrogen.  It  is  also  suggested  by  the  high  release  of  P  relative 
to  N  during  remineralization  of  both  particles  (Vaccaro,  1963;  Grill  and  Richards,  1964; 
Harrison,  1980)  and  organic  matter  (Jackson  and  Williams,  1985;  Smith  et  al.,  1985; 
Smith  et  al.,  1986;  Qrrett  and  Karl,  1987).  In  addition,  P  is  released  rapidly  from 
biogenic  particles  upon  collection  (Collier  and  Edmond,  1984;  Cbllier,  1991).  Overall  it 
is  not  clear  whether  the  high  C/P  ratios  measured  in  sinking  particles  are  meaningful 
features,  i.e.  faster  recycling  of  P  than  N  or  C,  or  artifacts  of  collection,  i.e.  loss  due  to 
cell  lysis. 

m  RESULTS  AND  DISCUSSION 

A  seasonal  study  of  the  distributions  of  ^  ^  P  and  ^  ^  P  in  different  particulate 
fractions  was  pursued  in  the  mixed  layer  at  the  U.S.  JGOFS  Bermuda  Atlantic  Time- 
series  Study  (BATS)  site  (31°  50’N,  64°  lO'W),  80  km  southeast  of  Bermuda,  during 
1991  and  1992. 

The  specific  activities  of  ^  ^  P  {i.e.,  ’  ^  P/P)  and  ^  ^  P  (i.e.,  ’  ^  P/P)  were  deter¬ 
mined  in  surface  suspended  matter  in  July  and  November  1991  and  in  February  and  April 
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1992  (Figure  4.1)  and  in  suspended  maner  collected  at  three  depths  in  April  1992 
(Figures  4.2  and  4.3).  The  activity  ratios  ’  ^  P/*  ^  P  were  derived  for  surface  and  depth- 
integrated  plankton  tows,  for  surface  suspended  matter  (Figures  4.4a  to  4.4e  and  4.5)  and 
fOT  the  depth  profile  of  ^  ’  P^  ^  P  in  suspended  matter  in  April  1992  (Figure  4.6).  The 
samples  were  collected  during  6  cruises  on  R/V  Weatherbird  n  and  R/V  Endeavor  at  the 
BATS  site.  The  results  for  4-5  March,  29  April-1  May,  15-19  July  (BATS  34)  and  1 1-15 
November  (BATS  38)  of  1991  and  for  24-26  February  (BATS  41A)  and  12-27  April  (EN 
235)  of  1992  are  presented  in  tables  4.1-4.6.  The  samples  collected  in  March  and  May  do 
not  correspond  to  BATS  cruises,  and  hydrographic  data  from  the  closest  BATS  (BATS 
29A  and  BATS  31)  cruises  are  presented  in  figures  4.7  and  4.8.  The  hydrographic  data 
for  BATS  34,  38  and  41A  arc  presented  in  figures  4.9, 4.10  and  4.1 1.  For  the  Endeavor 
cruise  (EN  235)  T  and  NOs  profiles  are  presented  in  figtues  4.12a  and  4.12b  for  April 
25. 

Suspended  particulate  matter  was  typically  collected  by  filtering  2,(X)0  to  7,000 
liters  of  water  through  Millipore  filters  (3  or  8  pm  mesh  size).  Presumably,  particles 
smaller  than  3  or  8  pm  were  collected  as  well  because  clogging  occurred  rapidly  during 
filtration.  The  Millipore  filters  were  changed  as  soon  as  the  flow  rate  started  to  decrease 
about  every  hour.  Large  particles  were  harvested  with  plankton  nets  (67,  150,  300  and 
500  pm  mesh  size).  These  large  particles  contained  macrozooplankton  and  also  large 
phytoplanktonic  species,  such  as  cyanobacteria  (Jrichodesmium).  An  attempt  was  made 
to  determine  ^  ^  P  and  ^  ^  P  on  large  particles  collected  by  filtering  5,000  1  of  water 
through  a  53-pm  Nitex  filter  in  February  1992.  The  amount  of  P  on  the  Nitex  filter  was 
very  low  (i.e.,  0.3  nM)  preventing  analysis  of  ^  ^  P  and  ^  ^  P.  In  April  92,  comparison  was 
made  between  PTP  for  a  surface  sample  collected  on  GF/F  (0.7  pm)  and  on  Millipore  (3 
pm  mesh  size).  It  was  found  that  20  %  of  the  suspended  particles  >0.7  pm  were  in  the 
fraction  0.7-3  pm.  This  percentage  is  a  lower  limit  estimate,  because  particles  smaller 
than  3  pm  were  probably  collected  on  the  3-pm  filter  due  to  clogging. 
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Figures  4.4^  b.  Activity  ratio  of  P/  P  in  suspended  pamcles  and  plankton  to' 
(67, 3(X)  pm  mesh  size)  during  May  and  July  1991  at  BATS  station 


November  1991 


1  I  Suspended  particles  (>8  \im) 

J5S3  Plankton  tow  67  um 
mu  Plankton  tow  300  um 

rx/yi  Plankton  tow  67  um,  depth-integrated  (0-60m) 

POOQ  Plankton  tow  300  um,  depth-integrated  (0-60  m) 

Figure  4.4.c.  Activi^  ratio  of  in  suspended  particles  and  plankton  tows  (67, 
3(X)  pm  mesh  size)  during  November  1991  at  BATS  station 
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LxWl  Plankton  tow  67  um 


Suspended  particles  (>3  um) 


Plankton  tow  500  um  f'X/Vl  Plankton  tow  150  um 


Figures  4.4.dy  e.  Activity  ratio  of  ^P/^P  in  suspended  particles  and  plankton  tows  (67, 
150, 5(X)  pm-mesh  size)  during  d)  February  and  e)  April  1992  at  BATS  station. 
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Figure  4^.  Monthly  activity  ratio  in  rainwater  front  March  1991  to  March  1992,  and  activity  ratio  in  surface 
suspended  particles  and  plankton  tows  for  May,  July  and  November  1991,  and  February  and  April  1992. 
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Figure  4.6.  Profile  of  the  activity  ratio  **p/*P  in  suspended  particles  (>3  pm)  during 
EN  235, 25-26  April  1992,  at  BATS  station. 
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Both  ^  ^  P  and  ^  ^  P  were  measured  in  suspended  particles  and  plankton  tows 
for  the  first  time.  Our  ability  to  determine  both  isotopes,  as  explained  before,  is  due 
mainly  to  the  very  low  stable  P  concentrations  in  the  Sargasso  Sea.  As  expected,  the 
q)ecific  activities  in  the  particulate  matter  are  one  or  two  orders  of  magnitude  highn  than 
those  encountered  in  the  South  California  Bight,  in  the  Celtic  Sea  and  in  the  Bedford 
Basin  (Lai  et  al.,  1988;  Lai  and  Lee,  1988). 

The  major  features  of  the  data  set  that  will  be  discussed  and  interpreted  are:  (1) 
the  specific  activities  in  suspended  particles  are  highest  in  November  1991;  (2)  the  ratio 
3  3  py3  2  p  in  piaijrton  tows,  which  is  either  equal  to  or  significantly  higher  than  the  ratio 
in  rain,  is  highest  in  February  1992;  and  (3)  the  phosphorus  radioisotopes  were  not 
homogeneously  distributed  in  zooplankton  in  March  1991,  allowing  a  labile  and  a  refrac¬ 
tory  pool  to  be  identified.  In  addition,  the  residence  time  of  P  in  macrozooplankton  and 
the  in  situ  grazing  rate  and  defecation  rate  of  zooplankton  are  estimated. 

32  33 

!•  Specific  activities  of  P  and  P 

The  specific  activities  range  from  1  to  3  dpm/mg  P  for  ^  ^  P  and  ^  ^  P  in  the 
suspended  matter  and,  from  0.3  to  3  dpm/mg  P  for  ^  ^  P  and  0.9  to  5  dpm/mg  P  for  ^  ^  P  in 
plankton  tows  (Tables  4.1  to  4.6).  The  activities  in  suspended  particles  were  usually  0.3- 
0.9  dpm/m^  for  ’  ^  P  and  0.3- 1.3  dpm/m^  for  ’  ^  P  (Table  4.7).  PTP  concentrations 
(collected  on  the  3  or  8-pm  filter)  were  found  to  range  from  7  to  19  nM  and  average  10 
nM  (Table  4.7).  This  is  in  agreement  with  the  few  existing  measurements  of  PTP  (filters 
of  0.45  pm  mesh  size)  at  Hydrostation  S  and  at  BATS  (Qrmaza-Gonzalez,  1990;  J. 
Ammerman,  personal  communication  at  BATS  workshop  1992)  suggesting  that  the 
Millipore  3  and  8- pm  filters  collected  particles  smaller  than  their  nominal  pore  sizes.  The 
highest  specific  activities  were  obtained  in  November  1991,  with  values  ranging  from  2 
to  5  dpm/mg  P  for  ’  ^  P  and  3  to  6  dpm/mg  P  for  ^  ^  P  (Table  4.4).  More  typical  values 
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for  the  period  sampled  are  around  0.2-1  dpm/mg  P  (Tables  4.2, 4.3, 4.5  and  4.6).  In  April 
1992,  three-point  depth  proHles  of  ^  ^  P  and  ^  ^  P  in  suspended  particulate  matter  were 
obtained,  and  both  show  similar  exponential  decrease  with  depth  (Figures  4.2  and  4.3). 
The  specific  activities  observed  at  BATS,  in  suspended  matter  and  plankton  tows,  are  one 
(H*  two  orders  of  magiutude  higher  than  the  speciHc  activities  of  ^  ^  P  in  plankton  tows 
reported  previously  (Lai  et  al.  1988;  Lai  and  Lee  1988).  The  authors  observed  specific 
activities  of  ^  ^  P  ranging  from  0.02  to  0.2  dpm/mg  P  in  the  Bedford  Basin  and  from  0.01 
to  0.1  dpm/mg  P  in  the  California  Bight.  The  higher  specific  activities  observed  in  the 
Sargasso  Sea  are  due  to  the  extremely  low  DIP  concentrations,  typically  below  the  detec¬ 
tion  limit  of  the  molybdenum  blue  method,  i.e.,  20-30  nM  and  as  low  as  1-2  nM  (Ormaza- 
Gonzalez,  1990).  In  contrast,  the  DIP  concentrations  in  the  California  Bight  are  0.3  to  2 
pM  (Lee  et  al.,  1991). 

Off  Bermuda,  the  major  features  of  the  upper  250  meters  are;  (1)  the  break¬ 
down  of  the  seasonal  thermocline  in  January/February;  (2)  the  productivity  peak  in 
March/mid-April  following  the  deep  winter  mixing  event;  and  (3)  the  development  of  the 
seasonal  thermocline  fiom  April  to  the  fall  (Schroeder  and  Stommel,  1969;  Menzel  and 
Ryther,  1960,  1961;  Michaels  et  al.,  1992).  During  the  period  of  development  of  the 
seasonal  thermocline,  new  production  is  usually  very  low.  In  the  absence  of  removal 
processes  (vertical  mixing  at  the  base  of  the  mixed  layer  and  net  export  of  particulate 
matter  out  of  the  mixed  layer),  the  activities  of  ^  ^  P  and  ^  ’  P  in  TDP  build  up  in  the 
mixed  layer.  At  steady  state  the  activities  can  be  simply  determined  from  the  balance 
between  tiie  supply  and  the  decay  of  the  inventories  of  ’  ^  P  and  ^  ’  P.  Accumulations  of 
^  ‘  °  Pb  and  trace  metals  in  the  mixed  layer  have  been  observed  in  the  summer  off 
Bermuda  (Boyle  et  al.,  1986;  Jickells  et  al.,  1990).  The  high  specific  activities,  i.e.,  2  to  4 
dpm/mg  P  for  ^  ^  P  and  3  to  6  dpm/mgP  for  ^  ^  P,  encountered  in  suspended  particles  and 
plankton  tows  in  mid-November  might  be  in  part  due  to  the  accumulation  of  ^  ^  P  and 


109 


’  ’  P  in  the  mixed  layer.  In  July  1991,  February  1992  (BATS  34,  41A)  and  April  1992 
(EN  235),  the  specific  activities  of  ^  ^  P  and  ^  ^  P  in  suspended  particles  in  the  mixed  layer 
were  lower,  because  they  were  controlled  more  by  vertical  mixing  and/or  new  production 
than  by  radioactive  decay.  In  February  1992  (BATS  41A).  the  mixed  layer  depth  was 
210  m  (Figure  4.1 1  a)  and  PO4  was  detectable  at  40  m  and  below  by  the  nxrlybdenum 
blue  method  (Figure  4.1 1  b).  These  conditions  imply  both  supply  of  stable  P  into  the 
phosphate-depleted  upper  ocean  and  loss  of  radiophosphorus  from  the  uppo*  ocean, 
leading  thus  to  lowest  specific  activities.  In  April  1992  (EN  235)  the  specific  activities  of 
both  ^  *  P  and  ^  ^  P  were  highest  at  the  surface  and  decreased  with  depth  (Figures  4.2, 
4.3).  April  1992  corresponded  to  the  onset  of  stratification,  with  the  seasonal  thermocline 
forming  and  being  destroyed  on  time-scales  of  a  few  days  during  the  12-day  long  cruise 
EN  235.  The  profrles  of  the  specifre  activities  of  ^  ^  P  and  ^  ^  P  reflect  these  phenomena 
integrated  over  a  few  weeks  and  thus  reflect  the  decrease  with  depth  observed  in  the  time- 
average  T  profiles  during  the  EN  235  cruise. 

2-  Activity  ratio  ^^P/^^P  in  plankton 

In  this  section  qualitative  interpretations  of  the  activity  ratio  will  be 

given,  while  in  sections  4  and  5  simple  models  will  be  developed  to  give  quantitative 
interpretations.  The  activity  ratios  ^  ^  p/  ^  p  measured  in  suspended  particles  and  in 
plankton  tows  are  found  always  to  be  either  equal  to  or  higher  than  the  ratio  in  rain, 
which  averages  0.96  on  an  annual  basis  (Figure  4.5).  This  is  explained  by  an  increase 
with  time  of  the  ratio  ^^p/^^p  due  to  the  differential  decay  rates  of  ’  ^  P  (14.3  d)  and 
^  ^  P  (25.3  d).  The  ratio  is  expected  to  increase  in  higher  trophic  levels  relative  to  the 
ratio  in  rain  because  the  doubling  time  of  plankton,  to  a  first  approximation,  is  positively 
correlated  with  the  particle  size  and  ranges  from  hours  for  phytoplankton  to  hundreds  of 
days  for  certain  species  of  zooplankton  (Sheldon  etaL,  1972;  Tranter,  1976). 
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The  turnover  time  of  P  deduced  from  the  time-dependent  curve  of  the  ratio 
3  3p^2p  zooplankton  is  model-dependent  Two  end-member  nuxlels  will  be 
developed  for  zooplankton  in  section  four.  Ihe  first  model  predicts  that  the  ratio  in¬ 
creases  with  time  until  the  activities  of  ^  ^  P  and  ^  ^  P  reach  a  steady  state,  at  which  time 
the  curve  reaches  a  plateau.  This  end-member  nnodel  represents  a  situation  where 
zooplankton  behave  as  an  "open  system",  continuously  grazing  on  phytoplankton  and 
regenerating  a  fraction  of  that  food.  The  other  end-member  model  represents  the  case 
where  zooplankton  is  behaving  as  a  "closed  system",  neither  grazing  or  regenarating  P. 
This  second  model  predicts  that  the  ratio  will  increase  exponentially  with  time. 

In  May  1991  the  activity  ratio  ^  ^  P^  ^  P  in  plankton  tow  (>300  pm)  averaged 
1.4.  This  ratio  is  significantly  higher  than  the  average  ratio  of  1.02+0.03  in  rain.  In  July 
1991  the  activity  ratio  in  surface  suspended  particles  (>3  pm;  presumed  to  be  mostly 
phytoplankton)  was  0.75+0.2  (Figure  4.5).  That  is  lower  than  the  activity  ratio  of  1. 5+0.3 
in  the  plankton  tow  (>67  pm;  presumed  to  be  mostly  macrozooplankton).  This  observa¬ 
tion  indicates  a  measurable  aging  of  P  in  macrozooplankton  relative  to  phytoplankton.  In 
contrast,  the  activity  ratio  in  suspended  particles  is,  within  the  errors,  identical  to  the  ratio 
in  rain.  The  similarity  of  the  ratios  suggests  that  suspended  particles,  presumably 
phytoplanktonic  in  origin,  have  a  very  short  P  turnover  time,  estimated  to  be  less  than  a 
few  days,  which  is  consistent  with  the  doubling  time  of  phytoplankton,  i.e.,  hours  to  days 
(Sheldon  et  al.,  1972).  Furthermore,  Harrison  and  Harris  (1986)  estimated  a  turnover 
time  of  P  in  phytoplankton  on  the  order  of  1-2  days  in  the  Sargasso  Sea.  As  pointed  out 
by  the  authors,  the  tracer  technique  used  for  the  estimation  tends  to  underestimate  the 
uptake  rates,  and  it  is  likely  that  the  turnover  time  of  P  is  even  shorter  than  1-2  days. 

In  November  1991,  there  was  no  difference  in  the  activity  ratios  of  1.6  between 
the  surface  suspended  particles  (>8  pm)  and  the  surface  plankton  tows  collected  with 
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either  a  67  or  300  (im  mesh  net  (Figure  4.4c).  The  ratio  was  also  1.6  in  the  depth- 
integrated  (0  to  60  m)  plankton  tows  collected  with  nets  of  mesh  size  67  pm  and  3(X)  pm. 
The  ratio  of  1.6  in  the  surface  suspended  particulate  sample  can  be  explained  in  part  by 
the  accumulation  of  ^  ^  P  and  ^  ^  P  in  the  shallow  mixed  layer  from  June  to  the  fall.  In 
other  words,  ^  ^  P  and  ^  ^  P  are  being  supplied  in  a  ratio  of  0.9-1.0  from  June  to  November 
(ISO  days)  and  are  being  removed  mostly  by  decay  from  the  noixed  layer.  Similar  conclu- 
sicHis  woe  derived  for  ^  Pb  at  station  S  where  accumulation  of  ^  °  Pb  was  found 
throughout  the  June-September  period  (Boyle  et  eU.,  1986).  Similarly,  some  trace  ele¬ 
ments  accumulate  in  the  shallow  mixed  layer  during  the  summer  months  (Jickells  et  aJ., 
1990).  In  November,  Trichodesmiwn  (blue-green  algae)  was  found  in  abundance  in  two 
tows  and  was  present  in  all  the  tows  (>67  pm  and  >300  pm).  Trichodesmium, 
presumably,  has  a  very  short  turnover  time  of  P  like  other  phytoplankton,  and  their 
presence  might  account  for  the  similarity  of  the  ratios  in  plankton  tows  and  in  suspended 
matter. 

In  February  1992  (BATS  41A)  the  activity  ratio  was  0.9+0.3  in  suspended 
particles  (>3  pm)  and  was  on  average  3.3±0.5  in  plankton  tows  (>67  pm  and  >500  pm), 
as  seen  in  figure  4.S.  There  were  no  significant  differences  between  the  ratios  in  the  67- 
pm  plankton  tow  and  the  ratio  in  the  S(X)-pm  plankton  tow  (Figure  4.4d).  In  an  attempt 
to  separate  species  of  zooplankton,  material  from  a  large  tow  was  poured  into  a  container 
filled  with  surface  seawater,  and  the  fast  sinking  fraction  was  separated  from  the  rest  of 
the  tow  which  consists  of  mostly  radiolarians.  The  samples  were  referred  to  as  ZSOO-Bl 
and  ZS00-B2  (fast  sinking  large  particles)  and  21S00-T  (rest  of  the  tow)  (Table  4.5). 
Those  two  samples  did  not  exhibit  different  ratios,  suggesting  a  homogeneous  residence 
time  of  P  in  the  different  species  and  types  of  zooplankton.  The  most  striking  feature  of 
the  February  cruise,  compared  to  the  other  cruises,  was  the  great  difference  in  the  activity 
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ratio  in  the  two  size-fractions  (Figure  4.S).  This  observation  suggests  that  macio- 
zooplanktcm  was  quite  old,  with  low  turnover  times,  relative  to  phytoplankton  and  that 
part  or  all  of  the  P  in  zooplankton  was  not  being  exchanged  with  external  sources  of  P. 
This  can  be  explained  if  the  grazing  activity  of  the  animals  had  partially  shut  down  during 
the  winter,  leading  to  considerable  decay  of  ^  ^  P  and  ^  ^  P  already  present  in  zot^lankton. 
It  would  explain  the  high  ratios  of  3  or  4  measured  in  macrozoc^lankum  while  the 
primary  producers  exhibit  ratios  of  only  0.9.  It  is  possible  that  the  plankton  collected  in 
February  1992  in  67  pm  and  5(X)  pm  mesh  size  nets  were  large  interzonal  migrants  living 
in  the  mesopelagic  zone  (Longhurst  and  Harrison,  1988). 

In  April  1992  an  opposite  trend  was  found.  The  activity  ratio  in  the  suspended 
particles  (>  3  pm)  averaged  1.74:0.2,  while  the  plankton  tows  (>  ISO  pm)  exhibited  ratios 
of  1.()4^.2  on  average.  No  significant  difference  was  found  between  the  activity  ratio  in 
night  and  day  plankton  tows  Z150-N  and  Z150-D  (Table  4.6),  and  both  had  ratios  lower 
than  suspended  particles.  It  suggests  that  possibly  a  fraction  of  the  suspended  particle 
pool  consisted  at  that  time  of  detritus  with  a  slow  turnover  time  of  P.  Detritus  could 
possibly  originate  from  old  wintertime  zooplankton  producing  particles  via  excretion  or 
ejection  of  fecal  matter.  If  this  is  the  case,  a  simple  calculation  with  two  end-member 
mixing  terms  (February  and  April  terms)  shows  that  a  fraction  of  1  to  10  %  of  wintertime 
zooplankton  organic  matter  characterized  by  activity  ratios  ’  ’  P/*  ^  P  of  3-4  released  in 
the  suspended  particles  could  increase  the  ratio  in  the  suspended  pool  to  1.7  in  April. 
Then  the  question  is  why  the  ratio  is  as  low  as  1  in  the  plank  on  tows,  1  being,  within  the 
error  of  the  measurements  the  ratio  in  rain.  The  ratio  in  plankton  tows  was  1.0,  suggest¬ 
ing  very  rapid  turnover  times  of  P  on  time  scales  of  less  than  a  day.  In  April  1992,  BATS 
42B  (11  April)  showed  maximum  Chi  a,  indicating  the  occurrence  of  the  spring  bloom. 
Peak  concentrations  in  the  pigments  (fucoxanthine,  diadinoxanthine  and  diatoxanthine) 
were  observed  during  that  cruise  as  well,  possibly  indicating  the  presence  of  abundant 
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diatoms  during  the  bloom.  It  is  hypothesized  that  large  diatoms  might  have  been  present 
in  the  surface  tow  of  25  April  (EN  235),  giving  rise  to  the  obsoved  low  ^  ^  p/^  ^  p  ratios. 

3-  L4ibile  and  refractory  P  in  zooplankton 

The  zoq)lankton  analysed  in  March  1991  revealed  that  radioisotopes  of  P  were 
not  homogeneously  distributed  in  zooplankton,  and  two  operationally-defined  P  compart¬ 
ments  were  identified  (Table  4.1).  Freshly  collected  plankton  tows  were  left  in  surface 
seawato'  at  4  **  C  fw  a  few  days  in  the  dark.  A  large  amount  of  phosphate  (about  50  %  of 
die  total  amount  of  P)  was  released  into  solution  during  storage  and  was  detectable  with 
the  standard  molybdenum  blue  method  (Murphy  and  Riley,  1962).  The  amount  released 
was  far  greater  than  typical  phosphate  concentrations  in  the  surface  waters  (i.e.,  <20-30 
nM).  The  three  samples  were  filtered  through  GF/F,  and  the  filtrates  were  combined. 
Phosphate  was  extracted  on  alumina  and  processed  as  explained  in  Chapter  One.  The 
analysis  of  the  filtrate  revealed  that  the  activity  ratio  ^  ^  ^  P  in  the  combined  three 

samples  is  1.7±0.3.  That  ratio  is  significantly  higher  than  the  ratio  in  the  rain.  Two  of  the 
fractions  remaining  after  leaching  had  very  high  ratios.  One  sample  exhibits  a  ratio  as 
high  as  18,  indicating  a  very  long  residence  time  of  P  in  the  refractory  pool.  The  errors 
<Mi  the  determination  of  the  ratio  in  die  3  refiactory  fractions  were  large  because  of  the 
extremely  low  count  rate  of  ^  ^  P  in  those  samples  (almost  undetectable  in  Z67-3).  It  is 
thus  clear  that  those  samples  were  particularily  old,  and  even  if  the  precision  of  the  ratio 
determinaticMi  is  not  good,  the  slope  of  the  total  net  count  rate  is  a  qualitative  indication  of 
a  voy  high  fraction  of  ^  ^  P  in  the  sample.  The  experiment  was  repeated  in  February 
1992  and  April  1992.  In  February,  unfortunately,  the  filtrate  was  lost  during  processing. 
The  refractory  fraction  exhibited  an  activity  ratio  of  3.7±|0.8  which  was  not  very  different 
from  the  ratios  observed  in  the  bulk  of  the  other  tows,  suggesting  that  the  filtrate  or  labile 
fraction  had  a  ratio  similar  to  2-4.  This  result  indicates  that  zooplankton  had  slow  P 
turnover  times.  In  April  1992,  both  the  filtrate  and  the  refractory  fractions  were  analysed 
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(Table  4.6).  The  filtrate  was  coprecipitated  with  Fe(OH)3  and  further  puriHed  with  the 
procedure  described  in  Chapter  One.  The  anK>unts  of  stable  PO4  in  each  subsatr^le 
Z150-R  and  Z150-L  were  identical  within  the  2%  analytical  error  of  the  colorimetric 
method,  indicating  that  the  two  subsamples  are  of  equal  size  in  terms  of  stable  P.  The 
activity  ratios  ^  ^  ^  p  of  both  samples  were,  within  the  1  o  error,  identical  to  the 

average  activity  ratio  in  rain,  suggesting  very  fast  turnover  of  P  in  the  plankton  present  in 
die  tows.  This  result  can  be  explained  by  the  probable  presence  of  diatoms  in  the  tows 
during  the  spring  bloom  (April  1992). 

The  difference  in  the  ratios  in  the  operationally  deHned  labile  and  refractory 
fractions  observed  in  March  1991  suggests  that  zooplankton  encompassed  two  compart¬ 
ments  of  phosphorus,  one  labile  and  one  refractory.  Two  distinct  pools  of  P  have  been 
found  previously,  characterized  as  labile  and  stable,  with  very  different  P  turnover  times, 
i.e.  0.4  days  and  13  days,  respectively  (Conover  et  ai,  1961).  Also  it  has  been  shown 
that  zooplankton  have  the  ability  to  store  P  as  fat  in  the  males  and  in  the  reproductive 
syst^  in  the  fenudes  (Comer  and  Davies,  1971).  The  ratio  of  1.7i0-3  observed  in  the 
labile  pool  suggests  that  P  turns  over  on  a  time  scale  of  a  few  days  to  a  week.  That  pool 
of  P  is  probably  the  labile  pool  of  P  identified  in  earlier  plankton  studies  which  showed 
substantial  and  rapid  release  of  P  subsequent  to  the  collection  of  plankton  tows  (Collier, 
1984).  It  also  appears  that  P  is  excreted  efficiently  by  zooplankton  relative  to  nitrogen, 
while  nitrogen  is  enriched  in  fecal  pellets  and  thus  in  sinking  particles  (Howard,  1988), 
suggesting  that  there  is  a  pool  of  labile  P. 

The  operationally-denned  refractory  pool  of  P  might  contain  (1)  phospholipids, 
which  are  ubiquitous  in  the  membranes  of  organisms  (Wefer  et  al.,  1982;  Cembella  et  al., 
1984);  (2)  phosphonates  (Cembella  and  Antia,  1986;  Liebezeit,  1991),  or  (3)  nucleic 
acids  (Holm-Hansen,  1966).  The  labile  pool  is  likely  to  be  partly  composed  of  (1)  phos¬ 
phate  monoesters,  like  AMP,  ADP  and  ATP,  which  are  very  important  P  compounds  in 
the  energy  system  of  cells;  (2)  PO4  groups  which  arc  end-products  of  many  enzymatic 
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itactims  (Westheimer,  1987);  and  (3)  PO4  present  in  the  guts  of  the  animals. 
Qtaracteiization  of  the  P  compounds  in  a  single  species  of  zooplankton  has  been  pursued 
by  Comer  (1973),  and  the  P  conqwunds  found  were  an  acid-soluble  fraction  (50%), 
phoqiholipids,  nucleic  acids  and  phosphoproteins  (phosphate  ester). 

In  April  1992,  no  significant  difference  was  found  between  the  easily  leachable 
fitacticMi  and  the  remaining,  supposedly  more  refractory,  fraction.  Both  fractions  had  low 
3  3  py3  2  p  suggesting  rapid  turnover  of  P.  April  corresponded  to  the  spring  bloom 
period  with  high  Oil  a  in  the  euphotic  zone.  The  ratio  in  zooplankton  suggests  that 
zooplankton  was  regenerating  P  at  high  rates.  It  is  consistent  with  Comer  and  Davis's 
(1971)  study,  which  showed  a  direct  link  between  food  availability  and  increased  excre¬ 
tion  of  N  and  P  by  zooplankton. 

4-  Residence  time  of  P  and  grazing  rates  of  macrozooplankton 

In  this  section,  I  quantify  the  residence  time  of  P  in  macrozooplankton,  based 
on  the  activity  ratios  ^  ^  p/^  ^  p  measured  in  suspended  particles  and  plankton  tows  when 
the  tows  were  dominated  by  zooplankton.  As  a  consequence,  the  grazing  rate  can  be 
derived  from  the  specific  activities  of  ^  ^  P  and  ^  ^  P  in  phyto-  and  zooplankton.  The 
calculation  of  the  residence  time  is  model-dependent.  The  simplest  model  that  is  ex¬ 
amined  is  a  steady-state  model  describing  the  mass  balances  of  ^^P  and  ^^P  in 
zooplankton.  The  mass  balance  for  ^  ^  P  in  zooplankton  is  given  by  the  following  equa¬ 
tion: 

(4.1)  g(”P)  =X32(’*P)  +se^P) 

p  z  z  z 

whoe  (^  ^  P)  :  concentration  of  ^  ^  P  in  phytoplankton  (dpm/m^ ); 

P 

C  ^  P)  :  concentration  of  ^  ^  P  in  zooplankton  (dpm/m’ ); 
z 

X32  :  decay  constant  of  P  (s  ); 

g:  zooplankton  grazing  rate  (s  '); 
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r.  regeneration  rate  of  P  from  zooplankton  (s  ’); 

s:  removal  rate  of  P  by  zooplanktonic  sinking  particles  (fecal  pellets  and  mat¬ 
ter,  molts,  etc...)t  s 

A  similar  equation  is  written  for  ^  ^  P.  It  is  assumed  in  equation  (4.1)  that  zooplankton 

acquire  ^  ^  P  and  ^  ^  P  from  the  grazing  of  phytoplankton  and  assimilate  the  isotopes  in 

the  proportions  occurring  in  phytoplaidaon.  Thus  the  activity  ratio  ’  ’  P/*  ^  P  in 

zooplankton  is  the  following: 

(4.2)  R  =  R  (X3  2  +  r  +  s)  /  (X3  3  +  r  +  s) 
z  p 

Where 


R  =  (^  ^  P)  /C  ^  P)  :  activity  ratio  ’  ’  P/®  ^  P  in  zooplankton 
z  z  z 

R  =  (®  ®  P)  /(®  ®  P)  :  activity  ratio  ®  ®  P/®  ^  P  in  phytoplankton 
P  P  P 

From  equation  (4.2)  one  can  calculate  r  +  s  as  follows: 

(4.3)  r  +  s  =  (X32 -X33  *R/R  )/(R/R -1) 

z  p  z  p 

This  model  predicts  that  the  ratio  in  zooplankton  ranges  from  Rp  to  a  maximum  value  of 

1.77  R  according  to  the  magnitude  of  r  and  s  relative  to  X3  2  and  X3  3 .  The  residence 
P 

time  of  P  in  zooplankton,  t,  is  then  simply  derived  as  follows: 

(4.4)  1  =  1  /  (r  +  s) 

From  equation  (4.1)  one  can  calculate  the  grazing  rate  as  follows: 


(4.5)  g  =  (^^P)  /(”P)  (X32  +l/t)  and, 

z  p 

(4.6)  g  =  (®®P)  /(®®P)  (X33  +l/t) 

z  p 

In  July,  the  average  activity  ratio  in  zooplankton  was  1.31  ±  0.19.  I  did  not 
include  the  high  value  of  1.79  ±  0.44  in  die  calculation  of  the  average  because  that  high 
ratio  was  measured  less  precisely  than  the  two  others  (Table  4.3).  The  measured  activity 
ratio  of  0.75±0-23  in  the  suspended  particles  was  found  to  be,  within  the  error,  similar  to 
the  ratio  in  the  rain  (Figure  4.5).  The  average  activity  ratio  in  rain  for  the  three  months 
that  preceded  the  July  1991  cruise  (BATS  34)  was  0.92±0.04.  Therefore  a  value  of  0.92, 
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was  chosen  for  the  activity  ratio  in  suspended  particles.  A  residence  time  of  t  =  45  ^ 
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days  is  estimated  for  the  plankton  tow  in  July  1991,  assumed  to  be  iiK>stly  con^)osed  of 
macrozooplankton.  This  residence  time  is  similar  to  the  estimated  residence  times  of  P  in 
zoq)lankton  of  30-40  days  or  more  in  the  California  Bight  cU.,  1988).  The  es¬ 

timated  X  is  also  within  the  range  of  turnover  times  of  single  species  of  zooplankton 
which  have  been  estimated  previously  based  on  the  ratio  of  the  production  rate  to  the 
biomass.  The  results  of  some  of  these  studies  are  summarized  by  Tranter  (1976)  and 
Sheldon  et  al.  (1972).  The  turnover  time  of  zooplankton  ranges  from  a  few  days  to 
hundreds  of  days  and  is  to  a  first  approximation  positively  correlated  with  the  diameter  of 
the  particles  (Sheldon  et  al.,  1972).  Large  boreal  species  exhibited  turnover  times  of 
about  100  days,  while  temperate  smaller  species  had  turnover  times  of  about  10  days. 

In  suspended  particles,  the  specific  activities  were  2.99  ±  0.50  dpm/mg  P  for 
^  *  P  and  and  2.24  +  0.52  dpm/mg  P  for  ^  ^  P.  In  macrozooplankton,  the  specific  activities 
were  0.90  ±  0.07  dpm/mg  P  for  ^  ^  P  and  1.21  +  0.14  dpm/mg  P  for  ’  ^  P.  Grazing  rates 
of  0.021  ±  0.006  d"*for  ’  ^  P  and  0.027  +  0.009  d  ’for  ^  ^  P  were  determined  utilizing 
equations  (4.5)  and  (4.6);  or  an  average  grazing  rate  of  0.024  ±  0.006  d  *.  Particulate 
organic  carbon  (POC)  concentrations  yield  a  value  of  30  pg  C/l  in  the  150-m  euphotic 
zone,  for  BATS  34  (15-19  July  1991).  An  integrated  value  of  4,500  mg  C/m^  was 
calculated  for  the  euphotic  zone  for  that  cruise.  The  flux  of  carbon  grazed  by  macro¬ 
zooplankton  was  thus  estimated  to  be  110  ±  30  mg  C/m^/d.  In  July  1991  primary 
productivity  was  500  mg  CVm^  /d  at  BATS  station  (Michaels  et  al.,  1992);  thus  the  es¬ 
timated  grazing  flux  corresponded  to  20  ±  5  %  of  primary  productivity.  This  result  is 
within  the  range  of  the  percent  of  primary  productivity  grazed  by  a  copepod  community 
CMi  the  outer  shelf  off  Long  Island  (Dagg  and  Turner,  1982)  and  on  the  outer  and  mid 
shelf  of  the  Bering  Sea  (Dagg  et  al.,  1982). 

Assuming  that  zooplankton  have  assimilation  efficiencies  ranging  from  30  to 
70  %  (Angel,  1984;  Dagg  et  al.,  1982),  I  estimate  a  flux  of  30  to  80  mg  C/m^  /d  out  of  the 
euphotic  zone  due  to  defecation  of  zooplankton  in  July  1991.  Although  assimilation 
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efficiencies  of  zooplankton  vary  with  the  constituent  and  the  food  type  (Comer  and 
Davies,  1971),  70%  is  often  considered  to  be  a  reasonable  assimilation  efficiency  (Angel 
1984;  Dagg  et  al.,  1982).  In  that  case,  I  estimate  that  30  ±  10  mg  CJvc?  /d  leaves  the 
euphodc  zone  by  defecation  of  zooplankton.  This  second  estimate  is  in  very  good  agree¬ 
ment  with  the  flux  of  particulate  carbon  of  35  mg  CJtc^  /d,  determined  by  sediment  traps 
deployed  at  ISO  m  at  BATS  station  in  July  1991  (Michaels  et  al.,  1992).  If  the  sediment 
tn^  flux  is  accurate  my  results  suggest  that  fecal  pellets  or  fecal  maner  constitute  a  major 
component  of  the  particulate  flux  in  the  upper  ocean  at  BATS.  It  has  been  recognized 
that  fecal  material  is  a  major  source  of  particulate  matter  in  the  ocean  (Pomeroy  and 
Diebel,  1980;  Fowler  and  Knauer,  1986;  Small  et  ai,  1987),  and  some  studies  have 
shown  that  it  made  up  the  bulk  of  the  particulate  flux  in  certain  areas  like  the  Panama 
Basin  and  the  equatorial  Atlantic  Ocean  (Bishop  et  al.,  1977;  Bishop  et  al.,  1986b). 

In  May  1991  the  ratio  ^  ^  P/  ^  P  in  plankton  tows  was  1.4  +  0.1.  The  ratio  in 
the  suspended  particles  was  not  measured.  It  is  suggested  hrom  figure  3.5  that  the  ratio  in 
suspended  particles  must  range  from  a  minimum  value  of  1.0,  the  ratio  in  rain,  to  a  maxi¬ 
mum  of  1.4,  the  ratio  in  the  plankton  tows.  If  a  ratio  of  1.02±0.03  is  assumed  (the 
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average  ratio  in  the  rain  in  the  spring),  a  residence  time  of  40  days  is  estimated  for  May 
1991.  This  estimate  is  in  agreement  with  previous  estimates  made  in  the  California  Bight 
(Lai  et  al.,  1988)  and  consistent  with  turnover  times  of  various  species  of  zooplankton 
(Sheldtm  et  al.,  1972  and  Tranter,  1976). 

In  November  1991  the  ratios  in  both  suspended  matter  and  zooplankton  were 
the  same.  Therefore  a  lower  estimate  of  a  few  days  was  estimated  for  the  residence  time 
of  P  in  plankton.  In  November  1991,  Trichodesmium  was  very  abundant  in  two  tows 
(i.e.,  Z67-10  and  Z300-10)  and  was  present  in  all  the  tows.  The  presence  of 
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Trichodesmium  in  the  tows  might  explain  the  similarity  of  the  ratios  between  the 
suspended  matter  and  the  plankton  tows  in  November  1991. 

In  February  1992,  the  activity  ratio  in  suspended  matter  was  0.93+0.3  and,  on 
average,  about  3.3+0.5  in  plankton  tows.  The  previous  end-member  model  can  not  ex¬ 
plain  the  high  ratios  observed  in  plankton  tows.  It  can  only  account  for  ratios  in 
zooplankton  ranging  from  0.9- 1.0  to  1.5- 1.8.  An  extreme  scenario  is  presented  in  which 
zooplankton  acquire  P  eariy  in  their  growth  cycle  and  store  it  Radiophosphorus  would 
decay  as  zooplankton,  assumed  to  be  a  "closed  system",  become  older.  The  ratio  can  be 
predicted  by  the  following  relationship: 

(4.7)  R  =R  exp((X3  2 
z  p 

The  residence  time  given  by  equation  (4.7)  represents  the  age  of  zooplankton.  Equation 
(4.7)  yields  a  residence  time  of  P  in  zooplankton  of  60  ±  12  days  for  Februaiy  1992.  It 
means  that  for  two  months  prior  to  the  cruise  (December  and  January)  radiophosphorus  in 
zooplankton  was  decaying  witiiout  significant  exchange  with  the  environment.  In  winter, 
when  food  is  scarce,  feeding  activities  are  likely  to  be  reduced,  and  zooplankton  are  likely 
to  retain  phosphorus.  This  is  consistent  with  the  evidence  of  a  direct  link  between  food 
availability  and  excretion  (Comer  and  Davies,  1971).  The  residence  time  of  60  days  is 
probably  an  upper  limit,  because  it  is  unlikely  that  zooplankton  were  completely  isolated 
from  the  environment  for  two  months. 
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5>  Mass  balance  of  P  and  P  in  the  euphotic  zone 

In  this  section,  the  inventories  of  ^  ^  P  and  ^  ’  P  in  the  suspended  particles  are 
determined.  The  decay  of  the  inventories  in  suspended  matter  is  then  compared  to  the 
wet  deposition  rates.  For  most  '‘f  the  cruises,  except  April  1992,  the  depth  profiles  of  the 
activities  of  ^  ^  P  and  ^  ^  P  in  suspended  matter  were  not  obtained.  The  profiles  of  ^  ^  P 


and  ^  ^  P  in  suspended  particles  for  April  1992  (Figure  4.2  and  4.3)  were  fitted  with  em¬ 
pirical  exponential  functions  given  as  follows: 

(4.8)  ’  *  P  =  1.00  exp(-  z  /  67.7) 

(4.9)  ”  P  =  1.68  exp(-  z  /  74.2). 

Where 

C  ^  P)_^*  1.00  dpm/m’  and  (”  P)  =  1.68  dpm/m^ 

Each  expcMiential  was  then  integrated  from  the  surface  to  150  m,  the  depth  of 
the  euphoric  zone,  to  calculate  the  activities  in  the  euphoric  zone.  It  is  observed  that  the 
depth  profiles  of  the  specific  acriviries  of  ^  ^  P  and  ^  ^  P  in  suspended  particles  have 
decreasing  trends  with  depth  which  are  very  similar  to  one  another  and  similar  to  the  T 
profile  (Figure  4.12  a).  For  July  1991  and  November  1992  (BATS  34  and  38),  I  assumed 
tha*  ^  ^  P  and  ^  ^  P  activities  in  suspended  matter  at  any  depth  in  the  mixed  layer  were 
identical  to  the  activities  at  the  surface.  I  also  assumed  that  the  activities  of  ^  ^  P  and  ^  ^  P 
in  suspended  matter,  below  the  mixed  layer,  decreased  exponentially.  In  February  1992 
(BATS  41A)  there  was  intense  vertical  mixing  within  the  upper  250  m  of  the  ocean  and 
constant  activities  of  ^  ^  P  and  ^  ^  P  in  the  mixed  layer  (210  m)  could  not  be  assumed.  No 
inventory  computation  was  pursued  for  that  cruise.  The  inventories  of  ’  ^  P  and  ^  ^  P  in 
the  euphoric  zone  are  reported  in  table  4.8  for  July  1991,  November  1991,  and  April 
1992. 

The  inventories  of  ^  ^  P  and  ^  ^  P  in  suspended  particles  are  then  compared  to 
the  amount  of  ^  ^  P  and  ^  ’  P  delivered  by  precipitation’.  The  wet  deposition  rates  of  ^  ^  P 
and  ^  ’  P  can  depart  significantly  from  the  annual  average  rates  by  factors  of  up  to  three. 
The  deposition  rates  were  computed  from  the  average  activities  and  rainfall  rates  over  a 
poiod  of  three  months  preceding  the  cruises.  Three  months,  or  four  times  the  half-life  of 
die  longer-lived  radiophosphorus  isotope,  corresponded  to  the  period  prior  to  which  the 
activities  of  ^  ^  P  and  ^  ^  P  in  rain  could  not  remain  detectable  in  particulate  matter  at  the 
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tune  of  sampling.  Since  May  1991  was  such  a  dry  nKMUh,  it  was  assumed  as  a  first  ap¬ 
proximation  that  events  prior  to  June  1991  would  not  be  detectable  in  mid-July  1991 
(BATS  34)  and  were  thus  not  taken  into  account  in  the  calculation.  For  April  1992  (EN 
235)  the  wet  deposition  rates  used  are  the  annual  averages  since  no  rain  data  were  col¬ 
lected  after  March  26,  1992  (one  month  befcue  the  sampling  time).  The  wet  deposition 
rates  of  ^  ^  P  and  ^  ^  P  so  computed  are  listed  in  table  4.8. 

The  differences  between  the  wet  deposition  rates  and  the  decay  of  the  inven¬ 
tories  of  ^  ^  P  and  ^  ^  P  are  computed  (Table  4.8).  It  is  found  that  the  decays  of  the 
inventories  of  ^  ^  P  and  ^  ^  P  in  suspended  particles  are  in  balance  with  the  wet  deposition 
rates  of  ’  *  P  and  ^  ^  P  in  July  1991,  November  1991,  and  April  1992  in  the  euphotic  zone. 
This  result  suggests  that,  (1)  removal  of  ^  ^  P  and  ^  ^  P  onto  sinking  panicles  (new 
production)  and  removal  due  to  vertical  mixing  at  the  base  of  the  euphotic  zone  must  be 
very  small,  and  (2)  most  of  the  inventories  of  ^  ^  P  and  ^  ^  P  are  in  the  suspended  par¬ 
ticles.  The  balance  between  the  deposition  rates  of  ’  ^  P  and  ’  ^  P  and  the  decays  of  the 
inventories  indicates  that  the  residence  time  of  P  in  the  euphotic  zone  must  be  longer  than 
4  or  5  half-lives  of  ^  ^  P  and  ^  ^  P  or  about  100  days.  The  estimated  lower  limit  for  the 
residence  time  is  consistent  with  the  range  of  0.7  to  2  years  estimated  by  Oimaza- 
Gonzalez  (1990). 

The  results  indicate  that  ^  ^  P  and  ’  ^  P  inventories  must  be  very  low  in  the 
dissolved  pools  of  P,  i.e.  DIP  and  DOP.  These  findings  contrast  with  previous  studies  in 
the  Southern  California  Bight  showing  that  most  of  the  inventory  of  ’  ^  P  is  in  the  TDP 
(Lai  and  Lee,  1988;  Lai  et  al.,  1988).  In  the  Southern  California  Bight,  DIP  is  usually  the 
major  pool  of  dissolved  P  with  a  range  of  concentartions  of  0.1  to  2  pM  while  DOP 
ranges  from  0.1  to  0.8  pM  (Lee  et  al.,  1991;  Lee  et  al,  1992).  These  values  of  DOP  are 
consistent  with  the  DOP  range  in  the  Pacific  ocean,  i.e.  0.1 -0.4  pM  (Oirett  and  Karl, 
1987;  Jackson  and  William,  1985;  Ridal  and  Moore,  1992).  In  contrast,  at  Hydrostation 
S,  DIP  is  a  very  small  pool  with  concentrations  ranging  from  below  1  nM  to  30  nM, 
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while  DOP  is  the  major  dissolved  pool  of  phosphorus  with  100-250  nM  P  (Ortnaza- 
Gonzalez,  1990).  The  DOP  values  found  off  Bermuda  are  also  consistent  with  the  values 
found  by  Ridal  and  Mooie  (1990)  in  a  Gulf  Stream  station.  If  most  of  the  inventories  of 
^  ^  P  and  ^  ^  P  are  in  PTP  at  BATS  station,  as  indicated  by  the  mass  balance,  then  it  sug¬ 
gests  that  the  residence  time  of  P  in  DOP  is  long  relative  to  the  half-lives  of  ^  ^  P  and 
^  ^  P.  A  lower  limit  for  the  residence  time  of  P  is  estimated  to  be  four  or  five  tinoes  the 
half-life  of  ^  ^  P,  i.e.,  1(X)-120  days.  These  estimates  agree  with  previous  studies  which 
showed  that  only  a  small  fraction  of  the  ^  ^  P  activity  was  found  in  DOP,  i.e.,  <0-0.09 
dpm/m^  (Lee  et  al.,  1992),  leading  the  authors  to  conclude  that  the  mean  time  period  of 
DOP  utilization  is  much  longer  than  6  weeks  in  the  Southern  California  Bight 

Cosmogenic  ^  ^  P  and  ^  ^  P  studies  give  a  somewhat  different  picture  of  DOP 
availability  than  previous  studies,  because  they  allow  computation  of  rates  and  turnover 
times.  There  is  evidence  that  DOP  is  a  labile  pool,  more  so  than  DON  (Jackson  and 
Williams,  1985;  Smith  et  al.,  1986,  Smith  et  al.,  1*^85;  Orrett  and  Karl,  1987)  and  that  all 
DOP  is  available  to  the  microbiological  community  in  the  subarctic  Pacific  (Ridal  and 
Moore,  1992).  In  addition,  it  has  been  well  documented  that  marine  bacteria  and 
phytoplankton  have  evolved  cell-surface  enzymes  that  are  capable  of  cleaving  organic  P 
compounds  (Perry,  1972;  Azam  and  Hodson,  1977;  Cembella  et  al.,  1984;  Ammerman 
and  Azam,  1985;  Rivkin  and  Swift,  1985).  One  of  these  enzymes,  alkaline  phosphatase, 
has  been  detected  at  Hydrostation  S  in  activities  high  enough  possibly  to  hydrolyse  all  the 
DOP  present  in  the  euphotic  zone  (Ormaza-Gonzalez,  1990).  The  important  question 
regarding  the  role  of  DOP  in  contributing  to  phytoplankton  growth  and  new  production  is 
the  rate  of  turnover  of  DOP.  Ckismogenic  ^  ^  P  and  ^  ’  P  provide  insight  regarding  that 
question  on  time  scales  of  days  to  weeks. 
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TaUe  4.7.  Activities  of  and  concentrations  of  P  in  suspended  natter  (>  3-8  Jim) 


at  BATS  station. 


Sample 

Cruise 

Initial 

amount  of 

Pin  the 

sample 

mgp 

Volume 

of  water 

pumped 

t 

Stable  P 

suspended 

maH^r 

nM 

dpm/m^ 

”P 

dpm/m^ 

Ml 

BATS  34 

0.64ij0.01 

2.230±30 

9.2±p.3 

0.86 

0.64 

M-Sl 

BATS  38 

1.08±002 

3.97Q+10 

8.8±0.2 

035 

0.89 

M-S2 

BATS  38 

0.78ij0.01 

3.04Q±40 

8.3±0.2 

0.62 

1.13 

M2 

BATS41A 

1.48ii0.03 

5.100!±25 

9.4±D.2 

0.29 

0.27 

M-0 

EN23S 

2.93±0.06 

7.400±150 

12.8+03 

0.81 
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M-35 

EN23S 

3.01±p.06 

5.36()±40 

18.1±Q.9 

0.60 

1.18 

M-70 

EN23S 

322±p.06 

5.52Qt80 

19±1 

0.44 

0.77 

32  33 

Table  4.8.  Activities  (A),  wet  deposition  rales  (W)  and  radioactive  decay  (\A)  of  P  and  P,  and  dif- 

32  33 

ferences  between  the  wet  deposition  rates  and  the  decay  of  Pand  P  (i.e.,  W-X.A)  in  thie  euphotic  zone. 


Cruise 

XA 

W-Xa 

”a  ”a 

“w 

33w 

”P 

dpm/in^ 

dpmAn^  /d 

dpm/m^  /d 

dpm/m^  /d 

BATS  34 

45.3 

36.4 

1.7±i0.3 

1.4i03 

23 

1.0 

0 

0 

BATS  38 

47.8 

719 

23±0.5 

23±0.4 

2.3 

2.0 

0 

0 

EN235 

613 

108 

2.4±q.4 

2.1t0.4 

3.0 

3.0 

0 

0 
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m  CONCLUSIONS 


The  specific  activities  ’  *  P  and  ’  ^  P  and  the  activity  ratio  ^  ^  P/®  ^  P  in  rain¬ 
water,  suspended  particles  and  plankton  tows  allowed  determination  of  the  in  situ 
turnover  time  of  P  in  macrozooplankton  and  the  in  situ  grazing  rates.  The  turnover  time 
of  P  in  macrozooplankton  was  found  to  range  from  40  days  to  60  days  at  BATS  station. 
In  some  cruises,  the  plankton  tows  were  presumably  dominated  by  large  phytoplanktonic 
species  as  evidenced  by  the  low  ratios.  The  estimated  turnover  times  of  P  in  plankton 
tows  ate  consistent  with  the  estimates  of  doubling  times  of  species  of  zooplankton. 

An  in  situ  grazing  rate  of  0.024  +  0.006  d  '  was  determined,  and  it  was  es¬ 
timated  that  30  ±  10  mg  C/m^  /d  left  the  euphoric  zone  by  zooplankton  defecation  in  July 
1991  at  BATS.  The  estimated  particulate  flux  matchs  the  bulk  of  the  particulate  carbon 
flux  determined  by  sediment  traps  deployed  at  150  m  at  BATS,  suggesting  that 
zooplankton  fecal  matter  made  up  the  bulk  of  the  particulate  flux  at  BATS. 

The  decay  of  the  inventories  of  ®  ^  P  and  ®  ®  P  in  suspended  particles  was  found 
to  balance  the  input  of  ®  ^  P  and  ®  ®  P  by  precipitation,  suggesting  that  turnover  of  DOP, 
which  is  the  major  pool  of  dissolved  P,  is  long  relative  to  the  half-lives  of  ®  ^  P  and  ®  ®  P. 
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Figure  4.12 


Table  4.1.  Specific  activities  of^^Paad^^Pand  aoivity  ratio  ^  ^  P/’  ^  P  in  plankton  tows  coUecied 
with  a  67-^  planktoo  net  in  Match  1991. 


Sanqtk 

Date 

DqNh 

Type 

Amount  of 
stable  P 
mgP 

Specific 

activity 

of^»P 

dpm/mgP 

Specific 

activity 

of’’P 

dpmMgP 

Ratio 

3  3  2  p 

dpm/dpm 

Lf 

3^1 

sivface 

tow 

/ 

/ 

/ 

1.7±03 

267-Rl* 

3M/91 

surface 

tow 

0.99±0.02 

130±0.35 

52±1.2 

4.Qtl.2 

267-R2* 

3/4/91 

surface 

low 

1.09iP.02 

1.7ii0.9 

1.12i0.30 

0.67i0.19 

Z67-R3* 

3/4/91 

surface 

tow 

0.9Q±0.02 

03Q±0.6S 

7.5±1.1 

18.7±12.5 

Z67 

3/31/91 

surface 

tow 

1.66±j0.03 

0.81±0.09 

0.49ii0.08 

O.6I1OI6 

labik  ficaction  of  the  three  tows  combined.  The  initial  amount  of  phosphate  is  not  kr»wn,  thus 
making  it  impossible  to  detennine  the  yield  in  that  sample. 

*:  refractory  fractions. 


Table  4J1.  Specific  activities  of  ^  ^  Pand^  ^  P  aitd  activity  ratio  ^  ^  P/^  ^  Pin  plankton  tows 
collected  with  a  300>iim  mesh  net  in  29  April- 1  May.  The  samples  ate  uiplicates  from  the  same  tow. 


Sample 

Date 

Depth 

Type 

Amount  of  Specific 
stable  P  activity 
mgP  of’^P 

dpm/mgP 

Specific 

activity 

Of”P 

Ratio 

33p^2p 

dpm/dpm 

Z300-2 

5/1/91 

surface 

night  tow 

2.57±q.05 

0.79±|0.07 

1.05i0.13 

133±D.17 

Z300-3 

5/1/91 

surface 

night  tow 

4.74±Q.09 

0.71ii0.07 

l.lQf0.11 

1.55±0.22 

Z30(M 

5/1/91 

surafce 

n^ttow 

6.01±0.12 

0.67±0.07 

0.90±0.12 

1.34110.19 

132 


Table  4.3.  Specific  activities  of  ^  ^  P  and  ^  ^  P  ami  activity  ratio  ^  ^  P^  ^  P  in  suspended  particles  (>  8 
^m)  and  in  plankton  tows  (67-|Xm  mesh  size  net)  collected  during  BATS  34, 15-19  July  1991. 


Sample 

Date 

Depth 

Type 

Amount  of 
stable  P 
mgP 

Specific 

activity 

of’*P 

dpm/mgP 

Specific 

activity 

of”P 

dpmAngP 

Ratio 

3  3  2  p 

c^Mn/dpm 

Ml 

7/16/91 

surface 

pump 

0.64iQ.01 

2.991030 

2.241032 

0.7510.23 

Z67-10* 

7/18/91 

surface 

night  tow 

3.1010.06 

0.961011 

lJ2i0.21 

1.37±0.30 

Z67-11* 

7/18/91 

sutfime 

night  tow 

3.1010.06 

0.87±0.13 

1.09±0.22 

1.2510.25 

Z67-12* 

7/18/91 

surface 

night  tow 

1.2710.03 

0.89±0.12 

1J9±0.30 

1.79i0.44 

Sg-1+ 

7/16/91 

surface 

macroalgae 

3.2210.06 

1.6810.24 

2J810.42 

13510.37 

Sg-2+ 

7/16/91 

surface 

macroalgae 

2.1410.04 

267±0J5 

3J310.42 

1.3210.22 

*:  triplicates  of  the  same  tow. 

duplicates  of  a  sample  of  sargassum. 


Table  4.4.  Specific  activities  of’*Pand^’Pand  activity  ratio  ’  ’  P/*  *  P  in  plankton  lows  and 
particulate  (>  3  (im)  samples  collected  during  BATS  38, 11-15  November  1991.  Z67  and 
Z300  refer  to  plankton  tows  collected  with  a  67-iim  and  300-iim  plankton  net  respectively.  Each  tow 
is  split  in  two  samples  except  Z67-10.  Z300-10  contained  Trichodesmium  and  marine  snow  mostly. 


Sample 

Date 

Depth 

Type 

Amount  of 

Specific 

Specific 

3  3  2  p 

stable  P 

activity 

activity 

mgP 

of’^P 

of’’P 

dpm/dpm 

dpm/mgP 

dpm/mgP 

M-Sl 

11/12/91 

surface 

pomp 

1.08±0.02 

2.0410.19 

3371032 

1.6010.26 

M-S2 

11/14/91 

surface 

pump 

0.78±0.01 

2.421037 

4.411.0 

13410.90 

Z67-10 

11/12/91 

surface 

night  low 

1.24±0.02 

3.841038 

4.7110.67 

13210.21 

267-20* 

11/14/91 

0-60m 

day  tow 

2.84±p.06 

2.19±038 

3.561036 

1.621039 

267-21* 

11/14/91 

0-60m 

day  tow 

1.95±0.04 

3.58±031 

5.6710.88 

1381036 

267-30 

11/14/91 

surface 

night  tow 

2.08±0.04 

3.35±032 

4.6910.67 

1.4010.32 

267-31 

11/14/91 

surface 

night  tow 

3.12±0.06 

2.001031 

3.1210.41 

1.56±0.22 

2300-10 

11/12/91 

surface 

night  tow 

3.24±0.06 

2.751038 

4.0010.65 

1.46l0.20 

2300-11 

11/12/91 

surface 

night  tow 

2.16±0.04 

23510.43 

53010.76 

1.8210.25 

2300-20*  11A4/91 

0-60m 

day  tow 

3.15±0.06 

2.4810.37 

4.3010.65 

1.7310.22 

2300-21*  11/14/91 

0-60ta 

day  tow 

3.15±D.06 

2.77±0.43 

4.95i0.73 

1.79i032 

•:  tows  collected  by  continually  yo-yoing  the  net  on  the  hydrowirc  from  0  to  the  base  of  the  mixed 
layer  at  60  m. 
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Table  4.5.  Specific  acnvities  of  ^  ^  P  and  ^  ^  P  and  activity  ratio  ^  ^  p/  ^  P  in  plankton  and 
suspended  panicles  (>  3  lim)  in  samples  collected  during  BATS  41A.  24-2S  February  1992.  ZSOO  and 
Z67  refer  to  samples  collected  with  plankton  nets  of  mesh  size  S00-|tm  and  67-iim  respectively.  One  of 
the  Z500  was  separated  into  two  fractions  21S00-T  and  2^00-B.  ZSOO-T  is  the  fraction  that  stayed 
suspended  in  seawater.  It  contained  mostly  radiolarians.  ZSOO-B  is  the  fraction  that  sank  and  contained 
mostly  marine  snow  and  cnistacians  (ZSOO-Bl  and  ZS00-B2  are  duplicates).  ZSOO-R  is  the  refracUMy 
fraction  of  a  tow,  the  labile  fraction  was  lost  during  processing. 


Sample 

Date 

Depth 

Type 

Amount  of 
stable  P 
mgP 

Specific 
activity 
of'^P 
dpm/mg  P 

Specific 
activity 
Of”P 
dpm/mg  P 

33p^Zp 

dpm/dpm 

M2 

02/25/92 

surface 

pump 

1.48±0.03 

1.01±0.19 

0.94±0.23 

0.93±0.31 

Z500-R 

02/25/92 

surface 

night  tow 

1.77±0.O4 

0.34±0.05 

1.25±0.23 

3.66±0.W 

ZSOO-Bl 

02/25/92 

surface 

night  tow 

3.49±0.07 

0.39±0.04 

1.16±0.18 

2.95±0.55 

Z500-B2 

02/25/92 

surface 

night  tow 

3.49±0.07 

0.34±0.05 

1.36±0.23 

4.0±1.0 

Z500-T 

02/25/92 

surface 

night  tow 

2.16±0.04 

0.49±0.09 

0.99±0.25 

2.03±0.51 

267-1 

02/25/92 

surface 

night  tow 

2.04±0.04 

0.34±0.06 

1.41±0.38 

4.1±1.3 

Z67-2 

02/25/92 

surface 

day  tow 

1.37i0.03 

0.21+0.08 

1.36+0.65 

6.5±4.8 

Table  4.6.  Specific  activities  of  ^  *  P  and  ^  ^  P  and  activity  ratio  ^  ^  P/*  ^  P  in  suspended  particles 
(>3  lim)  and  in  plankton  tows  collected  during  EN  235  in  April  1992  at  BATS  station.  ZISO  refer  to  a 
tow  collected  with  a  150  flm  net  Z150-L  and  R  refer  to  the  labile  and  refiractory  fractions. 

Sample 

Date 

DqKh 

Type 

Initial 
amount  of 
stable  P 
mgP 

Specific 

activity 

of’^P 

dpmAngP 

Specific 
activity 
Of”P 
dpm/mg  P 

Ratio 

33p/32p 

dpnt/t^Mn 

M-0 

4/25/92 

surface 

pump 

193±0.06 

2.05±020 

339±0.45 

1.6610.20 

M-35 

4/25/92 

35  m 

pump 

3.01±0.06 

1.07i0.31 

2.1 110.79 

1.9710.81 

M-70 

4/25/92 

70  m 

pump 

3.22±0.06 

0.75i0.12 

U2l0.25 

1.76±0.30 

Z150-N1 

4/26/92 

surface 

night  tow 

3.8710.07 

1.2510.26 

1.0910.23 

0.8710.17 

Zr50-Dl 

4/26/92 

surface 

day  tow 

3.88±0.07 

1.37i0.18 

1.66±0.21 

I.2I1O.I6 

Z150-L 

4/26/92 

surface 

night  tow 

4.06±0.08 

0.8410.08 

0.9310.11 

1.1010.28 

Z150-R 

4/26/92 

surface 

night  tow 

3.33±0.06 

0.72±0.08 

0.5810.08 

0.8210.13 

I 
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CONCLUSIONS 


I  PRESENT  FINDINGS 

32  33  -3 

The  wet  deposition  rates  of  P  and  P  at  Bermuda  were  detennined  to  be  5  10 
atom/cm^  /min  for  and  8  10  ^  atom/cm^  Anin  for  ^^P.  The  activity  ratio  ’  ’  P^  ^  P 
delivered  in  rainwater  was  found  to  be  close  to  0.96  throughout  the  year  (March  1991  to 
March  1992).  The  study  of  ’  ’  P  and  ^  *  P  in  rainwater  showed  that  the  ratio  ’  ’  P/*  ^  P  is  a 
powerful  tracer  of  atmospheric  processes.  The  variations  in  the  ratio  ’  ’  P/*  ^  P  can  be 
utilized  to  determine  the  residence  time  of  tropospheric  aoosols.  A  residence  time  of  27 
days  was  found.  The  ratio  exhibits  variations  that  are  probably  seasonal  due  to  seasonal 
changes  in  stratospheric  input  and/or  seasonal  changes  in  mixing  between  the  lower  and 
upper  troposphere. 

The  study  of  ’  ^  P  and  ^  ^  P  in  Cl,  K  and  S  targets  has  allowed  the  determina¬ 
tion  of  the  in  situ  production  rates  of  ’  ^  P  and  ^  ^  P  in  the  upper  ocean  due  to  neutrons  at 
all  latitudes.  The  production  rates  were  estimated  to  be  2.5  10*^  atom/cm^  /min  for  ^  ^  P 
and  8  10  ^  atom/cm^  /min  for  ^  ^  P  at  Bermuda.  In  situ  production  was  found  to  account 
for  5%  for  ^  ^  P  and  1%  for  ^  ^  P  of  the  wet  deposition  rates. 

The  combined  study  of  the  sources  of  ^  ^  P  and  ’  ’  P  to  the  ocean  and  the  dis¬ 
tribution  of  ^  ^  P  and  ^  ^  P  in  suspended  particles  and  plankton  allowed  the  determination 
of  die  turnover  time  of  P  in  zooplankton.  Residence  times  of  P  of  40  to  60  days  were 
found  in  zooplankton,  based  on  the  activity  ratio  ^  ^  P/^  ^  P  in  rainwater,  suspended  matter 
and  plankton. 

Grazing  rates  of  0.024  d  ^  were  determined,  and  it  was  estimated  that  30  ±  10 
mg  C/m^  /d  leaves  the  euphoric  zone  at  BATS  as  fecal  matter  or  pellets.  The  estimated 
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carbcMi  flux  suggests  that  fecal  pellets  and  matter  make  up  the  bulk  of  the  particulate  flux 
at  BATS. 

Inventoiy  calculations  of  ^  ^  P  and  ^  ^  P  in  the  euphotic  zone  suggest  that  nx>st 
oi  the  activities  of  ^  ^  P  and  ^  ’  P  delivered  by  rain  are  found  in  the  suspended  matter.  It 
suggests  that  DOP,  which  is  the  major  pool  of  P  at  BATS,  is  turning  ova*  slowly  relative 
to  the  half-lives  of  ’  ^  P  and  ’  *  P. 

n  FUTURE  WORK 

The  study  of  ^  ^  P  and  ^  ^  P  in  rainwater  has  indicated  the  complexities  in¬ 
volved  in  the  inteipretation  of  the  ratio  in  rainwater.  The  study  has  shown  the  potential 
usefulness  of  cosmogenk  ^  ^  P  and  ^  ^  P  and  investigations  could  be  done  in  the  following 
areas: 


1-  The  tropospheric  production  rates  of  ^  ^  P  and  ’  ’  P  need  to  be  better  known  in  order  to 
determine  the  residence  times  more  accurately.  It  is  especially  important  in  the  case  of 
cosmogenic  ’  ^  P  and  ^  ^  P,  since  the  expected  range  of  ratios  is  small  (Ro  to  Ro  *1.77) 
due  to  the  similarity  in  the  half-lives. 

2-  The  vertical  proflles  of  ’  ^  P  and  ^  ’  P  in  sub-micrometer  aerosols  in  the  troposphere 
shmild  be  determined.  One  of  the  assumptions  of  the  model  of  residence  time  is  that  the 
troposphere  is  well  mixed.  There  is  evidence  that  this  assumption  is  incorrect  and  the 
discord  should  be  assessed. 

3-  Another  way  of  investigating  the  extent  of  mixing  in  the  troposphere,  particularly  with 
respect  to  mixing  between  the  stratosphere  and  the  troposphere,  is  to  study  ’’  Be,  ^  ^  P  and 
^  ’  P  in  rainwater.  The  observed  ratios  ’  Be/’  ’  P  and  ”  P/”  P  can  possibly  be  used  as 
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tracers  oi  mixing  by  comparistxi  with  Ac  expected  ratios  ’  Be/*  *  P  and  ’  ^  P/’  *  P,  assum¬ 
ing  no  mixing. 


4-  The  efficiency  of  Ae  removal  of  sub-micrometer  aerosols  during  waAout  needs  to  be 
determined.  The  model  of  residence  time  relied  on  Ae  assumption  Aat  Ae  efficiency  is 
100%  and  Aat  needs  to  be  assessed. 

The  study  of  ^  ^  P  and  ^  ^  P  m  marine  particulate  matter  has  mAcated  Ae  usefuAess  of 
Aese  tracers  for  studying  zooplankton  dynamics.  Avestigations  could  be  pursued  A  Ac 
followAg  areas. 

1-  DepA  profiles  of  ^  ^  P  and  ^  ^  P  A  DIP,  DOP  and  suspended  particles  need  to  be 
measured  to  assess  correctly  Ac  inventories  of  ^  ^  P  and  ^  ^  P  in  Ae  mixed  layer  and  in  Ac 
euphotic  zone. 

2-  One  of  Ac  most  important  questions  wiA  respect  to  phosphorus  cycling  in  Ae  upper 
ocean  concerns  Ae  availability  of  DOP  to  phytoplankton.  The  ratio  ^  ^  P/^  ^  P  in  DOP 
can  be  used  as  an  AAcator  of  Ae  lability  of  DOP  by  proviAng  a  means  of  estimating  Ae 
residence  time  of  P  in  DOP. 

3-  An  important  question  in  Ae  cycle  of  P  is  Ae  extent  of  remineralization  of  P  in  sink- 

33  32 

Ag  particles  below  Ae  euphotic  zone.  The  ratio  P/  P  in  sinking  particles  at  Afferent 
depAs  can  provide  estimates  of  Ae  regeneration  rates  of  P  from  sinking  particles. 
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Appendix  A 


Table  A.1.  Tool  net  count  rates  (N )  and  count  rates  with  absorber  (N  )  for  the  rain  samples  collected  at 
Woods  Hole  (fust  three  rows)  and  at  sf^uda  during  1990. 1991  and  19^. 


Sample 

date 

Total  net  count  rate 
of^^Pand^’p 

N 

t 

Net  couiu  rate 
of^*Pand’’P 
with  absorber 

N 

a 

1 

a 

cpm 

a 

1/d 

X 

d 

b  I 

incpm 

27-DCC-90 

13.34±]0.07 

0.0430+0.0005 

16.1 

6.42±0.03 

16-Jan-91 

4.33±0.041 

0.0424±0.0008 

16.3 

2.12±0.02 

21-Apr<91 

4.44i0.07 

0.041±0.001 

16.9 

2.33±0.06 

S-Mar-91 

3.49i:|0.09 

0.0402±0.002 

17.3 

1.95±0.05 

lO-Mar-91 

13.74±P.l 

0.044 1±0.0008 

15.7 

6.55±0.06 

lO-Mar-91 

16.97±q.2 

0.043±0.001 

16.1 

8.30±0.20 

20-Apr-91 

13.24±.0.1 

0.041  li:0.0009 

16.9 

6.59±0.09 

20>Apr<91 

20.06±02 

0.0415+0.0008 

16.7 

9.59±p.06 

22-Apr<91 

9.27±0.1 

0.044±0.001 

15.7 

4.32±0.07 

22-Apr-91 

3.51±0.1 

0.045±0.001 

15.4 

1.78±0.03 

1-May>91 

8.60!±0.3 

0.042±0.003 

163 

4.32±0.06 

ll-jun'91 

2.12±0.07 

0.046±0.003 

15.1 

123±0.02 

ll-Jun-91 

1.78±0.08 

0.044+0.005 

15.8 

0.91±0.02 

20-Jun-91 

5.75±0.06 

0.044+0.001 

15.8 

2.69±0.03 

27-Jun-91 

2.66±0.06 

0.043±0.002 

16.1 

134±0.01 

7-Jul*91 

6J3+02 

0.040+0.001 

173 

3.01±0.06 

8-Jul-91 

lJS9±p.05 

0.043±0.002 

16.1 

0.80±0.06 

18-JuI'91 

7.15±0.1 

0.0416±.0.0008 

16.7 

3.16±U.06 

19-Jul-91 

6.95±p.l 

0.042±0.001 

163 

337±0.07 

7-Aug-91 

6.38±0.1 

0.040+0.001 

173 

3.06±0.04 

26-Aug-91 

9.13±0.09 

0.041±0.0007 

16.9 

4.18±0.03 

6-Sq>-91 

4.6SdM0^ 

0.041±0.009 

16.9 

234i0.04 

ll-Sep-91 

5.86±0.07 

0.042±0.0009 

163 

237±0.03 

14-Sep.91* 

0.77±0.04 

0.044 

15.8 

032±0.02 

16-Sep-91 

2.18±p.04 

0.041±0.002 

16.9 

1.08±0.03 

29-Sep-91 

5J7±P.05 

0.042±0.001 

163 

239±0.05 

28-Oct-91 

1032+0.2 

0.044+0.001 

15.8 

5.63±0.09 

28-Oct-91 

333±0.1 

0.042±0.006 

16.5 

1.98±D.05 

7-NOV-91 

5.56i0.1 

0.047±0.002 

14.7 

2.53±0.04 

17-NOV.91 

12.01±03 

0.040±0.002 

17.3 

6.02±0.09 

16-Jan<92 

1.9Q±0.04 

0.042±0.002 

16.5 

0.94±0.02 

21-Jan<92 

17.44±0.08 

0.0416±0.0005 

16.7 

7.71±0.04 

2S-Jan-92 

6.63±0.1 

0.040±0.002 

173 

3.28±0.06 
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31-Jan-92 

8.36±0.2 

0.039±0.004 

17.8 

4.03±0.04 

2-Feb-92 

5.05i0.07 

0.046±0.002 

15.1 

23Q±0.02 

14-Feb-92 

6.95±0.10 

0.042±0.001 

163 

3.18±0.06 

17-Mar-92* 

U7±D.06 

0.044 

1S.8 

0.68±0.02 

26-Mar-92 

7.53±0.2 

0.042±0.002 

163 

3.47±0.18 

26-Mar-92* 

3.57±0.06 

0.044 

15.8 

1.79±0.04 

The  net  count  cue  with  absorber  Na  and  without  absorber  Nt  are  fitted  by  the  following  relationships: 

N  ■a*exp(-a*t)andN  -b^expC-ua^t/M^) 

X«ln2/a.  tisihe  apparent  haft-life  of  the  total  activity  of  ^  ^  P  and  ^  ^  P.  The  parameters  a,  a  and  b 
resulting  fitom  the  least-squares  fit  are  given  with  the  standard  errors. 

*  The  fit  of  the  total  net  count  rate  gave  a  slope  cotie^ionding  to  x  <  143  days.  A  slope  of  0.044  was 
imposed  instead  for  the  fit 
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Table  A  Net  count  rates  of  the  source  with  abs(»ber  (N  }  and  without  (N  )  as  a  function  of  time  for 
the  particle  samples  collected  at  the  Bermuda  Atlantic  Time-^es  (BATS)  siat^  in  1991  and  1992. 


Sample 

Net  count  rate  Net  count 
of”pand”p 

N 

t 

first 

rate  with 
absorber 

N 

a 

day  of 
counting 

a 

a 

b 

L 

0.57±0.01 

0.042±0.003 

0.36±0.03 

3/9/91 

Z67-1R 

1.25±0.03 

0.024±0.002 

0.45±0.04 

3/11/91 

March 

Z67-2R 

0.77±0.04 

0.039±0.004 

O.S6±iO.OS 

3/11/91 

1991 

Z67-3R 

1.18±0.03 

0.017±0.001 

036±0.04 

3/10/91 

Z67-4 

0.61±0.03 

0.040±0.005 

0.46±0.02 

4/6/91 

Z300-2 

0.99±0.04 

0.035±0.005 

05S±fim 

5/5/91 

May 

Z300-3 

1.73±0.04 

0.049±0.004 

1.02±0.02 

5/5/91 

1991 

Z300-4 

1.62±0.03 

0.048±0.004 

1.00±0.02 

5/5/91 

Ml 

0.43±0.02 

0.06Q±0.004 

0.30±0.03 

7/23/91 

Z67-10 

1J0±0.07 

0.036±0.005 

0.88±0.04 

7/23/91 

July 

Z67-11 

1.50±0.04 

0.043±0.003 

0.87±0.03 

7/23/91 

1991 

Z67-12 

0.71±0.03 

0.054±0.005 

034±0.02 

7/23/91 

Sg-1 

l.%±0.1 

0.032+0.004 

1.12+0.06 

7/23/91 

Sg-2 

2.58±0.04 

0.043±0.001 

1.22±0.03 

7/24/91 

M-Sl 

1.01±0.04 

0.047±0.005 

0.53±0.01 

11/19/91 

M-S2 

1.0Q±0.07 

0.047±0.005 

0.57±0.10 

11/18/91 

Z67-10 

1,89±0.05 

0.042±0.003 

1.00±0.03 

11/19/91 

267-20 

3.12±0.2 

0.042+0.02 

1.63±0.03 

11/18/91 

Nov. 

267-21 

3.58±0.11 

0.048±0.003 

1.87±0.05 

11/21/91 

1991 

267-30 

3.64±0.03 

0.050±0.002 

2.14±0.16 

11/18/91 

267-31 

2.10±0.06 

0.042±0.003 

1.03±0.02 

11/25/91 

2300-10 

3.95±0.06 

0.040±0.001 

1.94±0.03 

11/21/91 

2300-11 

3.45±0.09 

0.046±0.002 

1.45±0.02 

11/19/91 

2300-20 

5.01±0.10 

0.046±0.002 

2.42±0.04 

11/20/91 

2300-21 

4.19±0.11 

0.044+0.0(» 

2.11±0.03 

11/20/91 
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M2 

OJ7±0.O4 

0.045±0.005 

035±p.03 

3/2/92 

ZSOO-R 

0.44i:0.02 

0.074±0.005 

0.19+0.02 

3/1/92 

ZSOO-Bl 

1.06±0.0S 

0.055±0.005 

032±0.01 

2/28/92 

Feb. 

ZS00-B2 

0.99±0.04 

0.059i:0.004 

0.43±0.03 

2/29/92 

1992 

Z500-T 

OJ6±0.03 

0.050±0.005 

039±0.02 

2/28/92 

Z67-1 

OJ8i!0.05 

0.084^:0.01 

0.19±0.02 

2/28/92 

Z67-2 

0.27±0.03 

0.10±0.05 

0.06±0.02 

2/28/92 

M-0 

233±0.05 

0.039±0.001 

139±0.01 

5/1/92 

M-3S 

038±0.02 

0.048±0.004 

0.30±0.01 

5/11/92 

M  70 

0.71±0.02 

0.051±0.002 

0.42±0.02 

5/1/92 

April 

ZlSO-Nl 

1.42±0.04 

0.044±0.002 

0.93±0.04 

5/2/92 

1992 

ZlSO-Dl 

2.02i0.04 

0.043±0.001 

1.16±0.02 

5/3/92 

Z150-L 

133±0.03 

0.040±0.002 

0.77±0.07 

5/4/92 

Z150-R 

034±0.02 

0.028±0.001 

038±0.01 

5/1/92 

N  (t)  and  N  (t)  are  flaed  with  exponential  functions  using  a  least-squares  method.  The  equations  used 
fa*  the  fit  are  the  following: 

N  ■  a  •  exp(- a*t) 

N*»b*exp(- In  2*1/14.28) 

The  parameters  a,  a  and  b  and  their  standard  errors  are  given. 
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Table  A3.  Thickness  (d),  yield  (Y)<  die  counting  efficiencies  ^  e),  volume  of  nin  water  (V)  and 
net  count  rates  at  the  time  of  collection  ^  N  and  ^  ^  N)  of  rain  samples  collected  at  Woods  Hole  (first 
3  rows)  and  at  Bermuda.  The  counting  efficiency  of  ^  ^  Pis49±2%forall  thicknesses  encountered. 


27-Dec-90 

16-Jan-91 

21-Apr-91 

4.SQ!tO.OS 

2.50!±0.03 

1.29±0.0I 

5()±2 

30±1 

14±1 

18.9i:0.8 

21.7±P.9 

23.6±0.9 

103C10.05 

3.7210.03 

3.7810.10 

4.8110.04 

2.0510.04 

1.4810.07 

5211 

6211 

6211 

05-Mar-91 

1.90±P.0S 

39±2 

22.810.9 

2.9110.06 

13310.06 

3712 

lO-Mar-91 

1.89±0.02 

20±1 

22.6±0.9 

11.6510.15 

5.9310.11 

6212 

lO-Mar-91 

l.lfi±.0.01 

14±1 

23.8±1.0 

14.70i0.32 

6.2410.06 

6212 

20-Mar-91 

3.4Q±p.03 

38±2 

20.4±0.8 

1 1.710.23 

5.2710.12 

6011 

20-Mar<91 

3.49±|0.03 

39±2 

20.2±0.8 

17.79l0.14 

7.9710.14 

6011 

22-Apr-91 

2.13±0.02 

24±l 

22.3±0.9 

7.2910.09 

3.4510.09 

6111 

22-Apr-91 

1.44±0.01 

16±1 

23.5±1.0 

2.6010.03 

1.52±0.07 

6011 

Ol-May-91 

4.22±0.04 

41±2 

19.3±0.8 

8.0410.11 

3.0610.15 

5311 

lMun-91 

2.05±0.02 

19±1 

22.4±0.9 

3.3810.06 

1.4710.06 

6012 

ll-Jun-91 

221±0.02 

21±1 

22.Q10.9 

23910.06 

0.8310.06 

6112 

20-Jun*91 

4.84±0.05 

54±2 

18.510.8 

6.1710.06 

23710.05 

4011 

27-Jun-91 

3.47±q.03 

40±2 

20.010.9 

5.9010.04 

2.4610.06 

6211 

07-Jul-91 

4.80±p.05 

54±2 

18.510.8 

22.91035 

6.1010.32 

6211 

08-Jul-91 

2J3±0.03 

25±1 

21.710.9 

9.9510.07 

2.910.13 

6211 

18-Jul-91 

5.45±0.05 

61±2 

17.810.8 

13.01031 

4.91032 

6211 

19-Jul-91 

4.87i:0.05 

55±2 

18.410.8 

6.4910.11 

23210.09 

4011 

(r7-Aug-91 

4.20±0.04 

47±2 

19310.9 

10.91037 

3.510.13 

4511 

26-Aug*91 

4.80±q.0S 

54±2 

18310.9 

12.57i0.08 

4.6210.08 

6112 

OS-Sep-91 

2.22ct0.02 

24±1 

22.110.9 

7.4810.13 

23110.12 

6112 

ll-Sep-9l 

331±0.03 

36±2 

20310.9 

7.1210.06 

2.6210.05 

6012 

14-Scp-91 

5.24±0.05 

59±2 

18.010.8 

4.7410.05 

1.7410.09 

6212 

16-Sep-91 

4.53±0.05 

51±2 

18.910.9 

133510,06 

4.1910.07 

6011 

29-Se{>-91 

4.95±0.05 

56±2 

18.410.8 

6.4210.14 

2.4810.07 

5011 

28-Oct-91 

4.76±0.05 

53±2 

18.610.8 

12.8l0.23 

3.810.12 

6211 
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28-Oa-91 

4J5±0.05 

51±? 

18J±0.8 

3.9±P.12 

1.2510.09 

6Qll 

e/7-Nov.91 

4.84±0.0S 

54±? 

10.2110.2 

3310.1 

4911 

17.NOV-91 

4.95±0.05 

55±2 

18.4±0^ 

14.2S±0-2 

4.910.2 

6011 

16-Jan-92 

2.38±j0.03 

29±1 

21.6±0.9 

in±fi.0S 

1.2210.05 

6212 

21-Jan-92 

S.0Q±0.QS 

S6±2 

18J±P.8 

21.110.2 

8.9Q10.1 

6112 

2S-Jan-92 

6.0Q±p.06 

67±? 

17.1±p.8 

7.0510.1 

2.6610.09 

6Qll 

31Jan-92 

S.62±j0.06 

63±? 

17Ji!0.8 

9.5510.09 

3.7310.1 

4612 

02<Feb-92 

6.0S±p.06 

68±? 

17.1dtPA 

S.1810.0S 

2.1810.0 

6111 

14-Fd>-92 

4.44±0.O4 

49±2 

19.0±0.9 

103±0.2 

3.710.1 

6212 

17-Mar-92 

4.1(bd0.04 

46i::2 

19.4±0.9 

2.1210.06 

0.7910.06 

6212 

26-Mar-92 

4.67±0.O4 

52±2 

18.7±0.8 

8.6110.4 

3.3610.2 

6212 

26-Mar-92 

2.45±0.03 

27±2 

21.7±p.9 

4.05±0.09 

1.6610.06 

4612 

Table  A.4.  Thickness  (d),  yield  (Y),  counting  e^iciencies  ^  e  and  ^  ^  e)  and  net  count  rates  at  the 
time  of  collection  ^  N  and  ^  ^  N)  in  plankton  tow  samples  collected  with  a  67-p^m  plankton  net  in 
March  1991. 


Sample 

dale 

depth 

type 

d 

mg/cm^ 

Y 

% 

% 

: 

% 

cpm 

cpm 

L* 

3/4/91 

Suiface 

tow 

2.86±0.03 

/ 

49 

21.1±1.8 

03710.03 

03710.01 

Z67-R1.«.  3/4/91 

Surface 

low 

3.5910.04 

92 

49 

20.Q±1.8 

038i0.0S 

0.96±0.10 

Z67>R2.i.  3/4/91 

Suiface 

tow 

3.8410.05 

89 

49 

19.811.8 

0.79±0.04 

03210.03 

Z67-R3+  3/4/91 

Surface 

tow 

3.2810.04 

91 

49 

20.511.8 

0.12i0.08 

135±0.13 

Z67 

3/31/91 

Suiface 

tow 

5.9310.06 

95 

49 

17.211.8 

0.63ifi.02 

0.1310.01 

*  Operationally  defined  leCractory  fractions  of  thtee  individual  tows. 

*  Labile  fractions  of  all  three  tows  combined. 


Table  A^.  Thickness  (d).  yield  (Y).  counting  efficiencies  ^  ^  e  and  ^  ^  e)  and  net  count  rates  at  the 
time  ot  collection  ^  N  and  ^  ^  N)  in  pimiktoo  tows  collected  with  a  300>iim  net  on  29  April-1  May 
1991.  The  samples  are  triplicates  firom  the  same  tow. 


Sample  Date 

DqNh 

Type 

d 

Y 

m^m^ 

% 

% 

% 

cpm 

cpm 

Z300-2  5/1/91 

Surface 

night  low 

5.93±p.l 

77 

49 

0.7710.01 

03610.01 

Z300-3  5/1/91 

Surface 

night  low 

13.1±q3 

90 

47 

11.4iP5 

1.43i0J>3 

03410.02 

Z30(M  5/1/91 

Surface 

night  low 

74 

47 

11.0»t0.4 

1.41i0.02 

0.4410.02 

Table  A.6.  Thickness  (d).  yield  (Y),  counting  efficiencies  (^^eand^^e)andnet  count  rates  at  the 
time  of  collection  C^Nand^^N)in  suspended  particles  and  in  plankton  tows  collected  with  a  67-p.m  net 
during  BATS  34. 15-19  July  1991. 


Sample 

dale 

dqxh 

type 

d 

Y 

”N 

mg/cm^ 

% 

% 

% 

qim 

q>m 

Ml 

7/16/91 

surface 

pump 

2.I81O.O2 

48 

49 

22.110.9 

0.4510.04 

0.1510.02 

267-10*  7/18/91 

Surface 

night  tow 

8.2910.08 

87 

48 

14.9103 

12410.06 

03310.05 

267-11*7/18/91 

Surface 

night  tow 

7.3310.07 

82 

48 

15.810.6 

1.0710.04 

0.4410.03 

267-12*7/18/91 

Surface 

night  tow 

3.8710.04 

87 

49 

19.710.8 

0.4810.03 

03510.03 

Sg.l+ 

7/16/91 

Surface 

macroalgae  7.84±0.08 

78 

48 

15.3103 

IOI1O.IO 

0.9910.1 

Sg-2+ 

7/16/91 

Surftce 

macroalgae  6.93i0.07 

87 

48 

16.210.6 

2.3810.07 

1.0710.05 

*  Triplicates  of  a  plankton  tow. 

+  Diq>licates  of  a  sample  of  sargassum. 
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Table  A.7.  Thickness  (d),  mid  (Y).  coundng  efficiencies  ^  e  and  ^  ^  e)  and  net  count  rates  at  the 
time  of  ctdkctioo  ^  N  and  ^  N)  in  plankton  tow  and  suspended  particulate  samples  (3-lim  mesh  size) 
collecied  during  BATS  38, 11-lS  November  1991.  Z67  and  Z300  refer  to  plankton  tows  collected  with  a 
67-iim  and  300-|im  plankton  net  Each  tow  was  split  in  two  except  Z67-10. 


Sample 

dale 

depth 

type 

d 

Y  e 

mg/cm^ 

%  % 

% 

cpm 

cpm 

M-Sl 

11/12/91 

surfKe 

pump 

2.6210.03 

69  49 

214i0.9 

0.7510A2 

0.5310.03 

M-S2 

11/14/91 

surface 

pump 

3.31±0.O3 

68  49 

204±0-8 

0.6310.11 

0.4810.12 

Z67-10 

11/12/91 

surface 

night  tow 

Z91±0.03 

67  49 

2I.I1O.8 

I.S61O.O5 

0.8310.05 

Z67-20* 

11/14/91 

(F60m 

day  tow 

7.76±D.08 

83  48 

IS.410.6 

2.3210.04 

I.2I1O.I4 

Z67-21* 

11/14/91 

0-60m 

day  tow 

7.45±0.07 

90  48 

15.710.6 

3.02t0.09 

lJ6l0i)9 

Z67-30 

11/14/91 

surface 

night  tow 

9.95±0.12 

96  47 

13.5i0^ 

2.9510.2 

I.I81O.12 

Z67-31 

11/14/91 

surface 

night  tow 

5.09±0.05 

66  48 

18.210.7 

l.f>  710.04 

I.I61O.OS 

Z300-10 

11/12/91 

surface 

night  tow 

10.8i^.l3 

95  47 

13.Q10.5 

3.9710.07 

I.6Q1O.OS 

Z300-11 

11A2/91 

surface 

night  tow 

6.47±0.06 

91  48 

I6.61O.7 

2.6910.04 

1.7010.07 

Z300-20*  11/14/91 

0-60m 

day  tow 

8.62i0.09 

9S  48 

I4.61O.S 

3.S61O.OS 

1 .8810.07 

Z300-21*  11/14/91 

0-60m 

day  tow 

8.2410.08 

79  48 

14.910.5 

3.3110.04 

1.8410.07 

*:  tows  collected  by  yoyoing  the  net  on  the  hydrowire  from  0  to  the  base  of  the  mixed-layer  at  60  m. 


Table  A,8.  Huckness  (d),  yield  (Y),  counting  efficiencies  (^^eaiid^^e)andnet  count  rates  at  the 
time  of  ctdkction  ^  N  and  ^  N)  in  plankton  tows  and  suspended  particulate  samples  (Millipore  3  itm) 
samples  collected  during  BATS  41A,  24-2S  February  1992.  ZSOO  and  Z67  refer  to  samples  collected  with 
plankton  nets  of  mesh  size  SOO  fim  and  67  itm  te^)ectively.  One  of  the  ZSOO  was  separated  into  two 
fractions  ZSOO-T  and  ZSOO-B.  ZSOO-T  is  the  fraction  that  stayed  subtended  in  seawater.  21S00-B  is  the 
fraction  that  sank.  ZSOO-Bl  and  ZSOO-B2  ate  duplicates.  ZSOO-R  is  the  refractory  fraction  of  a  tow,  the 
labile  fraction  was  lost  during  processing. 


Sample 

date 

depth  type 

d 

Y 

”N 

mg/cm^ 

% 

% 

% 

cpm 

cpm 

M2 

02/25/92 

surf. 

pump 

3.6410.04 

70 

49 

20.Q10.8 

O3I1O.OS 

0.1910.03 

ZSOO-R 

02/25/92 

surf. 

night  tow 

4.4410.04 

71 

49 

19.Q10.8 

03110.01 

03010.04 

ZSOO-Bl 

02/25/92 

surf. 

night  tow 

9.4010.09 

94 

47 

13.910.5 

O.6I1O.O2 

0.5310.03 

Z500-B2 

02/25/92 

surf. 

night  tow 

IO.QIO.1 

93 

47 

13.510.4 

03110.04 

0.59i0.06 

Z500-T 

02/25/92 

surf. 

night  tow 

5.2Q10.05 

73 

48 

15310.5 

0.3710.02 

0.2410.03 

Z67-1 

02/25/92 

surf. 

night  tow 

4.8210.05 

72 

48 

18310.7 

0.2410.02 

0.3810.07 

Z67-2 

02/25/92 

surf. 

day  tow 

2.5I1O.O3 

52 

49 

21.710.9 

0.0710.02 

03110.08 
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Table  A^.  Thidmess  (d)t  yield  (Y),  counting  efiBctencies  ^  e  and  ”  e)  and  net  count  tales  at  the 
tiine  of  coUeciion  ^  N  and  ^  ^  N)  in  suspended  panicles  and  in  idankton  tows  collected  during  EN  235  in 
April  1992  at  BATS  station.  ZISO  refer  to  a  tow  collected  with  a  ISO-itm  net  ZISO-L  and  R  refer  to  the 
labile  and  refiraciory  ftactioos. 


Sample 

Dale 

depth  type 

d 

Y 

mg/cm^ 

% 

% 

% 

qxn 

cpm 

M-0 

4/2S/92 

SUfl 

pump 

5.7510.06 

63 

48 

17.410.7 

1.81i0.02 

1.O910.O3 

M-35 

4/25/92 

35  m 

pump 

3.5S10O4 

36 

49 

20.110.9 

OJ7iO.ll 

0.46±0.12 

M-70 

4/2S/92 

70  m 

pump 

4AS10.0S 

44 

49 

18.710.7 

0J210.02 

03510.02 

ZlSO-Nl 

4/2d/92 

surf. 

night  tow 

7.0010.07 

59 

48 

16.1±0.7 

1.37±0O6 

O.4O1O.O3 

ZlSO-Dl 

4/26/92 

surf. 

day  tow 

8J710O8 

74 

48 

14.910.5 

1J910.03 

0.7110.03 

ZISO-L 

4/26/92 

suif. 

night  tow 

8.9310.09 

76 

48 

143±0.5 

I.2S1O.12 

O.41±0.O4 

ZISO-R 

4/26/92 

surf. 

night  tow 

6.49±0.06 

68 

48 

I6.61O.7 

0.7810.02 

0.22±0.01 

Table  A.10.  Thickness  (d),  yield  (Y).  counting  efficieocies  ^  £  and  ^  ^  C)  and  net  count  rates  at  the 
time  of  collection  ^  N  and  ^  ^  N)  of  the  samples  exposed  to  cosmic  rays  in  Woods  Hole  (S3  °  N).  Ml 
Wadiington  (55.7  °  N)  and  1'  AiguiUe  du  Midi  (47  °  N). 


Sample 

site 

target 

element 

d 

mg/cm^ 

Y 

% 

% 

% 

cpm 

- 

cpm 

Cl-WH 

Woods  Hole 

a 

3.62 

81 

20 

49 

03110.03 

O.Q34±0.008 

a-MW 

MtWashingtan 

a 

2.18 

50 

22 

49 

0.5510.09 

0.04510.01 

CU-MW  Ml  Washington 

a 

2.05 

48 

22 

49 

0.3610.02 

O.O6Q1O.OI 

K-MW* 

hft.  Washingioo 

K 

9.89 

90 

/ 

47 

O.I61O.O3 

/ 

S-MW* 

Ml  Washington 

S 

11.9 

80 

/ 

47 

O.lSlO.03 

1 

Cl-AM 

Aiguille  du  Midi 

a 

3.93 

78 

19.6 

49 

O.SOlO.03 

0.04910.004 

*  The  activity  of  P  could  not  be  detected  in  that  sample. 
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